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ABSTRACT 


Comparative results of two magnetometer surveys, one in December, 1927, before 
the discovery of oil, and one in January, 1932, after production had reached maturity, 
show that the vertical component of the earth’s magnetic field has been changed in de- 
tail by the presence of production equipment. The generalized magnetic pattern 
remains similar, however, for both surveys. 


The first problem outlined at the initiation of an exploratory pro- 
gram of earth magnetics 5 or 6 years ago was, almost without excep- 
tion, the definition of the vertical component pattern across well 
known and highly developed producing structures. The results of 
such surveys did not, in general, fit into the preconceived pattern as 
anticipated by the geologist-geophysicist. There were, however, such 
noteworthy exceptions as the buried extension of the Wichita Moun- 
tains in the Texas Panhandle and the southern portion of the Granite 
ridge in the vicinity of Eldorado, Kansas. With the knowledge that 
the magnetic pattern over these types of structure was so readily 
reconciled with the configuration of the granite, attempts were made 
at reconciling the known geological features in other areas with the 
measured magnetic anomaly, but these attempts proved, too fre- 
quently, unsuccessful. The geologist-geophysicist faced with the 
necessity for an explanation of phenomena with which he was little 


1 Read before the Association at the Oklahoma City meeting, March 25, 1932. 
* Geologist-geophysicist. 


1171 


/ 
/ 


1172 O. C. CLIFFORD, JR. 


acquainted, but concerning which he had definite convictions, ex- 
pounded such hypotheses as the hysteresis induced by buried casing, 
the possibility of polar reversibility under the impact of cable tools, 
et cetera ad infinitum, being limited only by an inadequate. vocabulary. 


” 


F 


a 


147 


7 


/ 


VERTICAL MAGNETIC INTENSITY 


CONTOUR INTERVEL $0 SCALE OF MILES 
‘ 


OC 


Fic. 1.—Map showing vertical magnetic intensity, 1927. 


During the past 6 years the study of micro-geomagnetics has pro- 
gressed beyond the limits of early credulity and, with the publica- 
tions of Heiland, Wilson, Barret, Lynton, Spraragen and others, 
many of these problems are nearing solution. Though other geo- 
physical methods have rightfully supplanted magnetic methods 
throughout most of the Mid-Continent and Gulf Coast areas, mag- 


\ | 

MN 
SUE ai 

wi 


MAGNETIC RESURVEY OF OKLAHOMA CITY FIELD 1173 


netometers must still be considered of exploratory value for special- 
ized structural types and for rapid, inexpensive surveys of regional 
character. 

There has been no material published, with one exception, on the 
measured change in the vertical component of the magnetic field in 
an area surveyed both before the discovery of production and after 
development had reached maturity. Barret! has published an interest- 
ing summary of such work in the Sligo field, Bossier Parish, Louisiana, 
which must be considered inconclusive because of the small number of 
control points. With the permission of the management of the Indian 
Territory Illuminating Oil Company, the results of two magnetic 
surveys in the vicinity of Oklahoma City are presented. 

In December, 1927, a magnetic survey of the vertical intensity 
anomaly was completed in the vicinity of Oklahoma City, with the 
use of a Schmidt vertical intensity variometer manufactured by the 
American Askania Corporation. The results of this survey contoured 
on a 50-gamma interval are shown in Figure 1. The dominant feature 
is an area of negative relief centered around the west half of Sec. 12, 
T. 11 N., R. 3 W. Here a 110-gamma closure is the maximum ob- 
served. No attempt is made to interpret this feature with respect to 
the existence of Permian, Pennsylvanian, Ordovician, or pre-Cambrian 
structure. The structure was undoubtedly there and the magnetic 
picture appeared as shown. 

The discovery well of the present Oklahoma City field was spudded 
in near the center of the SE. 4 of the SE. 4 of Sec. 24, T. 11 N., 
R. 3 W., June 12, 1928, and was completed December 4, 1928. An in- 
tensive drilling campaign followed and by January 1, 1932, a total of 
867 wells had been completed, of which 5 had been abandoned. This 
represents, according to Charles,’ a development of 80 per cent of all 
acreage which is thought to have possibilities of becoming produc- 
tive. The structural conditions of Oklahoma City have been de- 
scribed by Charles* and others. 

In mid-January, 1932, a resurvey of the vertical component of 
the earth’s magnetic field was made. The results are shown in Figure 
2, also contoured on a 50-gamma interval. 


1 William M. Barret, “Magnetic Disturbance Caused by Buried Casing,” Bull. 
Amer. Assoc. Petrol. Geol., Vol. 15, No. 11 (November, 1931), pp. 1387-88. 


2H. H. Charles, verbal communication. 


3H. H. Charles, “Oklahoma City Oil Field,” Bull. Amer. Assoc. Petrol. Geol., Vol. 
14, No. 12 (December, 1930), pp. 1515-33. 
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There appears, at first glance, to be a marked discrepancy be- 
tween the work of 1927 and that of 1932. This becomes less noticeable 
on closer scrutiny. The dominant feature of the later survey, also, is 
an area of negative relief, but this has been shifted from the west half 
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Fic. 2.—Map showing vertical magnetic intensity, 1932. 


of Sec. 12, T. 11 N., R. 3 W., to the SE. 4 of Sec. 10, T. 11 N., R. 
3 W. The maximum closure appears to be approximately the same, 
though it may be greater in the later work. In the recent survey, the 
contours lie closer together and the general pattern practically dupli- 


cates that of 1927, but the plane of the feature appears to have been 
tilted toward the west. 
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The base point, SE. corner, Sec. 36, T. 11 N., R. 3 W., is identical 
for both surveys and is thought to be in an undisturbed state, as no 
well is closer to it than 34 mile. Between the base and a point two 
miles farther east, also in an undisturbed condition, there is a marked 
accordance of anomaly, a difference of plus 34 gammas being measured 
in the first survey and a difference of plus 31 gammas in the second 
survey. The stations south and west of the base point, though they 
do not coincide, either in position or value, do roughly correspond. 
It appears, then, that in the undisturbed area on the south there is 
a pronounced accordance of anomaly value which gives a reasonable 
check on the comparative accuracy of the two surveys. 

The most pronounced discordance between the two surveys ap- 
pears one mile north of the base station. In the original survey this 
point falls on a gradual northward decline of anomaly which con- 
tinues uninterrupted for 5 miles. In the more recent work, there is a 
sharp decline of intensity in the first mile, followed by an almost 
equally sharp rise in the succeeding mile. Though every effort was 
made to keep pronounced effects of extraneous material from the 
survey, it is thought that this discrepancy represents such effects. 
When the instrument gave readings off the field of view, it was moved 
until large effects appeared vitiated. 

The changed aspect of the vertical component map may be at- 
tributed to one or all of several factors. No attempt is made to evalu- 
ate the relative importance of these. They are: (1) the recent survey 
was more compactly observed and for that reason the intensity 
changes may appear more accentuated than in the original survey; 
(2) the physical development produces major changes in the mag- 
netic field; (3) the actual extension of Oklahoma City proper, both 
laterally and vertically, has produced immeasurable changes in the 
magnetic field; and (4) secular variation causes slight differences. 

The results from one resurvey of this type can not be considered 
conclusive. It is reasonable to expect, however, that similar attempts 
in other areas will produce similar results. The results of the present 
survey tend to show that, with care in the selection of observation 
points, the general picture of an area will show little difference before 
and after the development of production; that, if the original mag- 
netic structure had any considerable relief, it will retain not only 
the same approximate closure, but also closure in the same direction. 
The early magnetic work in existing fields can thus be expected to 
have shown a rough approximation to the pre-development state. 
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It is hoped that similar surveys will be undertaken in areas which 
were magnetically defined prior to the discovery and development of 
production. A resurvey of Hobbs field, New Mexico, and various of 
the Permian basin fields of Texas should show interesting relation- 
ships. 
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ABSTRACT 


The local magnetic vectors at the United States Coast and Geodetic Survey stations 
in Louisiana, Texas, Arkansas, Oklahoma, Kansas, Mississippi, Alabama, and Cali- 
fornia have been-computed by deducting the normal values of the earth’s magnetic 
field from the absolute measurements, at each station, of the declination and of the 
vertical and horizontal magnetic intensities. 

The local magnetic vectors have been plotted as vector triangles at the respective 
stations on the vector maps of the different states. 

These vectors indicate the intensity and the direction in space of the magnetic 
lines of force, as due to local magnetic anomalies mainly within the first 15,000 feet of 
subsurface. 

As local magnetic anomalies are, with negligible exceptions, the result of geologic 
features, either of structural or petrographic character, a large amount of regional 
and local geologic information is obtained by a study of these maps. 

he main magnetic anomalies of the different states have been interpreted in 
terms of geology, though no attempt has been made to exhaustively interpret all 
vectors. 

The vector maps show at a glance which areas are of interest for magnetometer 
surveys, what size anomalies may be expected, how far a magnetometer survey has to 
be extended, or where it is best commenced to cover a certain area. 

As the vectors are of sufficient accuracy, they allow a checking and tying-in of 
scattered magnetometer surveys. 


INTRODUCTION 


Extensive studies of the anomalies of the magnetic vertical in- 
tensity and their correlation with geology were made some years ago 
in Europe,’ in the Mid-Continent and Gulf Coast areas* and in 
California.® 


1 Read in part before the Association at the Oklahoma City meeting, March 25, 
1932. Manuscript received, May 18, 1932. 


2 Geologist and geophysicist, 6241 Richmond Avenue. 


3 A. Nippoldt, “Karten der Verteilung des Erdmagnetismus und seiner oertlichen 
vertikalen Stoerungen in Europa,” Archiv des Erdmagnetismus, Heft 6 (Julius Springer, 
Berlin, 1927). 


‘ Magnetometer survey of seven states in the Mid-Continent and Gulf Coast 
om by L. Spraragen, series of articles published in the Oil and Gas Journal (1928- 
29). 

5 G. B. Somers, “Anomalies of Vertical Intensity Compared with Regional Geology 
for the State of California,” Colorado School of Mines Magazine (September, 1929). 
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The investigations were based on the measurements made by the 
surveys of the respective European states, and by the United States 
Coast and Geodetic Survey.' 

Though by these studies a definite connection between broad 
structural features—such as the West Texas Permian basins, the 
Central Mineral region, the main geological trends in California and 
Europe—and magnetic anomalies of the vertical intensity could be 
definitely established, such studies could only incidentally yield any 
information about the type of structure which is of special interest 
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to the commercial geologist, because the distances between the mag- 
netic stations, ranging from 20 to 50 miles, were far too large. 

It is quite evident that two, three, and more magnetic “highs” and 
“lows” may occur between stations placed at such distances and the 
lines of equal magnetic intensity, drawn on the base of such stations, 
as a rule veil, rather than set forth, the local magnetic anomalies 
which are due to small “structures.” 

Study of Figure 1 will make the previous statement clear. This 
cross section through two anticlines separated by a syncline demon- 
strates schematically the paths of the magnetic lines of force and the 


1 “United States Magnetic Tables and Magnetic Charts for 1925 by U. S. Depart- 
ment of Commerce,” U.S. Coast and Geodetic Survey Serial 453. 
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local magnetic vectors as due to the local “structures” only, after the 
normal magnetic field of the earth has been eliminated. The total 
local vector is shown to be resolvable into two components—the local 
vertical and the local horizontal magnetic vector. 

We have assumed that the material of the anticlines is of more 
than average magnetic permeability and the material of the syncline 
of less. Thus the magnetic “high” corresponds with a structural 
“high,” and the magnetic “low” with a structural “low.” This is the 
ordinary condition, for example, if two granite ridges are embedded 
in shales and limestones. Should the material of the anticlines be of 
less than average magnetic permeability and the material of the 
syncline of more, then the shape of the magnetic lines of force remains 
the same, but their direction is opposite to the direction shown in the 
picture. Thus a magnetic “low” would correspond with a structural 
“high” and vice versa. This condition exists in many places, for ex- 
ample in Kansas, where granite plugs or ridges pierce through the 
overlying highly magnetic schists. 

For a clearer understanding, let it be assumed that the magnetic 
effects of the two anticlines and the syncline are equivalent to the 
effects of two negative (—P) poles, and one positive (+P) pole, as 
shown in Figure 1. 

Above the negative poles the magnetic intensity is directed per- 
pendicularly downward toward the poles (positive vector). Above the 
positive pole the magnetic intensity is directed perpendicularly up- 
ward, away from the pole (negative vector). Along the surface, half 
way between the positive and negative poles, the two vertical tenden- 
cies compensate each other and only a horizontal component remains, 
directed from the positive toward the negative pole (neutral vector). 
Between the vertical and horizontal directions of the intensity there 
occurs a gradual change, as shown in the picture. As long as the vector 
is still directed downward, it may be called a positive vector; if the 
direction is upward, a negative vector. 

It might happen now that stations I and III had been occupied 
by the Coast and Geodetic Survey. As both stations have the same 
magnetic vertical intensity they lie on the same isogam, which would 
pass over the entire structure without taking the slightest notice of it. 
The same would happen, if also stations III and IX had been sur- 
veyed. Any other combination of the vertical intensity of two or 
three stations give widely varying indications of some anomaly. We 
should consider, however, that most of the “structures,”’ which create 
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the local anomalies, lie at depths ranging from 2,000 to 15,900 feet 
and that the magnetic effect is felt within a horizontal distance from 
the edge of a “‘structure” of only about twice its depth. Therefore it 
is reasonable to assume that, within a distance of 20-50 miles, only 
one station is occupied, as a rule, on a particular “structure” and there 
is little use in combining the measurements of the different stations, 
unless broad structural features are suspected. 

Contrary to the opinion of some magneticians,'? the writer holds 
that at such large distances between stations, much more detailed 
information may be obtained by studying the horizontal intensity 
than by studying the vertical intensity, and thinks that this results 
clearly from a study of the horizontal intensities in Figure 1, especially 
if it be realized that the anomalies are 3-dimensional. 

The best information is of course obtained by a combination of the 
horizontal and the vertical intensities of a magnetic vector in space. 

The Coast and Geodetic Survey has measured the declination and 
the horizontal and vertical intensities at the occupied stations. If 
the normal declination and the horizontal and vertical intensity be 
figured for each station and these normal values are deducted* 
from the actually measured values, the local deviation of the declina- 
tion and the local vertical and south-north horizontal component of 
the total local magnetic vector are obtained as the difference. From 
the deviation of the declination may be easily figured also the local 
west-east component of the total local magnetic vector. 

The average absolute accuracy of the components of these local 
magnetic vectors empirically lies at +40 gammas, the average ac- 
curacy relative to neighboring stations at + 20 gammas. 

The local magnetic vectors have been plotted at their respective 
stations on the maps of the principal oil-producing states, in the 
shape of vector triangles. The two dashed lines represent the hori- 
zontal and the vertical vectors; the hypotenuse is the total vector. 
If the vector triangle be turned through go° around the horizontal 
vector, which is of course the dashed line beginning at the station 
point, the direction in space of the total vector may easily be visual- 
ized and the direction marked by a pin, if the map is mounted on 
cardboard. 

1 A. Nippoldt, op. cit. 

2G. B. Somers, op. cit. 

3 A. Nippoldt, of. cit. 


4 L. A. Bauer, “Chief Results of a Preliminary Analysis of the Earth’s Magnetic 
Field for 1922,” Terr. Mag., Vol. 28 (1923), pp. 1-28. 
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In studying these local vector maps, it should be considered that, 
as shown in Figure 1, there exists somewhere along the prolongation 
of the horizontal vector either a magnetic “high” or “low,” which 
corresponds with some kind of geologic feature, either of structural 
or petrographic character. 

The great and irregular variations of the total intensity and of the 
directions of most of the horizontal vectors should make it clear that 
most of these features necessarily must be localized. 

This is the reason why no attempt has been made to connect the 
stations by certain curves, indicating, for example, equal total in- 
tensities, equal angles of inclination or equal vertical and horizontal 
intensities, though by help of the vectors the last two sets of curves 
could, of course, be much more intelligently computed than from the 
vertical or horizontal components alone. 

A few of the regional geologic “structures” can be quantitatively 
interpreted by the help of these vectors. 

Due to the long distances between stations, it has been necessary 
to use the most reasonable geological interpretation of the vector 
maps. Though two, three, or more local “highs” and “lows” are 
possible between two neighboring stations, it is probable that re- 
gional geologic “structures” extend at right angles to the horizontal 
magnetic components, where the majority of the latter are nearly 
parallel with one another over a large area, that is, a common geologic 
feature is probably responsible for a series of vectors, if they are all 
directed toward, or away from, the same point. 

Due to the long distances between stations, the interpretation 
was mainly restricted to regional magnetic anomalies, though it is 
expected it will be found that the main value of the vector maps lies 
just in these areas, where the magnetic anomalies are so localized, or 
so shallow that no relation between the individual vectors can be 
established. 

For all those areas where the vectors indicate regional magnetic 
tendencies, we may reasonably assume that the vectors are mainly 
due to deeply buried regional structure. But if the main magnetic in- 
fluence in these areas originates at great depths and if the total inten- 
sity of the vectors amounts to only a few hundred gammas, we must 
conclude that the average magnetometer surveys can not possibly 
yield much information with regard to commercial structure, unless 
the deeper structure is truly reflected within the shallower beds. The 
results of the average magnetometer surveys are comparable with the 
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results which we would obtain with torsion balances, if their accuracy 
were reduced to +10 Eétvés. If, however, the accuracy of the mag- 
netometer surveys is greatly increased, these surveys may yield im- 
portant information about structure at commercial depths also in 
these areas where the magnetic anomalies as located by the average 
magnetometer survey are unquestionably due to deeply buried re- 
gional structure.' 

As the vectors are of sufficient accuracy, they allow a checking and 
tying-in of scattered magnetometer surveys. These stations should 
prove of the same value to 3-dimensional magnetometer surveys as 
pendulum stations to gravimetric surveys. 

The vector maps show at a glance which areas are of interest for 
magnetometer surveys, what size anomalies may be expected, how 
far a magnetometer survey has to be extended or where it is best 
commenced to cover a certain anomaly. 

In the following study of the vector maps no attempt has been 
made for an exhaustive interpretation of all vectors. The reader may 
easily perceive an abundance of further detailed information, which 
may interest him, especially in connection with scattered magnetom- 
eter surveys and detailed geologic information. 


LOUISIANA” 


The vector map of this state (Fig. 2) gives a true picture of the 
known main structural trends and further indicates many distinct 
“structures.” 


Sabine uplift—The horizontal components of the two negative 
vectors at Benton and Shreveport concur at a point about 1o miles 
southwest of Shreveport, at which point therefore a magnetic “high” 
is suggested. A line connecting this “high” with the “high” indicated 
a few miles south of Many coincides well with the axis of the Sabine 
uplift. 

North Louisiana geosyncline—The north Louisiana geosyncline 
is expressed by the negative vectors at Ruston, Jonesboro, and 
Winnfield. 


1D. M. Collingwood, “Magnetics and Geology of Yoast Field, Bastrop County, 
Texas,” Bull. Amer. Assoc. Petrol. Geol., Vol. 14, No. 9 (September, 1930), pp. 1191-97. 

E. D. Lynton, “Some Results of Magnetometer Surveys in California,” Bull. Amer. 
Assoc. Petrol. Geol., Vol. 15, No. 11 (November, 1931), pp. 1351-70. 

A. Nippoldt, Verwertung magnetischer Messungen zur Mutung (Berlin, 1930), p. 58. 


2 The local magnetic vectors are based on the information given in ‘United States 
Magnetic Tables and Magnetic Charts for 1925,” op. cit. 
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Area east of geosyncline—East of the geosyncline are found two 
trends of magnetic “highs,” which parallel the axis of the Sabine 
uplift: (1) the Farmerville-Columbia-Harrisonburg trend and (2) the 
Bastrop- Winnsboro trend. 

On the second trend lie the gas fields of Monroe and Richland. 
The magnetic “high” indicated slightly north of Bastrop is interpreted 
as an uplift in the deep basement rocks, which uplift accounts for the 
Monroe gas field. From magnetometer work we know that half way 
between Rayville and Winnsboro a similar magnetic “‘high” and cor- 
responding uplift accounts for the Richland gas field. 

The magnetic “high” as indicated south of Winnsboro therefore 
suggests other possibilities of uplifts along the known trend toward 
the south. 

Central Louisiana.—A magnetic “high” south of Colfax is indi- 
cated. 

The concurrence of the horizontal vectors of Jena, Alexandria, 
and Marksville at a point about 10 miles northeast of Marksville 
suggests a magnetic and possibly structural “high” of large dimen- 
sions at that point. 

Other magnetic “highs” are suggested between Leesville and De 
Ridder and southeast of Cheneyville. 

At St. Francisville and New Roads there is a striking example of 
two stations which may reasonably be assumed to lie on the same 
“structure.” The magnetic lines of force come out of the ground at a 
low angle at St. Francisville, assume a horizontal direction between 
the two stations, and penetrate into the ground at a low angle at 
New Roads, directed toward a magnetically positive “‘structure” on 
the south. 

The distances are too great to say that this “structure” is the 
eastern extension of a west-east ridge or anticline as indicated be- 
tween Ville Platte and Opelousas, but this possibility is strongly 
suggested. 

Southern Louisiana.—There is a northwesterly trend of magnetic 
“highs” from Morgan City to St. Martinsville; this trend parallels 
a line through the Five Island domes, about 20 miles west, and may 
correspond with the so-called Iberian structural axis. 

No continuation of this trend beyond St. Martinsville can be 
preceived from the vectors. An east-west trend of magnetic “highs” 
from St. Martinsville to Lake Charles, however, is clearly noticeable. 
This trend passes south of Lafayette, northeast of Crowley, north 
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of Evangeline, south of Lake Charles, and corresponds with the well 
known structural axis in southern Louisiana along which, besides 
several unproved prospects, lie the domes and oil fields of Jennings, 
Roanoke, Welsh, and Iowa. The indicated “high” south of Lake 
Charles suggests the conclusion which has been drawn from recent 
seismograph work in this area, that the structural axis turns south- 
ward past the Iowa dome and continues south of Lake Charles 
toward the Lockport dome. 

The large negative vectors at DeRidder, Oakdale, and Marble 
Quarry and the horizontal vector at Oberlin suggest a west-east 
magnetic “‘high” trend south of Oberlin. Due to the large intensity of 
these vectors, it does not seem justifiable to connect them with the 
magnetic “‘high” trend of Lafayette-Lake Charles but it may be con- 
cluded that the magnetic “high” or anticlinal structure indicated be- 
tween Ville Platte and Opelousas extends westerly past the Texas 
border. A corresponding gravimetric high trend is indicated by 
torsion-balance work and is usually interpreted as due to a deep 
petrographic or structural regional feature. 


TEXAS! 
REGIONAL STRUCTURE 


Delaware basin—The magnetic center of the Delaware basin 
seems to lie between Pecos and Barstow. The negative vector at 
Pecos is directed toward the west, and at Barstow, toward the east 
or toward the Central Basin platform, which is indicated by the 
positive vector at Monahans. 

Main Permian basin.—In the southwest part of the main Permian 
basin the negative vectors from Odessa to Big Springs are directed 
toward the southeast, and the negative vectors at Andrews and 
Lamesa toward the northwest, which means that the magnetic axis 
of a basin lies between these two sets of stations. 

The Bend arch.—The Bend arch is indicated by the positive vec- 
tors at Brownwood and neighboring stations, with a north extension 
to Jacksboro. 

Central Mineral region —The horizontal components of the posi- 
tive vectors at Llano and Burnet and the neutral vector at Lampasas 
all concur at a point half-way between these three stations. This 


1 The local magnetic vectors are based on the information given in “United States 
Magnetic Tables and Magnetic Charts for 1925,” op. cit. 
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point should correspond with the center of a large magnetic “‘high,” 
due to the Llano-Burnet uplift. 

East Texas syncline—The East Texas syncline is indicated by 
many negative vectors in the northeastern part of the state. 

Northwest Texas.—A basin of considerable extent is indicated in 
northwest Texas, with its magnetic center half-way between Here- 
ford, Dimmitt, and Bovina. The rim of this basin probably extends 
from Bovina northward toward Endee, eastward and southward 
toward Canyon, Tulia, and Plainview, and westward north of Olton 
and Muleshoe back to Bovina. This basin may hold oil possibilities 
similar to the Delaware or Black Shale basins. 

Gulf Coast—Many of the horizontal magnetic components along 
the Gulf Coast are nearly perpendicular to the coast; we may there- 
fore expect regional “structures” parallel with the coast. 

Especially in southwest Texas such “structures” seem to be 
clearly indicated. We may further follow up the regional magnetic 
high trend between Ville Platte and Opelousas, as discussed in the 
last paragraph on “Louisiana,” south of Newton, Woodville, Cold 
Spring, Conroe, Hempstead, Bellville, Columbus. This trend may 
prove of interest in the light of recent developments south of Conroe. 


LOCAL STRUCTURE 


Anticlinal structure at Big Lake is indicated by a positive vector 
at Rankin and at Stiles. 

A southwest-northeast ridge is indicated by the positive vectors 
at Sherwood and San Angelo. 

A magnetic “high” is indicated between Ozona, Eldorado, and 
Sonora. With this “high” may correspond one “structure” or more, 
which, in the light of recent development west of Ozona, may prove 
of great interest. 

The great variety in the magnetic field along the coast should 
prove that magnetometer surveys there may be not only of a regional, 
but also of a local, structural interest in many places. 

The Refugio anticline and the Pierce-Junction and the Liberty 
domes are indicated by positive vectors. 

It is possible that these magnetic “‘highs” are the result of the 
uplifted strata at relatively shallow depths, though the possibility 
also exists that they are due to uplifts in the basement. In any case, 
more detailed magnetic investigations may prove of great value, es- 
pecially in southwest Texas. 


| 
| 
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ARKANSAS! 


Arkansas has been divided into four main geological sections: (1) 
the Gulf Coastal Plain, south and east of a line passing through 
DeQueen, Arkadelphia, Little Rock, and Pocahontas; (2) the Arkan- 
sas Valley, between a line from Van Buren to Heber Springs on the 
north and a line from Waldron to Little Rock on the south; (3) the 
Ouachita Mountain region, between the Arkansas Valley and the 
Gulf Coastal Plain; (4) and the Ozark region, north of the Arkansas 
Valley. 

Gulf Coastal Plain.—The southern part of the Gulf Coastal Plain 
from DeQueen to Lake Village is indicated as a large negative anom- 
aly, interrupted by a north-south trend of magnetic “highs” from 
Hamburg to Star City, which seems to be the northward continuation 
of the magnetic high trend from Winnsboro to Bastrop in Louisiana. 

The well known local anomaly at Rison is due to peridotite in- 
trusions with as much as 20 per cent magnetite content. 

The positive vectors from Stuttgart to Forest City seem to be- 
long to a large northwest-southeast magnetic high trend, toward 
which are also directed the negative vectors at Wynne and Helena. 
The large irregular positive and negative vectors farther north sug- 
gest very localized and relatively shallow sources of magnetism. 

Ouachita Mountain region—In the Ouachita Mountain region it 
is interesting to note that the two positive vectors at Mena and 
Mount Ida are almost horizontal and directed toward the Arkansas 
Valley, instead of south toward the axis of the Ouachita Mountains, 
as would be expected if the Ouachita Mountains were the result of 
uplift in the Basement complex. 

If it be assumed that the Basement complex is the cause of the 
magnetic anomalies in this region, the magnetic conditions suggest 
that this complex lies deeper below the Ouachita Mountains than 
below the Arkansas Valley, because the magnetic lines of force seem 
to come out of the ground at DeQueen and Murfreesboro, to assume 
a horizontal direction over the mountain range, and to penetrate into 
the ground in the Arkansas Valley. 

Thus the magnetic conditions seem to be in favor of the theory 
that the Ouachita Mountains are the result of a huge overthrust, as 
explained by W. A. J. M. van Watershoot van der Gracht.? Similar 


1 The local magnetic vectors are based on the information given in ‘United States 
Magnetic Tables and Magnetic Charts for 1925,” op. cit. 

2 W. A. J. M. van Watershoot van der Gracht, “(Permo-Carboniferous Orogeny in 
South-Central United States,” Bull. Amer. Assoc. Petrol. Geol., Vol. 15, No. 9 (Sep- 
tember, 1931), pp. 991-1057. 
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magnetic conditions exist in the region of the nappes of the European 
Alps. 

Arkansas Valley and Ozark Mountain regions.—The magnetically 
positive zones of the Arkansas Valley and Ozark Mountain regions 
are separated by a negative zone extending from Fayetteville to 
Marshall. 

OKLAHOMA! 


Eastern Oklahoma.—In eastern Oklahoma is to be observed the 
same arrangement of the magnetic anomalies as in western Arkansas; 
that is, the positive vector at Miami is indicative of the Ozark Moun- 
tains, the negative vectors near Westville are indicative of the mag- 
netic “low” between the Ozark uplift and the Arkansas Valley, which 
latter “‘structure”’ is represented by the positive vectors at Sallisaw, 
Muskogee, Stigler, and Poteau. 

The few stations in the general area of the Ouachita Mountains 
seem to confirm the interpretation given for the magnetic anomaly 
in the Ouachita region of Arkansas. 

On the basis of the west-east direction of the horizontal magnetic 
components, it seems possible to follow a positive magnetic trend 
from Nowata to Pryor and a second trend from Bartlesville through 
Tulsa, Sapulpa, and Okemah to Holdenville, along which trends lie 
the main oil “‘structures” of northern Oklahoma. 

Arbuckle Mountains.—The Arbuckle Mountains are indicated by 
a strong positive vector at Tishomingo. The positive vector at Durant 
might indicate the axis of the Arbuckle Mountains as extending south- 
east and northwest. 

Wichita M ountains.—The Wichita Mountains lie along a magnetic 
high trend beginning at Marietta and passing north of Waurika, 
Walters, and Hobart to Sayre. The oil fields of southern Oklahoma 
lie along this line. 

A magnetic high trend is indicated between Sulphur and Okla- 
homa City. 

Another high trend lies east of Rush Springs, Chickasha, and 
Minco. 

Another possible high trend lies between Alva and Carmen, be- 
tween Trail and Taloga, and east of Hammon. 

The southern extension of the Nemaha mountains might be indi- 
cated by the positive vectors at Medford, Pond Creek, and Enid. 


1 The local magnetic vectors are based on the informaticn given in “United States 
Magnetic Tables and Magnetic Charts for 1925,” op. cit. 
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KANSAS! 


Nemaha Granite ridge—The outstanding feature of eastern 
Kansas is the Nemaha Granite ridge, which extends from Seneca in 
a southwesterly direction toward Winfield. The top of the granite 
lies at +560 feet near Seneca and at —2,200 feet near Winfield, 
which indicates a southward dip of the ridge. From the Kansas State 
Geological Survey Bulletin 13 we quote: 

The earth movement, which elevated the granite appears to have taken 
place in the post-Mississippian and pre-Pennsylvanian time, for the upper 
part of the Ordovician and all of the Mississippian has been stripped away 
in the region of Eldorado and Augusta. North of Butler County, along the 


axis of the uplift, all of the sedimentary rocks older than the Pennsylvanian 
were removed so that the Pennsylvanian rests directly on the granite ridge. 
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If, therefore, the granite is embedded in lower Ordovician in the 
southern part of the ridge and in the Pennsylvanian in the northern 
part, we may expect this ridge to be indicated magnetically as a 
“low” on the south and as a “high” on the north if the Ordovician 
is of a higher and the Pennsylvanian of a lower magnetic permeability 
than the granite. 


1 The local magnetic vectors are based on the information given in “United States 
Magnetic Tables and Magnetic Charts for 1925,”’ op. cit. 
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But, though in central and western Kansas magnetic “lows” seem 
actually to correspond in some places with structural “highs,” the 
vector map seems to be in favor of the assumption that in the eastern 
part of the state magnetic “lows” correspond with structural “lows,” 
and magnetic “highs” with structural “highs.” 

The axis of the Nemaha ridge lies slightly east of Seneca, 5 miles 
west of Alma and Council Grove, slightly west of Cottonwood Falls 
and Eldorado, and between Winfield and Wellington. The ridge 
has a very steep eastern flank but slopes much more gently toward 
the west. Thus, the center of mass of the ridge lies considerably west 
of the structural axis, which would explain a corresponding shift of 
the magnetic axis. 

We think that the axis of the ridge is indicated by the positive vec- 
tors at Seneca, Council Grove, and Winfield-Wellington, and by the 
negative vectors at Alma and Eldorado, both of which point toward 
the west, so that the magnetic lines of force, emerging at Alma and 
Eldorado, may enter the positive ridge a few miles west of the sta- 
tions. 

Drilling has revealed a structural low trend along the east flank of 
the ridge and another extended “‘low,”’ with a west-east axis, between 
Eureka-Emporia and Yates Center-Burlington. The first trend is indi- 
cated magnetically by the negative vectors at Hiawatha, Holton, 
Alma, and Eldorado. In the second trend, it is clear that the negative 
vectors at Eureka, Yates Center, Iola, Burlington, Cottonwood 
Falls, and the positive vector at Emporia may all be assumed to 
emerge from the same magnetic “low.” 

We have pointed out already that in central and western Kansas 
both magnetic “lows” and “highs” may correspond with structural 
“highs.” 

The large and irregular vectors should yield considerable informa- 
tion, if the vector map is combined with magnetometer surveys and 
detail geologic maps. 

MISSISSIPPI! 


Many of the magnetic horizontal components of the vectors in the 
northern half of Mississippi have a northeast or a southwest direction, 
which may indicate a general northwest-southeasterly trend of the 
deeper “‘structure.”” We might suspect two magnetic high trends: (1) 
Butler-Philadelphia-Kosciusko-Greenwood-Rosedale on the south and 


1 The local magnetic vectors are based on the information given in “United States 
Magnetic Tables and Magnetic Charts for 1925,” op. cit. 
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(2) Carrollton-Columbus-Pittsboro- Water Valley, between Helena and 
Tunica on the north, and interpret them as indicative of the generally 
assumed connection of the Appalachian with the Ouachita Moun- 
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tains. The distances between the vectors and their irregularities are, 
however, too great to allow any definite statement, without addi- 
tional magnetic and geologic data. 
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The positive vector at Jackson indicates a magnetic “high” 
slightly southeast of the town. This “high” is due to an igneous plug, 
supposedly nepheline-syenite, and is in agreement with magnetometer 
surveys of the well known gas field. An interesting west-east trend of 
magnetic “highs” extends from Fayette to Ellisville. 

In general, there are hardly more than two vectors in Mississippi, 
which can be definitely attributed to the same structure. This circum- 
stance may be interpreted as due to local and relatively shallow 
sources of magnetism, which should be favorable for detailed investi- 
gations by the magnetometer. 


ALABAMA! 


From the northeastern corner of Alabama, a magnetic high trend, 
extending southwest to Greensboro, indicates the southern extension 
of the Appalachian Mountains. 

The vectors of the stations surrounding Greensboro are all positive 
and indicate by their direction that a large continuous positive anom- 
aly exists, covering that part of the state, which is explained as the 
“southernmost tip of the folded Appalachian Mountains” by D. R. 
Semmes.’ 

This “‘high” is an example of a large magnetic anomaly, corre- 
sponding with an extended structural feature, which could be cor- 
rectly interpreted by the study of the vertical intensity alone.® 

Another outstanding feature is the “low”? at Wetumpka. This 
“low” may lie along a low trend Center-Anniston-Talladega-Rock- 
ford-Wetumpka, with the high trend Guntersville-Greensboro on the 
west and another high trend on the east, which would pass near Ash- 
land, Dadeville, and between Tuskegee and Wetumpka. 

In the southern part of the state a low trend is indicated between 
Georgiana and Greenville, extending northwestward toward Camden 
and southeastward toward Elba. Another possible low trend lies east 
of Tuskegee, Union Springs, and Troy. 

The positive vectors at Opp, Geneva, and possibly also at Ozark 
and Abbeville appear to be due to localized “structures,” which might 
be of interest for oil and gas possibilities, because they lie in the Gulf 
Coastal belt. 

1 The local magnetic vectors are based on the information given in “United States 
Magnetic Tables and Magnetic Charts for 1925,” op. cit. 

2D. R. Semmes, Geol. Survey of Alabama Bull. 22. 

3 L. Spraragen, of. cit. 
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An interesting feature is the indicated high trend between D’Olive 
_Bay-Mobile-Lucedale on the south and Bay Minette-Citronelle- 
’ Leakesville on the north. This trend is parallel with the Hatchetigbee 
anticline, about 50 miles farther north. 

The Hatchetigbee anticline is not indicated by the vectors, be- 
cause the stations are not favorably located. It is parallel with another 
magnetic high trend about 30 miles farther north, extending from 
Butler toward midway between Grove Hill and Pine Hill. This last 
trend is known from magnetometer work in combination with the 
vector map. 

It seems, therefore, that there are three parallel pronounced anti- 
clinal trends, south of, and practically at right angles to, the Appala- 
chian trend from Guntersville to Greensboro. 

In the writer’s opinion, the three anticlinal trends in southwestern 
Alabama represent a mountain system distinctly different from the 
Appalachians. If the Appalachians extend farther than Greensboro 
and Linden, they are expected to disappear below the anticlinal trend 
at Butler. It is further suggested that the three anticlinal trends in 
southwestern Alabama possibly represent an eastern extension of the 
Ouachita Mountain system, inasmuch as the high trend at Butler 
might be followed through Philadelphia-Koskiusko-Greenwood-Rose- 
dale into the Ouachita Mountains, and the high trend at Mobile 
through Ellisville-Fayette-Winnsboro-Bastrop-Hamburg-Star City. 

In the northern part of Alabama we observe a pronounced mag- 
netic high trend Hamilton-Double Springs-Cullman. This trend seems 
to join the high trend Guntersville-Greenboro, and to extend east- 
ward between Stevenson and Fort Payne. In the writer’s opinion, the 
trend north of Fort Payne to Hamilton is representative of the main 
Appalachian trend, whereas the trend Guntersville-Greensboro is a 
southwestward extension, branching off the main trend at Gunters- 
ville. 

A possible westward continuation of the main trend of the Ap- 
palachians seems obscured by the magnetic influence of the Cin- 
cinnati uplift, which is indicated in the northeastern corner of Mis- 
sissippi by the “high” at Iuka. On the basis of magnetics only, we 
might suspect a high trend from Oxford toward a location midway 
between Helena and Tunica, with a westward continuation toward 
Brinkley and the Arkansas Valley. A continuation of the Arkansas 
Valley magnetic high trend into Oklahoma might be suspected south 
of Muskogee-Guthrie-Kingfisher-Watonga. 
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Additional magnetic and geologic information is needed, however, 
to clear up the possible existence and geologic meaning of this sus- 
pected trend. 

CALIFORNIA! 

It seems probable that the main magnetic effects in California are 
due to deep regional structure with a northwest-southeasterly trend. 
This interpretation explains the southwest and northeast directions 
of the most of the horizontal components, and is well established by 
the known regional geology. 

A magnetic high trend extends all along the California Valley from 
Red Bluff in the north through Stockton to Bakersfield in the south. 

A second parallel high trend may be noticed east of Bradley, San 
Lucas, Soledad, and Salinas. 

Farther south it seems possible to construct a high trend from 
Oceanside to Santa Monica and another parallel trend from San 
Clemente to San Nicholas. 

The magnetic stations are considerably denser in the valleys than 
in the mountainous regions of the Sierra Nevada and Coast Ranges 
or in the Mojave Desert. Inasmuch as the valleys represent sedi- 
mentary basins, whereas the mountainous regions are mainly com- 
posed of volcanic rocks, schists, gneisses, and highly metamorphosed 
sediments, it is surprising to note that the average vector in the val- 
leys is of about the same magnitude as the average vector in the 
mountainous regions. This circumstance may no doubt be explained 
in part by the geographic location and distribution of the stations, 
but it may possibly also have a deeper meaning, which the writer will 
try to explain. 

If the magnetic field is due to regional “structures” and if the cross 
sections in Figure 2 of ‘Decline of Great Basin’” are true, we should 
logically expect a magnetic high trend about 20 or 30 miles west of 
the San Andreas fault, over the Santa Lucia “blocks,” another high 
trend over the Sierra Nevada “block,” and a low trend in the Cali- 
fornia Valley, or sedimentary basin, between these two “‘blocks.” 

The cross section mentioned passes approximately through Sole- 
dad and slightly north of Madera, on our vector map. We notice a 
considerable difference between the actual magnetic data and those 
that would be assumed if based on the cross section. 


1 The local magnetic vectors are based on the information given in ‘United States 
Magnetic Tables and Magnetic Charts for 1925,” op. cit. 


2 J. Edmund Eaton, “Decline of Great Basin, Southwestern United States,” Bull. 
Amer. Assoc. Petrol. Geol., Vol. 16, No. 1 (January, 1932), p. 5. 
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The magnetic high trend indicated between Salinas and Hollister 
lies approximately along the San Andreas fault, that is, along the east 
flank of the Santa Lucia “blocks,” and the high trend in the California 
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Valley, a few miles northeast of Mendota, practically lies above the 
western edge of the Sierra Nevada “‘block,” as shown in the figure 
mentioned. 

There are not enough stations on the Sierra Nevada or Santa 
Lucia ‘‘blocks” to make possible any definite statement regarding 
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their magnetic behavior. From the magnetic vectors at Bradley, San 
Lucas, Soledad, and Salinas it seems, however, permissible to conclude 
that the Santa Lucia “blocks” would act somewhat neutral magneti- 
cally. If this should be the fact, then these “blocks” could not possibly 
reach great depths. 

It is interesting to note that the two pronounced magnetic high 
trends occur where they are least expected, that is, along the east 
flank of the Santa Lucia “blocks” and along the west flank of the 
Sierra Nevada “block.” If we interpret these high trends in the usual 
way, we should expect an anticlinal trend in the Basement complex 
below the California Valley, or along the west edge of the Sierra 
Nevada “block,” and another anticlinal trend in the Basement com- 
plex approximately below the San Andreas fault, or along the east 
edge of the Santa Lucia “blocks.” 

To interpret the magnetic data properly, we have to assume that 
the Santa Lucia “blocks”’ and the Sierra Nevada “‘block” are floating 
on a Basement complex, which is greatly depressed below the “‘blocks”’ 
and elevated along the edges of the “blocks.” This might be explained 
either by assuming that these two “blocks” are laccoliths, or by the 
assumption of huge overthrusts. 

The mechanics of the “progressive breaking down of the crystal- 
line Pacific borderland,” as explained by the arrows in Figure 2 and 


in the text of “Decline of Great Basin,” are not sufficiently free from 
possible objections, to speak plainly against the foregoing interpreta- 
tion.? 

In connection with a detail geologic map and scattered magneto- 
meter surveys, the vector map of California may prove of value also 
for the study of local anomalies. 


SUMMARY OF REGIONAL ANOMALIES 


Figure 10 summarizes the regional magnetic trends and areas of 
magnetic anomalies, which have been mentioned before in the dis- 
cussion of the respective states. We have discriminated between 
probable and possible anomalies. The probable anomalies are fairly 
well established by the magnetic vectors; the possible anomalies, how- 
ever, may upon more detailed investigations prove to be more local- 
ized and not as continuous as shown on the map, though the vectors 


1 J. Edmund Eaton, op. cit. 


? Friedrich Noelke, ““Geotektonische Hypothesen,” Sammlung geophysikalischer 
Schriften, No. 2 (Berlin, 1924), p. 70, b. 
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so far computed seem to be in favor of continuous trends, or at least 
do not disprove the trends based on regional magnetic or geologic evi- 
dence. 


HOO 


LIST OF REGIONAL ANOMALIES (FIG. 10) 


. High trend beginning between Stevenson and Ft. Payne and extending to Hamilton 
. High trend Guntersville-Greensboro 

. Low trend Center-Wetumpka 

. High trend Ashland to a location between Tuskegee and Wetumpka 

. Low trend Tuskegee-Troy 


Low trend Camden-Elba 


. High trend Butler to a location between Grove Hill and Pine Hill 

. Hatchetigbee anticline 

. High trend north of Mobile 

. Trend of “highs” from Ellisville to Fayette 

. Possible high trend Butler-Philadelphia-Greenwood-north of Rosedale-between 


De Witt and Varner 


. High area near Iuka 

. Possible high trend Oxford, Water Valley to a location between Helena and Tunica 
. High trend Marianna-Brinkley 

. High trend of Arkansas Valley 

. Low trend Marshall-Fayetteville-Tahlequah 

. “High” at Miami and Ozark region 

. Neutral trend parallel with axis of Ouachita Mountains from area between Hugo 


and Antlers toward area south of Mena and Mount Ida to Sheridan 


. High trend Hamburg-Star City 

. High trend Winnsboro-Bastrop 

. Possible high trend Farmerville-Harrisonburg 

. Low trend Ruston-Winfield 

. “Highs” south of Shreveport and Many, with connecting trends parallel with axis 


of Sabine uplift 


. Possible “high” northeast of Marksville 

. High trend between Ville Platte and Opelousas 

. High trend between Morgan City and Napoleonville 

. High trend Lafayette-Lake Charles 

. “High” south of Liberty, regional trend 

. “High” south of Houston, regional trend 

. “High” at Refugio, regional trend 

. Possible high trend Point Isabel-Corpus Christi 

. Possible high trend east of Edinburgh-Falfurias-San Diego 

. Possible high trend south of Newton-Conroe-Columbus 

. “High” at Lampasas, Llano, Burnet, corresponding with the Llano-Burnet uplift 
. “High” at Sonora, Ozona, Eldorado 

. “High” from Brady to Jacksboro, corresponding with the Bend arch 

. Possible high trend north of San Angelo 

. “High” indicated by the vectors at Rankin and Stiles, probably corresponding with 


the Big Lake anticline 


. “High” south of Monahans, with trend corresponding with Central Basin plat- 


form 


. “Low” near Pecos, corresponding with Delaware basin 

. Low trend, corresponding with the axis of main Permian basin 

. “Low” between Bovina, Tulia, Canyon, interpreted as basin 

. High trend Marietta~-Lawton-Hobart-Sayre, corresponding with Wichita Moun- 


tains 


. High trend Tishomingo-Durant, corresponding with Arbuckle Mountains 


. High trend Sulphur-Oklahoma City 
. High trend Rush Springs-Minco 


. Possible high trend Fort Smith-Muskogee-Watonga 
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. Pcssible high trend Nowata-Pryor 

. Possible high trend Bartlesville-Tulsa-Holdenville 

. Pessible high trend east of Medford-Enid, corresponding with possible south 
extension of Nemaha Mountains 

. Possible high trend east of Alva-Trail 

. Possible high trend corsponding with Nemaha Mountains 

. Low trend east of Nemaha Mountains 

. Low trend between Yates center and Burlington 

. “Low” at Ulysses 

. High trend along California Valley 

. High trend along Salinas Valley 

. High trend from Oceanside to Santa Monica 

. Possible high trend San Clemente-San Nicholas 


/ 

48 
49 
50) 
51 
52 
53 
54 
55 
56 
57 
58 
5 


BULLETIN OF THE AMERICAN ASSOCIATION OF PETROLEUM GEOLOGISTS 
VOL. 16, NO. 12 (DECEMBER, 1932), PP. 1204-1211, 4 FIGS. 


APPLICATION OF REFLECTION SEISMOGRAPH' 


SUGENE McDERMOTT? 
Dallas, Texas 


ABSTRACT 


This paper describes in detail a reflection seismograph survey of western Henderson 
and eastern Navarro counties, and covers the area just north of the Mexia-Powell oil 
tield. A number of faults are located and closure is indicated on two. A cross section of 
the area shows the relation of the reflection records to the subsurface structure. Re- 
flections were obtained throughout the area on the Pecan Gap and a basal member of 
the Austin chalk. Contour maps on these two reflection horizons are included. 


INTRODUCTION 


The area surveyed in this application of the reflection seismograph 
lies just north of the Mexia-Powell fault-line fields in eastern Navarro 
and western Henderson counties, Texas, and has long been considered 
a favorable area in which to explore for a possible extension of this 
line of fields. At least twenty dry holes have been drilled in this area. 


A seismograph survey such as the one under consideration would have 
made possible the selection of the best location, structurally, so that 
one, or at most two, wells would have thoroughly tested the area for 
the presence of oil. 

A total of approximately seventy-five depth determinations were 
made in a period of about two weeks, resulting in a cost per acre of 
somewhat less than ro cents. 

The author has attempted to present the data in the form of two 
contour maps, one of the top of the Pecan Gap chalk and the other 
of the basal member of the Austin chalk, a cross section with corre- 
sponding reflection records from which the depths to these two re- 
flecting horizons were made, an enlarged view of a typical record 
showing the relation between the reflection events and the section, 


and has endeavored to make these as nearly self-explanatory as pos- 
sible. 


1 Presented before the Association at the Oklahoma City meeting, March 25, 1932. 
Manuscript received, June 4, 1932. 


2 Geophysical Service, Inc., 1311 Republic Bank Building. 
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PROCEDURE 


As the trend of any faults that might exist was assumed to be the 
same as that of the field on the south, east-to-west lines of depth de- 
terminations were laid out, taking advantage of roads available. When 
the data indicated the presence of a fault, the layout was changed so 
as to increase the density of shot locations in the neighborhood of the 
fault and follow its trend. The program was essentially of a recon- 
naissance nature. 


CORRELATION OF REFLECTIONS 


Two dominant reflections appear on the records. The interval be- 
tween them corresponds to that between the top of the Pecan Gap 
chalk and the basal Austin. A standard Cretaceous velocity chart was 
used in the computation of the depth determinations. Although this 
could not be expected to give correct absolute depths beyond approxi- 
mately 100 feet, as it was constructed from well determinations in 
another part of East Texas, interval and relative-depth determina- 
tions made from it could be relied upon to be correct to 0.5 per cent. 
The two contour maps show the interval between the Pecan Gap and 
the Austin to be approximately constant but indicate some regional 
variation. 

A typical reflection record showing the relation of the reflection 
events to the geologic section is shown in Figure 1. The fact that a 
dominant reflection is obtained from near the base of the Austin chalk 
would indicate that the basal member of this body of chalk is consider- 
ably harder than the rest of the chalk section. Occasionally a reflection 
is obtained from the upper part of the chalk section, but it is never as 
pronounced as the basal reflection nor does it appear on the records 
with any regularity. In view of the dominant character of the two 
reflections and the regularity of interval between them, the possi- 
bility of miscorrelation is extremely small. 

The two contour maps, Figures 2 and 3, were constructed from 
reflection depth determinations entirely. No well data were used. Al- 
though the contour maps may not check similar maps made from 
well data in absolute value, they should check as regards relative data. 
Most of the wells drilled in this area are shown on the contour maps. 

A cross section showing two faults is presented in Figure 4. This 
cross section is along line AA indicated on the contour maps. The re- 
flection records from which the depth determinations were made are 
shown opposite the cross sections. Each record corresponds with the 
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depth determinations opposite it on the cross section. The two re- 
flections, Pecan Gap and basal Austin member, are marked on the 
records and the datum values shown. These datum values are shown 
to scale on the cross section. 

Due to the greater depth of the Austin, the reflection from it 
should be smaller in amplitude than that from the Pecan Gap. For a 
reasonably large Pecan Gap reflection the Austin reflection amplitude 


NAVADQO COUNTY - TEXAS 


Fic. 1.—Typical reflection record showing relation of reflections to geologic section. 


is almost too small to use. It is, therefore, general practice to take 
several records with different dynamite charges. In this way both 
large- and small-amplitude records are obtained for the two reflec- 
tions. The small-amplitude Austin reflections shown in Figure 4 have 
been further substantiated by larger amplitude records. 

Attention is called to the absence of the Austin reflection on record 
F.P.210. It so happens that the normal reflection point on the Austin 
falls in the Austin fault plane. In view of the extraneous disturbance 
in this record at the time of normal occurrence of the Austin reflec- 
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PECAN GAP CONTOURS 
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Fic. 2“—Reflection contours on top of Pecan Gap chalk. 
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CONTOURS ON TOP OF HARD MEMBER 
OF THE AUSTIN OWLK 


Fic. 3.—Reflection contours on top of basal member of Austin chalk. 
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tion, the absence of the reflection in this instance can hardly be called 
conclusive evidence of the presence of the fault. 

Due to the absence of the Pecan Gap reflection, record F.P.58 is of 
special interest. As the amplitude of the Austin reflection on this 
record is larger than on any of the other records, the Pecan Gap re- 
flection under normal conditions should be very large indeed. There 
can be little doubt that this must be the result of very abnormal sub- 
surface conditions. It is seen from the cross section and contour maps 
that the normal reflection point for the Pecan Gap lies in the fault 
plane of the Pecan Gap. This fault, it should be noted, has been lo- 
cated by several other determinations. In many cases it is impossible 
to get any reflections very close to a fault, due probably to minor 
fracturing in the proximity of the fault plane. The presence of the 
Austin reflection in this case would indicate that the fault must be a 
fairly clean break without much minor fracturing. 


CONCLUSION 


It is thought that the reflection seismograph may economically 
justify its use in many cases by permitting the location of a test well 
in the most advantageous position structurally and thus save the cost 
of many unnecessary wells. 
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Fic. 4.—Cross section AA and corresponding reflection records. Horizontal length of section, 30,500 feet. 
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USE OF RECORD CHARACTER IN INTERPRETING 
RESULTS AND ITS EFFECT ON DEPTH CAL- 
CULATION IN REFRACTION WORK! 


T. L. ALLEN? 
San Antonio, Texas 


ABSTRACT 


In seismograph exploration, the author has encountered acoustic horizons which do 
not transmit seismic energy readily. The first impulse on the seismogram has a short 
period with small amplitude, and requires considerable energy to make it apparent at 
long distances. This horizon is underlain by one of slightly higher velocity which trans- 
mits the energy readily and has a distinctly different wave form. On account of the 
difficulty of energy transmission, the horizon which has the small “forerunner” wave 
might be interpreted as a lens and not a continuous horizon. By making use of the 
character of the wave form on the seismograms, the two horizons may be definitely 
identified, and the one which is a poor conductor of energy definitely proved to be a 
continuous acoustic horizon. 


INTRODUCTION 


This paper covers a short discussion of a problem encountered in 
refraction work done in Venezuela during 1931. The work was done 
by the International Petroleum Engineering Corporation of New 
York, using electrical seismographs and radio equipment built by the 
Petty Geophysical Engineering Company of San Antonio, Texas. The 
writer was in charge of the party in the field. 

The problem to be considered involved an interpretation of re- 
sults in an area where an acoustic horizon transmitted the seismic 
energy poorly. Referring to profile A, Plate 1, it is evident that there 
are four acoustic horizons involved, with velocity values of approxi- 
mately 2,000, 2,160, 2,580 and 6,510 meters/second, the first three 
being sedimentary, while the fast one represents the igneous base- 
ment. There were several other profiles shot in the same general area 
as profile A and all of them had the same general velocity values. 
Profiles B, C, and D were approximately 50 kilometers from profile A, 
and the sedimentary horizons had become thinner in this distance. 
Several profiles were shot in the distance between the two areas, and 


1 Published by permission of International Petroleum Engineering Corporation, 
80 Maiden Lane, New York, N. Y. Read before the Association at the Oklahoma City 
meeting, March 24, 1932. 


Sie 2 Field manager, Petty Geophysical Engineering Company. Introduced by L. W. 
u. 
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their results showed a gradual thinning of the sedimentary beds and 
a gradual decrease in the velocity values. The seismograms shown in 
Plate 2 are the ones from which profile B is plotted, Plate 3 shows the 
seismograms from which profile C is plotted, and Plate 4 shows those 
from which profile D is plotted. 


PLaTE 2.—Seismograms for profile B. 


STUDY OF RECORDS 


A study of the three sets of seismograms reveals four distinct 
record characteristics. The first five seismograms of profile B, Plate 2, 
break sharply and plot on the surface velocity curve; the next four 
seismograms do not break sharply but have a small wave of short 
period for the first impulse followed by a wave with a sharp break 
and large amplitude. The small impulse and the sharp break were 


INTERPRETING RESULTS ON DEPTH CALCULATION 1215 


both considered and are plotted on the profile. The sharp break with 
large amplitude is the first impulse on records 2B and 2B’, while on 
the longest three records the sharp break is preceded by an absorbed 
impulse which is characteristic of the igneous basement. Several other 
profiles in the same general area as profile B all had a small wave of 
short period on the second horizon which required large dynamite 
charges to bring out. A curve drawn through the points represented 
by the first small impulse or “forerunner” does not intersect the curve 
drawn through the points representing the sharp breaks of the third 
horizon until after the basement velocity curve intersects both of 
them. Should the small wave or “forerunner” be interpreted as a thin 
lens-like formation overlying the horizon which gives the sharp break, 
or is it a separate acoustic horizon? If the former, it might be dis- 
regarded as far as depth calculation is concerned, but if it is the latter 
its thickness must be considered. 


THEORETICAL CONSIDERATION 


If the conditions shown in Figure 1 be assumed, it is possible to 
derive equations for the time-distance curves which these velocities 
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. 1.—Assumed conditions of depth and velocity. 
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and thicknesses will give. The distance of the receiving station from 
the shot point will be designated as x, and the equations for the time 
through each paih derived. For the surface layer, the general equa- 


tion for the time is, 


x 
Ti = 
1.980 


The general equation for the time through the second horizon is, 


= 
* cos 69°25’ X 1.980 2.115 


For the third horizon, it is, 


2X .330 2X .322 
= 


cos 60°03’ X 1.980 cos 67°46’ K 2.115 
x — 2X .330 X tan 60°03’ — 2 X .322 X tan 67°46’ 


2.285 


and for the fourth horizon, it is, 


T, = 


2X .330 4 2X .322 2X .981 

cos19°21’ X 1.980 cos20°44’X 2.115 cos 22°29’ XK 2.285 
I 
[x — 2X .330 X tan 19°21’ — 2 X .322 X tan 20°44’ 

5-975 
— 2X .981 X tan 22°29']. 
These are all straight line equations which may be reduced to the 
simple forms which follow: . 


T\ 


1.980 
= 
2.115 
3.990 
1.45387 
2.285 
x — 1.288 
T, = 1.608 + —————__ - 
5-975 
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The curves represented by these equations are shown in Figure 2. It 
is apparent that the wave from the third horizon would never be the 
first impulse, since the intersection of the 72: and 7; curves is past the 
point where the 7, wave arrives first. The only manner by which the 
third horizon velocity may be determined is to select the secondary 
waves on records shorter than 4.174 kilometers, which is the value of 
x where the 7; and 7, curves intersect in this particular case. 


Seconds 
Seconds = 


Kilometers 


4 


Fic. 2.—Time-distance curve for assumed conditions. 


This theoretical consideration seems to indicate that the small 
wave or “forerunner” represents a distinct acoustic horizon and not a 
lens-like formation. If this be true, it should always give the first 
seismic impulse on the seismogram, provided there is enough energy 
available to bring it out. In Plate 3, there was more energy available 
on account of the shorter distances and the amplitude of the first im- 
pulse is quite large on all the seismograms except 2D. The period of 
the wave is quite short, however, indicating that it is the impulse 
characteristic of the second horizon. On Plate 4, records 1A and 1A’ 
have a sharp break as the first apparent impulse, but a study of these 
seismograms under a glass shows a slight movement ahead of the ap- 
parent initial impulse. Record 2C was only about roo meters longer 
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than 1A’ but had four times as much dynamite, and the small fore- 
runner wave is very evident on this record. This proves conclusively 
that the formation which transmits the small forerunner type of wave 
is a continuous acoustic horizon which does not transmit the seismic 
energy readily. 


CONSIDERATION OF DIFFERENT PHASES OF THE PROBLEM 


The problem presented here might be thought of from two differ- 
ent angles. If all horizons had been of equal energy-carrying ability, 


PLATE 3.—Siesmograms for profile C. 


the third one would have been difficult to detect, due to the fact that 
the particular thicknesses and velocities are such that the third hori- 
zon impulse would never have been first. On the other hand, it would 
be easy to neglect the beds which transmitted energy poorly unless 
sufficient dynamite was used to make the first impulse large enough 
to be seen. In the area where this work was done, the dynamite charge 
was chosen that would give just enough amplitude to the second 
horizon wave to make it apparent, yet not so large that it obscured 
the large secondary impulse from the third horizon. In this manner, 


‘ 
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it was possible to choose a primary wave and designate it as the 


second horizon impulse and also choose a secondary wave and con- 
sider it as a primary wave from the third horizon. If the second hori- 


PLATE 4.—Seismograms for profile D. 


zon had been omitted from the calculations, the total depth to the 
basement would have been too small and the surface horizon would 
appear much thicker. Omission of the third horizon from the calcula- 
tion would also have made the total depth to the basement too small 
but would not have changed the thickness of the surface horizon. 


T. L. ALLEN 


REASON FOR POOR ENERGY-CARRYING CHARACTER 


The reason for poor energy transmission is probably due to al- 
ternating beds of hard and soft material, perhaps sandstones and 
clays. There is an area in the Gulf Coast where a poor energy-carrying 
horizon is found, and it is especially noticeable in Harris County, 
Texas. The movement on the seismogram is a small wave whose be- 
ginning is scarcely perceptible and is sometimes mistaken for a salt 
dome ‘“forerunner.”’ Recent work by the reflection method in the 
Gulf Coast region shows several reflections from the horizon which 
transmits the refracted wave poorly, indicating a change in density 
which would logically be explained by alternating beds. If the al- 
ternating beds which compose the acoustic horizon had slightly 
different velocity values, the refracted wave would gain a slight 
amount in the faster bed, and be gradually damped out. 


CONCLUSION 


The seismic wave accurately shows the subsurface conditions. The 
character of the wave is affected by the various strata through which 
it passes, and a study of the wave character is important in the cor- 
rect interpretation of seismic results. The possibilities of the seismo- 
graph as an accurate exploratory instrument have not been fully 
realized, and this paper is presented with the hope that further study 
along similar lines will broaden the scope of usefulness for this in 
strument. 
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SEISMOLOGICAL DISCOVERY AND PARTIAL DETAIL OF 
VERMILION BAY SALT DOME, LOUISIANA! 


E. E. ROSAIRE? and O. C. LESTER, JR 
Houston, Texas 


ABSTRACT 


The Vermilion Bay salt dome in Iberia Parish, Louisiana, was one of the first 
discovered by a geophysical survey of water-covered areas. The field organization and 
equipment of such a crew is described. The discovery and detail data are presented, 
along with the original interpretation of these data. The results of drilling are shown, 
and the agreement and disagreement of these tests with the original interpretations 
discussed to some extent. 


LOCATION 

Vermilion Bay is one of the major features of the Louisiana coast 
line. The Vermilion Bay salt dome is located in the southeastern 
corner of the bay, but in the extreme southern portion of Iberia 
Parish, Louisiana. It is best reached by boat along the Boston Canal 
from Abbeville, or by car to Cypremort Point, and thence by boat to 
the dome. 

DISCOVERY 

This salt dome was discovered on September 25, 1927, in a recon- 
naissance refraction exploration conducted for the Louisiana Land 
and Exploration Company by a Geophysical Research Corporation 
party* under the direction of the junior writer. 

The discovery of this dome was of importance beyond the usual, 
for several reasons. Among these may be listed the following. 

1. This salt dome was discovered as a result of the use of the 
seismograph in exploring large bodies of water. The Louisiana Land 
and Exploration Company was the first to launch a major exploration 
campaign with this definitely in view, using for this purpose two Geo- 
physical Research Corporation parties,® equipped for water surveys. 


1 Read before the Association at the Oklahoma City meeting, March 25, 1932. 
Published by permission of the Louisiana Land and Exploration Company, and the 
Geophysical Research Corporation. 


2 Consulting geophysicist, 2210 Esperson Building. 
3 Geophysical Research Corporation. 


‘ Within a month, this same crew had also discovered the Four Isle and Dog Lake 
salt domes, in southwest Terrebonne Parish. 

5 The other crew, under the direction of C. G. Rosaire, was the first to discover a 
salt dome during a water survey, since his discovery of the Calcasieu Lake salt dome 
antedated O. C. Lester’s discovery of the Vermilion Bay dome by several days. 
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By virtue of their being first to embark upon this virgin field, these 
two refraction crews staged what is probably the most spectacular 
campaign in the history of geophysical exploration, discovering no 
less than eleven salt domes in 9 months. These successes, starting 
similar campaigns by other companies, caused the intensive explora- 
tion (by seismograph crews) of southeastern Louisiana, and further 
resulted in the discovery of six additional salt domes in this same 
area, all within a year after the discoveries of the Calcasieu Lake 
and Vermilion Bay salt domes. 

2. This successful application of the seismograph to water ex- 
plorations undoubtedly brings part, at least, of the Gulf of Mexico 
within the available petroleum reserves of this geological province, 
if or when the value of crude oil justifies the solution of such practical 
and legal complications as are apparent. 


OPERATIONS 


Vermilion Bay varies in mean depth fron: 1o to 13 feet in an area 
of approximately 125,000 acres (Fig. 1). The operations in this in- 
stance, as in other water explorations of similar magnitude, were con- 
ducted entirely from boats. The crew was quartered in a houseboat 
with auxiliary barges for explosives and supplies, while the recording 
apparatus was mounted in fishing luggers 40 to 75 feet in length. A 
similar lugger was used in locating, planting, and firing the charges 
of explosives. 

The explorations were conducted by the usual method of refrac- 
tion fans. Communications and transmissions of time signals were 
effected by two-way radio operations on 180 meters, using 5-watt 
IC W transmitters. Distances from shot to recorders were determined 
from the air travel times of the explosive sound waves. 

Though the analysis of such exploration data is comparatively 
simple, its acquisition involves the solution of many practical diffi- 
culties. Adverse winds may carry the sound wave above the recorders, 
or may lower the water level to such an extent that operations of any 
kind are impossible over the subsequent mud flats; and may even 
maroon the recording boats for as much as a week. In the latter case 
the resulting diet of canned tomatoes and oysters becomes rather 
monotonous after a few days. On the other hand, moderately rough 
water is no great handicap, since seasick observers have taken excel- 
lent records with the seismographs planted in water and mud at 
depths of 15-20 feet. Also, as this coast is not infrequently visited 
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by tropical hurricanes, the field crews watch the Coast Guard Weather 
Reports (received by radio) and occasionally have to run for land, re- 
turning later to salvage gasoline drums and boxes of dynamite scat- 
tered over the shore by sudden and severe squalls. 

Even under such conditions, the explorations averaged 200,000 
acres per month; and as much as 15,000 or more acres have been sur- 
veyed in a single operating day. 


VERMILION BAY SEISMIC SURVEY 
ey 


GEOPHYSICAL RESEARCH CORPORATION 
10th te 


For the Lend & Exploration Co. 


Fic. 1.—Refraction fan layout for exploration of Vermilion Bay. 


In the survey of Vermilion Bay, three recording boats were used. 
Each was equipped with a radio transmitter and receiver, an oscillo- 
graph, and a velocity meter (geophone) recording through a vacuum 
tube amplifier. The timing device was an electrically operated 50- 
cycle tuning fork. The air wave was picked up by a carbon-button 
unit. The resulting records were developed immediately, thus guiding 
the observers in the size of the next shot. The instant of explosion was 
determined by the breaking of a wire wrapped around part of the 
charge of explosive, thus terminating the operation of the shot trans- 
mitter. In general, the charge of explosive was simply lowered to the 


\ 
> 


1224 E. E. ROSAIRE AND O. C. LESTER, JR. 


bottom of the water, or buried in a hole on shore. A smaller simul- 
taneously detonated surface charge was used for the air wave dis- 
tance determination. Travel times were read to hundredths of a 
second. 

Distance determinations were made from the travel times for the 
air waves, corrections being used for wind velocities and directions, 
determined by compass and hand anemometers at the time of each 
shot. For short-shot lengths, up to about 5 miles, these were relatively 
accurate to about 200 feet. (At greater distances over land, large errors 


Fic. 2.—Travel time data for Vermilion Bay. 


appeared in this method, as for instance, when an error of 2,000 feet 
appeared in a shot line 40,000 feet long.) 

Distance determinations and final travel-time counts were made 
at headquarters in the houseboat. The travel-time data (Fig. 2) were 
plotted, and the normal relationship indicated by drawing the obvious 
“normal” travel-time curve. Travel times appreciably shorter than 
normal were then considered as suggesting the existence of a salt 
dome, and bona fide time differences checked by a cross fan. 


INTERPRETATIVE PROCEDURE 


Long distance refractions.—The data collected in Vermilion Bay 
may be taken as typical of a shallow Gulf Coast salt dome. For con- 
venience of study, the time intervals are plotted, on each line, from 
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an arc indicating normal arrival times. If minus time differences are 
plotted out from the shot point, an approximation is made to the 
outline of the wave front at the recording points. In the case of the 
Vermilion Bay data, the resulting pattern of time differences un- 
mistakably indicated the presence of a previously unknown salt 
dome. 
The smooth curve used as the ‘‘normal”’ in Figure 1 is calculated 
from the travel-time relationship 
2b aT 
X = —sin h— 
a 2 
for a linear increase of velocity with depth. The constants b (surface 
velocity) and a (increase in velocity per foot of depth) were deter- 
mined by trial, and found to be 
b = 5,633 feet per second 
a = 0.5483 foot per second. 


The penetration, Z, for such a curve is 


Short refraction data.—A short negative profile was shot north 
from SP 1, to serve as a “‘normal’’ in mapping the top of the dome. 
This is shown in Figure 3. 

An interesting feature of this profile is that the extrapolated 
travel-time curve indicates a negative time for zero distances. This 
is to be expected for explosion-generated data, and indicates the ab- 
normally high velocity of sound in the immediate neighborhood of 
the explosion. However, in comparison with the greater number of 
such short refraction profiles, these data are exceptional, since the 
usual data indicate a plus time at zero distance, due to the general 
existence of a low-speed layer at the surface of the ground (weathered 
or aerated zone). In the case of the Vermilion Bay data, there is ap- 
parently no aerated zone,' so that there is opportunity to observe the 
effect of the abnormally high sound velocities in the immediate 
vicinity of the explosion.? 


1 This is not necessarily an obvious conclusion, since in some cases of water- 
covered marsh, definite geophysical evidence of an aerated zone has been found on the 
shorter travel-time curves. As Paul Weaver has pointed out in “Geophysics of the 
Soil,” this is to be expected in the case of thick layers of soil being formed from rotting 
vegetation. 


2 A. B. Wood, Textbook on Sound, p. 266. 


Z= 
4 a? a 
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This short profile can also be fitted by a travel-time curve such as 
that approximated by the longer refraction data. For the shorter data, 
the constants are 

b= 4,990 feet per second 
a= 1.25 feet per second per foot. 


The difference in the values of these “‘constants’’ for the two sets of 
data is significant, since it is an indication of the extent to which 
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Fic. 3.—Short refraction data for Vermilion Bay. 


these data can be approximated by a linear increase of velocity with 
depth. 

Partial refraction detail._—The location of the dome by refraction 
fans was recognized as too qualitative for drilling. Therefore, with the 
thought of securing more definite information as to the depth and 
location of the top and flanks, the following detail survey was carried 
out. 

Figure 4 shows the location of these refraction profiles, together 
with the resulting interpretation of the outline. The profile data are 
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shown in Figure 5 and Figure 6. The solid part of the travel-time 
curves are taken from the negative profile data described in Figure 3. 


Cypremort Point 


Cotcwintes depth ot point 700 


Marah ~ 


Fic. 4.—Refraction profile layout for partial detail of Vermilion Bay salt dome. 


All four of these profiles show the usual features of a refraction 
profile characteristic of a shallow salt dome. The most essential feature 
is the more or less sudden decrease in apparent velocity shortly after 
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the recording points pass well beyond the salt flank. The junior writer, 
at the time of this detail, made a depth calculation of 800 feet to cap 
or salt, at the intersection of the profiles. The agreement with the 
results of the well later drilled at this point was very satisfactory. 
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Fic. 5.—Southwest-northeast refraction profile data, Vermilion Bay salt dome. 


Fic. 6.—Northwest-southeast profile data, Vermilion Bay salt dome. 


The original interpretation of the outline of the dome (dotted) 
was made in the New York office, not in the field. The southeast flank 
location was good, as shown by the findings in the wells drilled on the 
line. The only other check of this outline was made on the northeast 


/ / 

4 | 
| SSS SSS SS SSS SS SSS SSS SSS SSS | 

lat | T 
| | | 
SSS SS SSS SSS SSS ARS SSS SS 
| 
+ —+ + + + + + + + 4+ + + + + — ++ + + + + —+-—+- — 
| t | 
+ + + ret 

ESS SSS SSS 

SSS 
SS SSR SS SSS SS VSS SS SS SSS SS SRS SS SSS 
T T | 
Sees et ee ees eevee ee eee 

in 
Profite 


SEISMOLOGICAL DISCOVERY OF VERMILION BAY 1229 


flank, with disappointing results. Well No. 5, supposedly a location 
similar to Well No. 2, was abandoned in rock at 6,500 feet. Well No. 
6, well inside this estimated edge, was abandoned at 5,824 feet with- 
out finding dome materials, when cap or salt was expected not deeper 
than goo feet. A recent revision of the interpretation on this flank was 
made by the junior writer, and led to the location of this edge as shown 
by the dashed line. As no further wells have been drilled since No. 6, 
the original interpretation of the northwest and southwest flanks, 
and the re-interpretation of the northeast flank remain unchecked 
to date. 

With the experience accumulated, since this detail, on several 
other domes, the writers are in a position to comment on the job. It is 
obviously incomplete, since four points, no matter how well deter- 
mined, can hardly furnish satisfactory control on the periphery of a 
salt dome like that in Vermilion Bay. Later practice was to locate 
eight points, reversing all profiles. These additional points, determined 
on stronger evidence, resulted, in the case of the Lake Barre salt 
dome, in commercial production by the first flank test, followed by a 
string of nine successive flank producers around the periphery of the 
dome. 

In addition, the profile data on the northeast flank somewhat 
definitely indicate that this flank slopes less steeply than any of the 
other three. Failure to realize the amount of this slope (as shown by 
the apparent velocities) undoubtedly played a major part in the 
erroneous location of this flank. 

Finally, these conclusions were reached on profile data with too 
low a survey density. Recorder intervals of 1,500 feet can hardly be 
expected to locate salt dome flanks to plus or minus 250 feet. In the 
later details, these recorder intervals were cut down to as little as 250 
feet, along the profile. 

In presenting the data for publication, attention is called to the 
fact that curved lines have been used to approximate the travel-time 
curves. At the time of this survey, field practice was to use succes- 
sive, intersecting straight lines, but for the past 3 years, curved path 
interpretations have entirely replaced the older, empirical, strong- 
arm methods of attack in almost all phases of refraction surveys. In 
view of this, the writers feel justified in restricting their case treat- 
ment to the earlier predictions, and so indicating, to one experienced 
in the art, the contrast between the old and the new interpretation 
methods. 
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SEISMIC WEATHERED OR AERATED 
SURFACE LAYER' 


O. C. LESTER, JR.? 
Tulsa, Oklahoma 


ABSTRACT 


The existence of a comparatively thin surface layer with a low velocity character- 
istic has been recognized in seismic work for several years, and has been generally 
referred to as the “weathered layer.” This low velocity surface characteristic has been 
found to be almost universally present regardless of the nature of the surface deposits, 
and does not conform to the geologic weathering of the area. 

The purpose of the writer is to offer as an explanation of this phenomenon, the 
mixture of air in a free state with surface materials. Theoretical calculations, if the 
earth is assumed to be a fluid, indicate that velocities less than that of sound in air 
should be obtained from such mixtures. This is borne out by experimental data, as far 
as available, indicating that this so-called weathered layer might properly be termed 
an “aerated” layer. 


It must be said at the outset that the conclusions of this paper are 
based on meager field data. The primary argument arises from theo- 
retical study and postulates certain effects. These effects, however, 
have been confirmed by experiment in so far as such experimental 
data are available. The purpose of the discussion is, first, to present a 
probable explanation of a phenomenon which has been recognized 
and dealt with—though unexplained; and second, as a result of such 
explanation, to suggest a more descriptive and less confusing name 
than that of “weathered.” 

The existence of a comparatively thin surface layer with a low 
velocity characteristic has been recognized in seismic work for several 
years, and it has been generally referred to as the “weathered layer.”’ 
This condition has been found to be almost universally present, re- 
gardless of the nature of the surface deposits, and has been recognized 
from its low velocity characteristic rather than by any physical or 
geological difference between it and the materials immediately below. 
The drilling of deeper shot holes, below the depth of this low velocity 
horizon, has failed to indicate any visible differentiation between it 
and adjacent materials. 


1 Read before the Association at the Oklahoma City meeting, March 25, 1932. 
2 Geophysical Research Corporation, Box 2040. 
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Since the recognition of this stratum, various methods have been 
applied to correct for the abnormally slow travel time through it 
(particularly in reflection work), but little or no attempt has been 
made to explain its nature, or to determine the cause of its important 
characteristic. 

Data taken in various places, over a wide area, from Kansas and 
Oklahoma to the gulf coasts of Mississippi, Louisiana, and Texas, 
have indicated that this layer may vary in thickness from practically 
zero to more than 150 feet, though the general average is approxi- 
mately 50 feet. The velocity of sound in this “weathered” material, 
though variable, is always comparatively low and a general average 
might be placed between 2,000 and 2,500 feet per second; while the 
average velocity of sound in the same materials immediately below is 
found to be generally about 5,600 feet per second. 

As this top layer is at the surface of the ground and exposed, as 
its low velocity characteristic exists almost universally regardless of 
composition, and as this low velocity characteristic could result from 
a broken, disintegrated or loosely bedded condition, the term ‘‘weath- 
ered layer” is perhaps natural. This term has been in use for some time 
for the sake of convenience and for lack of a better one, though it has 
been granted that this layer does not, in many cases, conform to the 
geologic term “weathered.” In the Gulf Coast region, for example, 
the geophysical weathered layer is exceptionally uniform and thin, 
though the geologic term “‘weathered” could be applied to the sedi- 
ments at considerable depths. 

The suggestion that the “weathered layer” might be associated 
with the regional ground-water table has often been made, but though 
in some areas the thickness of this material seems to be determined 
by the depth of the water table, this has not been found to be true in 
a sufficient number of cases to permit a general conclusion. 

The extremely low velocities (approaching, or even less than, the 
velocity of sound in air) that have been obtained in some areas, led 
to the idea of air inclusion in the surface materials as the cause.' 

With this idea in mind it was predicted that short (shallow) re- 
fraction profiles would show lower and lower velocities with shallower 
depth penetrations, approaching if not actually equaling the velocity 


The discussion here refers to air mixture in a free state—not in solution. As 


= . air or gases may be considered together, as opposed to fluids or solids, be- 


cause of the order of magnitude of difference in E and P for these two classes of ma- 
terials, E=elasticity: P=density. 
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of sound in air, as the air became present in larger and larger propor- 
tions near the surface. The few profiles in which velocities lower than 
that of sound in air had been indicated, were considered interesting 
but somewhat anomalous, and were strongly suspected of contain- 
ing some unseen error. A discussion of sound velocities in air-water 
mixtures,' which was pointed out by E. E. Rosaire, was the direct 
cause of the experiment here discussed. This experiment was carried 
out in the Gulf Coast region near Houston, Texas. 

On the assumption that Gulf Coast gumbos may be treated as 
fluids, within limits, calculation of theoretical velocities to be ob- 
tained from earth-air mixtures in varying proportions were carried 
out from the equations given in the article mentioned,” namely: 


{XEs + (1 — X)E:}{XPi + (x — 


Where: 


P,=density of air=.oo12 
E,=elasticity of air=1.2X 10° 
P.=density of earth, assumed as 1.9 
E,=elasticity of earth, computed as 5.58X 10" 
X =proportion air to total by volume 
1— X =proportion earth to total by volume 
V =velocity of sound in mixture 


If the values of the composite velocities be plotted against pro- 
portions of air to earth, a curve of the form in Figure 1 is obtained. 

Though this curve should perhaps be considered only qualitatively 
because of the assumption of earth as a fluid, it indicates the prob- 
ability of obtaining velocities from earth-air mixtures, whose magni- 
tudes may be actually less than that of sound in air alone. This pre- 
diction has been verified experimentally by shooting a short detailed 
refraction profile, extreme precautions being taken to muffle the shot 
point and hence eliminate any effect of the air wave. The recording 
was begun at a distance of 5 feet from the shot point and carried out 
on the surface to a distance which definitely indicated the ‘“un- 
weathered” velocity on the time-distance curve. This curve is shown 
as curve (A) in Figure 2. 


1 A. B. Wood, A Textbook on Sound, pp. 327-28. 
2 A. B. Wood, of. cit. 


SEISMIC WEATHERED OR AERATED SURFACE LAYER 1233 


As can be seen from the curve, a velocity within the weathered 
layer is indicated (slope of curve near origin) as low as 550 feet per 
second, while the velocity of sound in air alone is approximately 1,100 
feet per second. Assuming the low velocity to be due to the inclusion 


of air as postulated, we would naturally expect the velocity to be 
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lowest nearest the surface, where the proportion of air should be 
greatest. Referring this velocity (550 feet per second) to the curve 
(Fig. 1), the corresponding air concentration is 0.001 part by volume, 
which seems probable in the near-surface materials. A computation 
of the depth corresponding with the “‘break-point”’ of the time-distance 
curve indicates 714 feet as the thickness of the “aerated layer.” 
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After completion of this profile, a hole was dug at the shot point 
and the top of the water table was determined to be at a depth of 
slightly less than 8 feet. On the assumption that below the water table 
all air in a free state is displaced by water, whose density and elasticity 
are of the same order of magnitude as the earth under consideration, 
and hence has no material effect on the velocity of transmission of 
sound through such substance, both the shot point and recorders 
were lowered to the top of the water table. At this depth the shorter 
or “weathered” portion of the profile was re-shot with the result 
shown in curve (B), Figure 2. 
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On examination, this curve indicates an “‘unweathered” velocity 
(5,200 feet per second), which passes through the origin with only that 
curvature which is to be expected with a slight increase of velocity 
with depth. At the same recording distances from the shot point, the 
“weathered” portion of the time-distance curve has been eliminated 
with the elimination of the “aerated layer.’’ The obvious conclusion 
is that where the water table is sufficiently shallow the thickness of 
the aerated layer is determined by its position. However, if the water 
table be at considerable depth, the ‘‘weathered”’ thickness is deter- 
mined by the depth to which any appreciable amount of air penetrates. 

This explains, in a measure, the sometimes extreme variation of 
“weathered” thickness with different surface materials. It would also 
predict a change of weathered thickness with wet and dry seasons. 

Though these experiments have not been carried out to such an 
extent as to make the conclusions indisputable, it may be said that 
all the experimental data obtained favor the theoretical predictions 
resulting from the assumption of an “‘aerated”’ layer. 
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ACCURACY OF DETERMINATION OF RELATIVE 
GRAVITY BY TORSION BALANCE! 


DONALD C. BARTON? 
Houston, Texas 


ABSTRACT 

Determination of relative gravity by the Eétvés torsion balance has been shown by 
Oltay and others to be as accurate as the determination by the invariable pendulum. 
In the present paper, the probable error of the torsion balance determination of rela- 
tive gravity is calculated from the error of closure of 45 traverses comprising 2,800 
stations, most of them in the Gulf Coast region. The probable error of each individual 
observation in those surveys is calculated to have been +2.2 Eétvés units for the 
traverses taken together, or from +1.9 E for the traverses with the larger station 
interval to + 5 E for the traverses for the small station interval. But the greater error 
of the individual observations in the traverses which have the shorter station interval 
has been compensated, intentionally, by that shorter interval, and the probable error 
of the determination of relative gravity by those 45 traverses is approximately 0.4 
milledyne per 10 kilometers of traverse without regard to the magnitude of the probable 
error of the individual observations or of the station interval. The probable error of the 
determination of relative gravity between key points in good torsion balance surveys 
presumably is approximately + 2.5 to + 3.5 E per 10 kilometers airline distance between 
the two places. If pendulum determinations of relative gravity are used to supplement 
and to increase the accuracy of torsion balance surveys, there is a minimum interval 
at which the pendulum stations should be placed, for at lesser intervals the determina- 
tion of relative gravity by the torsion balance is more accurate than that by the pendu- 
lum. That minimum interval ranges from 100 to 200 kilometers for the pendulum ob- 
servations of the first quarter of the century to 8 to 50 kilometers for first class modern 
pendulum observations. 


The relative accuracy of the determination of relative gravity by 
good torsion balance surveys in favorable terrane has been shown 
by Oltay, Jung, Numerov, and others to be comparable with the 
accuracy of the determination of relative gravity by the invariable 
pendulum during the past quarter century. The differences between 
the relative gravity by the torsion balance and pendulum on six 
traverses in the Arad district of the Hungarian plain were respectively 
0.001, 0.001, 0.004, 0.000, 0.005, 0.000 cm. sec.-? The respective 
lengths of the traverses were 12, 9, 13, 20, 42, and 50 kilometers.* 
The probable error of the pendulum observations was calculated to 


1 Read before the Association at the Oklahoma City meeting, March 25, 1932. 
? Consulting geologist and geophysicist, Petroleum Building. 


3 Karl Oltay, “Die Genauigkeit der mit der Eétvésschen Drehwage Durchge- 
fuhrten Relativen Schwerkrafts-messungen,”’ Geoddtische Arbeiten d. Baron R. 
Edtvisschen Geophysischen Forschungen, No. III (Budapest, 1928). 
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be +0.001 cm. sec.~? In the surveys in the Ries area of Bavaria, the 
differences between the determination of relative gravity by the 
torsion balance and the pendulum for six different lines were 0.003, 
0.000, 0.000, 0.002, 0.003, 0.001 cm. sec.~* The lines ranged from 4 to 
6 kilometers in length.' In the Ural Emba district of Russia, the 
differences were 0.001, 0.002, 0.007, 0.004, 0.003, cm. sec.~? The length 
of the traverses ranged from 9.5 to 15 kilometers.’ In the Grosny dis- 
trict of Russia, the differences were 0.001, 0.002, cm. sec.~* and the 
lengths of the traverses ranged from 9.5 to 19 kilometers. As the prob- 
able error of the pendulum observations was +0.001 to +0.002, and 
as errors of three times the probable error should be expected to be 
fairly common, the torsion balance surveys gave as accurate a de- 
termination of relative gravity as did the pendulum observations. 
But from the data, it is impossible to tell whether the determinations 
of relative gravity by those torsion balance surveys were more ac- 
curate than those by the pendulum. 

The accuracy of the torsion balance determinations of the gradient 
has been checked by the writer by calculation of the normal north- 
ward (planetary) gradient from torsion balance observations and com- 
parison of that value with the theoretical value. The theoretical value 
was 7.14 £ north in comparison with the value of 7.4 E., N. 34°E., 
which was calculated from the torsion balance observations.* 

The present paper reports the determination of the effective ac- 
curacy of the (horizontal) gradient which is observed with the Eétvés 
torsion balance and of the accuracy of relative gravity which is cal- 
culated from the gradient. 

The data on which the study is based consist of the error of closure 
in relative gravity in forty-five closed traverses from commercial 
torsion balance surveys. The length of the traverses ranged from 3 
to 144 kilometers and in twenty-eight of the traverses, the lengths 
were evenly distributed over the range of 10 to 40 kilometers. The 
interval between stations on thirty-two of the traverses was 300-500 
meters, on five traverses, 150+ meters, and on eight traverses, 600-— 
800 meters. The instruments which were used on the surveys mainly 

1 Karl Jung, “Drehwagemessungen im Ries bei Nordlingen,”’ Zeitschrift fiir Geo- 
physik, Vol. 7 (1931), pp. 18-19. 

2 B. Numerov, ‘Results of the General Gravity Survey in the Emba District,” 
Zeitschrift fiir Geophysik, Vol. 5 (1929), pp. 268-70; and “Results of Gravitational Ob- 
servations in the Region of Grosny in 1928,” ibid., Vol. 5 (1929), pp. 271-75. 

* Donald C. Barton, “Gravity Measurements with the Eétvés Torsion Balance,” 


in “Physics of the Earth,” Vol. II, ‘““The Figure of the Earth,” Bull. National Research 
Council, No. 78 (February, 1931), pp. 186-87. 
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were the small type of the Siiss “Original Eétvés” torsion balances 
and, on a few traverses, large Askania (automatic) instruments. The 
calculation of relative gravity around each traverse was made by a 
method in which superposition of a transparent graph over a gradient 
arrow gives the spacing of the isogams at that station; the isogams 
are sketched forward toward the next station and are bent as the 
experience of the geophysicist and the exigencies of the adjacent 
gradient arrows seem to demand. Single abnormal gradient arrows 
were neglected or given reduced weight according to subjective esti- 
mation of their value. The location of most of the surveys was in the 
Texas-Louisiana Gulf Coast within 150 miles of the coast, but a few 
of the surveys were farther inland and some of the surveys were 
located outside of Texas and Louisiana. Most of the terrane was good, 
but a few of the surveys were made in areas in which caliche was 
present in the subsoil and other surveys were made in river bottoms. 
The magnitude of the gradient most commonly was less than 20 E. 
The station sites on a few of the surveys were prepared with the 
meticulous care of the European custom, but on most of the surveys, 
slightly less care was used. Probably two-thirds of the observations 
were made during the day-time and one-third at night. The general 
rate of occupation of stations was two stations per day per instrument 
and two stations per night per instrument if observations were taken 
at night, although in about a quarter of the surveys, the rate of ob- 
servation was one station during the day and two at night per in- 
strument. These surveys were made in commercial oil prospecting in 
which, generally, speed and moderate accuracy are regarded as more 
important than high scientific accuracy, and the purpose of the sur- 
veys was qualitative or only roughly quantitative rather than ac- 
curately quantitative. The surveys were made by many different 
torsion balance parties. Most of the observers were not highly trained 
scientists. 

The precision of scientific measurements of such quantities as 
gravity in America is stated in terms of what is called the “probable 
error’; in Continental Europe, the “mean square error” is used; the 
“probable error” is 0.6745 times the “mean square error.”’ The term 
“probable error” is a technical term and is defined as the value such 
that half of the errors which are likely to occur in a large number of 
measurements of the quantity will be larger and half smaller; that is, 
if the determination of gravity at a station is stated in an American 
publication to be 978.569 + 0.001 cm. sec.~*, it is understood that there 
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are 50 chances out of 100 that algebraically, the actual error is less 
than —o.oo1 or greater than +0.001 and that there are 50 chances 
out of 100 that algebraically, the actual error will be between —o.oo1 
and +0.001. According to the theory of the “‘probable error,” the 
actual error will be five times the “‘probable error” 1 time in 1,000, 
four times the “‘probable error” 7 times in 1,000, three times the 
“probable error’ 43 times in 1,000, twice the “probable error” 177 
times in 1,000. The “probable error,” therefore, affords a method of 
comparing the precision of different measurements and of determin- 
ing the probable range within which the actual error of a measure- 
ment of quantity may lie. 

The “probable error’ of a series of measurements depends upon 
the “probable error’ of each measurement, the number of measure- 
ments, and the character of the series of measurements. If PE is the 
“probable error’”’ of the series, if pe—1, pe—2-- - pe—n are the re- 
spective “probable errors’’ of m measurements numbered, 1, 2, - - - to 
n and if the series consists of repeated measurements of a single 
quantity, the “probable error’’ of the arithmetic mean of the measure- 
ments is: 


(1) PE = 


Vm 


if the pe’s are all the same. But if the series consists of a series of con- 
secutive applications of the measuring device, as for example, the 
consecutive use of a surveyor’s tape, in measuring the length of a line 
which is longer than the tape: 


(2) PE = (per)? + (pes)? + (pen)? 
or if the pe’s are all the same: 
(3) PE = pe 


The “‘probable error” (pe) of a mean in terms of the actual error 
(ae) of each measurement is given by the formula: 


(4) PE == 0.67454/ + + (den) 


n 


The accuracy of torsion balance surveys can be expressed by a 
probable error which can be calculated by the formula: 


1 Mansfield Merriman, A Textbook on the Method of Least Squares, 8th ed. (1911), 
pp. 66-79. 
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EC,\? EC,\? 
(5) PEus: = 0.6745 


n 


where: 

EC, EC2, - - - EC, are the respective errors of closure in Ag of n 
closed torsion balance traverses numbered from 1, 2, - - - ton; 

Ni, No, ---N, are the respective number of stations in those 
traverses; 

I, L2,---L, are the respective lengths of those traverses in 
centimeters; 

n is the number of traverses; and 

PEv;: is the probable error in the dimensions of the gradient. 

The formula is adapted from Merriman’s formula:* 


2 d,? 
6 
(6) +=) 


where PE is the probable error of each measurement with a surveyor’s 
chain. It is assumed that the length of () lines (1, 2, - - - m) has been 
measured with the chain and is expressed in chain lengths; and that 
two measurements of each line were made. The respective differ- 
ences between those duplicate measurements of each line are repre- 
sented by d:, do, - - - d, and the length of the respective lines in chain 
lengths by LZ, Le, : - - L,. PE, then, is the probable error of any one 
of the many single measurements of the chain. The probable error 
of the measurement of each line, then, can be calculated by formula 
(3) from the PE which has been obtained by formula (6), and from 
the number of individual consecutive measurements with the chain 
which were necessary to measure the length of the line. 

The measurement of Ag by the torsion balance around a series of 
closed traverses is comparable to that duplicate measurement of each 
of a series of lines with a surveyor’s chain. Each duplicate measure- 
ment of the length of the line may be regarded as a closed traverse 
which will consist of as many single measurements with the chain 
as the length of the line in chain lengths. A torsion balance traverse 
which is composed of NW successive stations may be regarded as NV 
consecutive measurements of a line comparable with those consecu- 
tive applications of the surveyor’s chain to obtain the length of those 
lines. The error of closure of Ag, (EC), of a closed torsion balance 

1 Op. cit., p. 103. 
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traverse is comparable with the difference, (d), between the duplicate 
measurements of the length of a line with the surveyor’s chain; and 
the number of torsion balance stations, NV, is comparable with the 
number of applications of the surveyor’s chain. 

The formula (6), therefore, may be re-written in the notation 
which we have assumed for the torsion balance surveys, to read: 


(7) = 0.6 4/2 —) 
eee 
145 n\ N, Nz 


The error of closure, EC, is in terms of Ag (that is, cm. sec.~*) and the 
probable error (PE), therefore, would be given in terms of Ag per 
unit of measurement which according to our preceding assumptions 
is the station interval. The latter rather commonly is fairly constant 
within a traverse, but is variable from traverse to traverse. The cal- 
culation of Ag over a station interval has the mathematical form of 
a gradient times distance, that is, times the length of the station 
interval; and, therefore, any error in Ag over a station interval will 
tend to be lineally proportional to the length of the station interval. 
The calculated “probable error” in terms of Ag and of the inconstant 
station interval would be useless. The “probable error,” however, 
may be reduced to the dimensions of a gradient and then will be in- 
dependent of the length of the station interval. All error in the cal- 
culation of Ag over a station interval may be assumed to be in the 
gradient, if, as we usually assume, the stations are close enough to- 
gether so that the variation of the gradient is linear from station to 
station. The probable error in Ag over a station interval, then, will 
vary linearly with the length of a station interval and may be ex- 
pressed as the product of the length of the station interval times that 
probable error in the gradient. Conversely, the probable error of the 
gradient may be expressed as the quotient of the probable error in 
terms of Ag over the station interval divided by the length of the 
station interval. If the terms 


EC? 
N 


of formula (7) are weighted in terms of the square of the reciprocal 
of the length of the station interval, S7, to read 


EC EC, 
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the errors of the torsion balance surveys, then, are expressed in terms 
of gradient and not of Ag and the probable error of the torsion balance 
measurements will be in terms of gradient. But 
L EC? 
SI =— and the other t 
yan other terms (sD? 
EC? N-EC? 
or 
L? L? 
N 


can be written: 


N?3 
EC? 
If the terms Y of formula (7) are replaced by 


= = formula (5) is 
obtained, which gives the probable error in the torsion balance ob- 
servations in terms of gradient; that is, in Eétvés units, although the 
observed errors are in terms of Ag; that is, of cm. sec. 

Formula (5) involves the two assumptions, first, that the error in 
calculation of Ag over a station interval is a linear function of the 
length of the station interval and of the error in the gradient and, 
second, that the error is independent of the magnitude of the gradient. 
Neither assumption is exactly true but in the present study no at- 
tempt has been made to analyze that more complicated variation of 
the error of torsion balance measurements. 

The original data and the calculated values of 

EC?-N 
are given in Table I. 

The “‘probable error’”’ of the determination of the gradient at those 
twenty-eight hundred odd stations, therefore, is of the general mag- 
nitude of +2.2 EZ, but the “probable error” of the determination of 
the gradient may be different in different areas. If the surveys of 
Table I are grouped by the length of the interval between stations, 
the “probable error” of the gradient varies considerably between 
groups. The various “probable errors” are given in Table II. They 
range from +1.5 E to +5.0 E. The corresponding “probable error”’ 
in the individual traverses must have a yet higher range. The figure, 
+2.2 E, therefore, itself has a large “probable error’ and only its 
general magnitude is of significance. Intuitively, the writer for years 
has estimated the common error of torsion balance observations of 


the gradient in the Gulf Coast region as of the magnitude of 1.5 to 
2.5 E. 
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TABLE I 
EC*N 
ORIGINAL DATA AND CALCULATED VALUES OF 
I? 
Error of Number Lengthof EC*N | Error of Number Lengthof EC*N 
Closure EC of Traverse Closure EC of Traverse 
in 10-*Cm. Stations in Stations 
Sec. N in Km, Sec. N i 
71 2. 
72 2. 
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TABLE II 


PROBABLE ERROR OF THE DETERMINATION OF THE GRADIENT IN SURVEYS WITH 
DIFFERENT STATION INTERVALS 


Station intervalin meters.... 150 300 400 500 600 700 800 
Number of traverses 5 12 8 12 6 2 I 
Probable error 35-4 22.8 £1.96 


TABLE III 


PROBABLE ERROR OF DETERMINATION OF Ag BETWEEN Paces BY Goop TorsION 
BALANCE SURVEYS IN 10~* DyNEs (0.1 MILLEDYNE) 


Distance between places in kilometers. 1 5 10 20 
Probable error of: single direct trav- 

erse; or single traverse distance 

measured along traverse; or two 

more or less indirect traverses 
Several traverses 
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11.4 37 19.5 12.6 
7-6 44 30-5 2.7 
11.4 45 36.0 5.8 201 105. 
11.4 51 29.6 7.6 75 41. 
7.6 56 28.7 4.0 102 23. 
7.6 64 36.6 2.9 97 21. 
° 123 37.2 ° 
3-8 65 19.5 0.2 
7-6 82 23.2 8.9 
7.6 73 22.6 8.3 
13.3 201 69 7 
36.1 164 65 5° 
30.1 171 53 76 
5-7 64 27 2 
4.8 5° I 6.4 I 
8.0 79 5 I 5 
4.8 26 13 
1.6 35 3 7 
1.6 38 4-7 
4.8 71 18.5 
2E 
40 80 100 
26 
= = = = +6 +8 +9 
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That “probable error,” +2.2 E, is the composite error of a con- 
siderable number of factors. It arises: (a) from all instrumental errors 
in the determination of the crude gradient at the station, (b) from all 
errors in applying corrections to the crude gradient, (c) from errors 
in plotting gradient arrows on the map and from the crudeness of 
gradient arrows on a scale of 1 mm.=£ as the starting point for the 
calculation of relative gravity, (d) from inaccuracy in the use of the 
writer’s method of calculating relative gravity, and (e) from error 
arising from the fact that the interval between stations is not suffi- 
ciently close to give a true picture of the variation of the gradient. 
It probably also varies with the magnitude of the gradient; the data 
are not sufficient to study that variation. An actual error five times 
the ‘“‘probable error” should occur once in one thousand times, ac- 
cording to the theory of the “probable error.’ Divergent gradients 
which depart vectorially much more than 5X + 2.2 E are moderately 
common. Such plainly aberrant gradient arrows were disregarded by 
the geophysicist in the calculation of relative gravity in the traverses 
from which the data of Table I were taken, and their weight, there- 
fore, was not felt in the calculation of that “probable error.” 

The ‘‘probable error” of the determination of relative gravity by 
good torsion balance surveys can be stated in a simpler and more 
practical form. The “probable error” of the determination of relative 
gravity by the torsion balance can be controlled within certain limits 
by variation of the station interval. Formula (3) can be rewritten: 
from: PE = V/n- pe 


(8) 


L 
to: PEs, = V/n-SI- peve: =—peves 
Vn 


for the special case of the torsion balance. If the number of stations 
in a traverse of constant length is varied as the square of the variation 
of the “probable error’”’ of the individual observations, the “‘probable 
error” of the determination of Ag for the traverse remains constant. 
In practice in good torsion balance surveys, the station interval be- 
tween stations is decreased in areas of irregular and erratic gradient 
arrows and is increased in areas of consistent gradient arrows. The 
“probable error” of good ordinary torsion balance surveys in practice, 
therefore, should be independent of the station interval and should 
vary as the square root of the length of the traverse. If (SI: pevs:) 
in formula (8) is maintained constant, then (m) will vary directly 
with the length of the traverse, and PE,, will vary as the square root 


/ 
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of the length of the traverse. As SI: pey,, somewhat commonly in 
practice is held approximately constant, the “‘probable error” of the 
determination of Ag by ordinary good torsion balance surveys can 
be expressed in terms of Ag per some standard length of traverse. 

That this is so, can be seen by analysis of the error of closure of 
those 45 traverses. The error of closure of each closure of Ag in each 
of those traverses can be scaled up or down by the formula 


(9) EC = = 


to the error of closure (EC\o), which the traverse would have had, if 
its length had been ro kilometers instead of (LZ) kilometers, and if 
there had been no change in the station interval. The resulting errors 
of closure per 10-kilometer traverses have been sorted by length of 
the station interval and are given in Table IV. The respective errors 
of closure and the respective “probable errors” of Ag for the traverses 
of 100, 200, 300, 400, 500, 600 meters station interval, can be seen to 


TABLE IV 
ERRORS OF CLOSURE SCALED TO A 10 KM. LENGTH OF THE TRAVERSE 


Station interval in meters 100 400 600 


8 


Error of closure of Ag 
in 10~* dynes 


PAW ACO HSU 
ON 
002000 Ow 


“Probable error’ of 
$3.6 23.4 24.6 +6. 


+ 
w 


“Probable error” for all 45 traverses: + 3.9 


be independent of the length of the station interval. That “probable 
error” ranges only from 0.32 milledyne for the traverses with a station 
interval of approximately 600 meters to 0.46 milledyne for the traverses 
of approximately 300 meters; and for the traverses of approximately 
100, 200, 500, and 600 meters length, that “‘probable error’ ranges 
only from 0.32 to 0.36 milledyne. The “probable error” of the indi- 
vidual observations, however, was considerably larger in the traverses 


/ / 
Fo 8.7 6.8 15.6 14.1 4-3 6.0 

4-9 6.8 5.0 5-0 2% 
3-4 43 4.7 
2.3. 3-4 41 3-4 
2.0 3.2 1.9 
3.7 1.6 
0.0 
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with the shorter station interval (Table II). That larger ‘“‘probable 
error” of the individual observations in the traverses with the shorter 
station interval has been compensated by the use of that shorter 
station interval; and the “probable error” of the traverse as a whole 
is no greater for those traverses with the shorter station interval than 
for those with the longer station interval. The “probable error’’ of 
determination of relative gravity by those 45 traverses is 0.4 milledyne 
per 1o kilometers of traverse without regard to the station interval; 
and the actual error of closure, with two exceptions, is less than 0.9 
milledyne and with only five exceptions, is less than 0.75 milledyne 
per to kilometers of traverse without regard to the station interval. 

The “probable error’ of the determination of relative gravity 
between two places by torsion balance surveys depends in part on 
the number of torsion balance traverses which connect the two places, 
decreasing inversely as the square of the number of equally good 
traverses. Key places in torsion balance surveys in general are con- 
nected by at least two traverses, and in many surveys are connected 
by a net of traverses. Different weights in general should be given 
the different traverses in the calculation of the relative gravity be- 
tween the two places and in the calculation of the “probable error’ 
of that determination of the relative gravity. But that “probable 
error” in practice probably will be 1/+/2 to 1/+/3 times the “probable 
error” of the determination of the relative gravity by a single trav- 
erse. The “probable error” of the determination of relative gravity 
between key places within the surveys of those 45 traverses, therefore, 
should be (1/+/2 to 1/+/3) times 0.4 milledyne or 0.30 to 0.25 mille- 
dyne for places 10 kilometers apart along the route of the torsion 
balance traverses. For an airline distance of 10 kilometers between the 
two places, the corresponding “probable error” in practice should 
be 0.7 to 1.0 times the “probable error’’ of the determination of Ag 
by an airline torsion balance traverse between the two places and, 
therefore, should be 0.3 to 0.4 milledyne. 

Shrewd adjustment of errors in those 45 traverses presumably 
has reduced the “probable error” of those determinations of relative 
gravity below the figures of the preceding paragraph. The error of 
closure in Ag in torsion balance traverses commonly is caused not 
only by the cumulative effects of small, or moderately small, errors 
at each set-up but also by a few large errors. A skilled interpreter 
commonly can recognize the weak spots in which those errors are 
mostly likely to have come in; by disregarding or adversely weighting 


/ 
/ 
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the effect of single seriously abnormal gradient arrows, or less com- 
monly, pairs of gradient arrows, he attempts to eliminate such errors 
in the first calculation of relative gravity; in the subsequent ad- 
justment of the traverse or traverses, somewhat commonly much 
of the error of closure may be eliminated by a slight change in his 
weighting of those few irregular gradient arrows. The variation of 
the gradient from station to station in many places is not linear and 
the geophysicist has a certain latitude in his subjective estimation 
of the probable variation of gravity between the two stations. The 
error of closure in many traverses can be almost eliminated by taking 
those slightly different alternative choices in the spacing and the run- 
ning of the isogams. The subjective weighting of observations, of 
course, is open to danger of error and in the hands of an inexpert 
geophysicist, leads to much more erroneous results than the auto- 
matic mathematical distribution and adjustment of the error. But 
the writer believes that the expert interpreter can obtain better 
accuracy by exercising subjective judgment in throwing as much of 
the errors of closure as possible into the weak places in lines rather 
than by automatic mathematical distribution of the errors. The 
“probable error” of the determination of relative gravity in those 45 
traverses presumably, therefore, is slightly less than +0.25 mille- 
dyne per 10 kilometers of double traverse and less than +0.35 mille- 
dyne between places 1o kilometers apart, airline distance. 

That “probable error,” 0.25 milledyne per 10 kilometers of double 
traverse or 0.35 milledyne per 10 kilometers of airline distance, pre- 
sumably is characteristic of most good torsion balance surveys in 
what a torsion balance operator would call good, fair, and slightly 
poor torsion balance terrane. Those 45 traverses lay in all types of 
terrane in the Gulf Coastal Plain of Texas and Louisiana and a few 
of them were on the llanos of Venezuela and elsewhere than in the 
Gulf Coastal Plain. Less than one-fourth of the traverses were in the 
Coastal Prairies of Texas and Louisiana, which are especially favor- 
able for a low “probable error” of observation. The 45 traverses, 
therefore, are fairly characteristic of torsion balance work in areas 
which are commonly regarded as good, fair, or slightly poor for 
torsion balance observations. Mathematically the same grade of 
accuracy could be maintained by use of a sufficiently close net of 
stations to compensate the “probable error” in the individual ob- 
servations, but the number of stations necessary increases as the 
square of the increase in the “probable error” of the individual ob- 
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servations, and, therefore, the cost and tediousness of too close a net 
of stations puts a practical limit to the maintenance of that accuracy 
of +0.35 milledyne between places 1o kilometers apart, airline dis- 
tance, or for the maintenance of any particular grade of accuracy. In 
reconnaissance in oil work, a station interval less than 250 meters 
commonly is impracticable, and if the “probable error” of the in- 
dividual gradient observations rises above +2 E, it becomes imprac- 
ticable to maintain an accuracy of +0.35 milledyne in the determina- 
tion of relative gravity between places approximately 10 kilometers 
apart, airline distance. In some areas in South Texas, in which the 
caliche lies close to the surface, the “probable error” of individual ob- 
servations probably is greater than +7 E; and in certain areas in 
East Texas, the ‘probable error” of individual observations probably 
is greater than +5 E. The “probable error” of the determination of 
gravity between two places approximately ro kilometers apart in 
such areas will be very much greater than +0.35 milledyne. 
Knowledge of the “probable error” of the determination of Ag 
by torsion balance surveys is important in certain types of combina- 
tions of torsion balance and pendulum surveys. According to a rather 
common suggestion from “‘pure” geophysicists, pendulum observa- 
tions of relative gravity at a net of key stations should be used to 
provide a framework of accurate gravity benchmarks to which the 
supposedly less accurate torsion balance surveys would be tied and 
adjusted. But the accuracy of the torsion balance surveys will be in- 
creased thereby only if the “probable error’ of the determination of 
relative gravity by the pendulum between two stations is less than 
that by the torsion balance surveys. The determination of relative 
gravity by the relative pendulum is independent of the distance be- 
tween stations; the value of Ag between two stations is simply the 
difference between the observed values of relative gravity at the two 
stations, and the “probable error” of Ag between the two stations is 
the square root of the sum of the squares of the respective “‘probable 
errors” of the individual observations. The “‘probable error’ of the 
determination of Ag between two places by good torsion balance 
surveys varies, in practice, with the square root of the distance be- 
tween the two places. The torsion balance surveys will measure Ag 
more accurately than the pendulum if the two stations can be con- 
nected by very short torsion balance traverses, and the pendulum 
will measure Ag more accurately than the torsion balance if the two 
stations can be connected only by very long torsion balance traverses. 
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If the pendulum stations are placed closer together than a certain 
minimum distance, the ‘probable error” or the pendulum determina- 
tion of Ag between adjacent stations will be greater than the “‘prob- 
able error” of the torsion balance determination of that Ag; and ad- 
justment of the torsion balance survey to the pendulum values of Ag 
will decrease the accuracy of the torsion balance survey. If the dis- 
tance between the adjacent pendulum stations is greater than that 
minimum distance, the “probable error’ of the pendulum determina- 
tion of Ag between adjacent stations will be less than that of the 
torsion balance determination; and the accuracy of the torsion balance 
survey will be improved by adjusting its results to the pendulum 
results. 

The accuracy of the determination of relative gravity by the rela- 
tive pendulum varies. The “probable error”’ of the pendulum observa- 
tions of the latter part of the past century presumably was, in general, 
many milledynes. The probable error of the observations of the re- 
cent past in general has been +1 milledyne, although a few poorer 
but acceptable stations have had a “probable error” of +2 mille- 
dynes. The better modern observations are reported to have a “prob- 
able error” of +o.3 to to.5 milledyne. The geophysical department 
of one oil company asserts that its pendulums are reoccupying sta- 
tions with a “probable error” of observation less than +0.1 mille- 
dyne. The “‘probable error’’ of the determination of the elevation or 
relative elevation of the pendulum station is relatively unimportant 
in connection with the pendulum observations whose “probable error” 
is +1 milledyne or more, but is, of course, of very considerable im- 
portance in connection with observations whose “probable error’’ 
is less than +0.3 milledyne, for the variation of (g) with elevation 
is +0.308 milledyne per meter. The “probable error” in feet of first 
class (but not precise) levelling is approximately +0.05 (number of 
miles), but only 1.5 miles of closed traverse, or 3 miles of unrepeated 
traverse can be run per day in levelling of that accuracy and, there- 
fore, most commercial levelling in connection with commercial geo- 
physical work probably will be speedier and less accurate and its 
“probable error” may have to be taken seriously into account. The 
“probable error” of the station determination of gravity will be the 
square root of the sum of the square of the “probable error” of the 
pendulum determination of gravity, plus the square of the “probable 
error” of the determination of the elevation. The “probable error” of 
the determination of Ag between two places will be the sum of the 
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squares of the respective “probable errors” of the two pendulum 
observations and of the two determinations of the elevation. The 
“probable error’ of the pendulum determination of gravity at each 
key station can be decreased by reoccupation of that station or by 
the occupation of additional closely adjacent stations which are tied 
together by torsion balance determinations of the value of Ag be- 
tween the key station and each of the subsidiary stations. The de- 
crease in the “probable error” will vary as the reciprocal of the square 
root of the number of stations. 

The minimum interval between pendulum observations which 
are to supplement ordinary, good torsion balance surveys in good, 
fair, or even'slightly poor torsion balance terrane is given in Table V. 


TABLE V 


MINIMUM INTERVAL AT WHICH PENDULUM OBSERVATIONS SHOULD BE USED TO TORSION 
BALANCE SURVEYS 


“Probable Error” “Probable Error” of Pendulum Determination of Ag 
of Torsion Bal- in 1X10~* Dynes 
lance Determina- 
tion of Agin ro~* 
Dynes for Places + +s +3 
10 Km. A part 
Airline Distance Interval in Kilometers 


A* B d A A 


1 
2 


Poorer ‘“‘Good” 
surveys...... 7 
Average “Good” 
surveys 
Better “Good” 
2 650 120 30 55 


* A. The two places are connected by a single direct or two indirect torsion balance traverses. 
B. The two places are connected by several traverses. 
In surveys primarily by the pendulum, or in pendulum surveys which 
are supplemented by sketchy torsion-balance surveys, the intervals 
between the pendulum stations very much less than the minimum 
intervals of Table IV may be used in order to increase the accuracy 
of the survey. 
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CURVATURE OF EQUIPOTENTIAL SURFACES! 


M. M. SLOTNICK? 
Houston, Texas 


ABSTRACT 


If one is to have a complete and satisfying understanding of the theory of the tor- 
sion balance, it is essential to have clearly in mind a physical picture of the quantities 
involved. Most authors assume either that their readers are well versed in the differen- 
tial geometry of surfaces in a Euclidean three-dimensional space or that a discussion of 
surface curvature is far beyond the reader’s ability to grasp. Too often is neither of 
these true. The writer attempts to present a mathematical discussion of this matter in 
which only the bare fundamental concepts of the differential calculus are needed. 


INTRODUCTION 

In this paper, the writer attempts to develop the mathematical 
theory of the curvature of the equipotential surfaces due to the 
gravitational field, with the purpose of producing a clear “physical’’ 
picture of this quantity. It is the writer’s. impression that such a de- 
velopment, although fundamental to the mathematician, is lacking 
in the usual torsion balance literature available to the geophysicist. 
Our goal is to find the curvature relations used in torsion balance 
work. The mathematical needs for this subject are not too great— 
and for further simplification, the demands of mathematical rigor are, 
at times, sacrificed in this paper. 


DEFINITIONS 


Consider an arbitrary point P on a surface S, which is “smooth” 
in the neighborhood of P in the sense that the equation of S may be 
expressed as a series when the surrounding space is referred to a 
convenient cartesian codrdinate system. The tangent plane to S at 
P is that plane in which all the lines tangent to S at P lie. The straight 
line perpendicular to this plane at P is the normal to the surface S at 
that point. 

COORDINATE SYSTEM 


We now proceed to choose, as is our privilege, a codrdinate system 
in such a manner that the mathematical work involved is materially 


1 Read before the Association at the Oklahoma City meeting, March 25, 1932. 


2 Geophysics department, Humble Oil and Refining Company. Introduced by 
L. W. Blau. 
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simplified. Let the normal to S at P, oriented in either of its direc- 
tions, be the z axis. In the tangent plane to S at P let any two mutu- 
ally perpendicular lines be chosen as the x and y axes. Both of these 
axes are obviously perpendicular to the z axis. Let us also point out 


Fic. 1 


that the tangent line at P to any curve on S passing through P lies 
in the (x, y) plane. The situation thus far described is schematically 
shown in Figure 1, where, for convenience, only part of the surface 
S in the neighborhood of P is indicated. 


EQUATION OF S 


Suppose that the equation of S, referred to the chosen coérdinate 
system, is 
(1) 2 = 2(x, 9), 


and that this expression is expanded in series about P; that is, 
(2) 2=at (px + qy) + W(rx? + asxy + 


where 
=) =) 
Ox Oy 


r=—, s=—; t=- 
Ox? Oxdy dy? 


the values of these derivatives being taken at the point P, whose 
coérdinates are, of course, (0, 0, o). The series (2) has been written 
only to the second order terms, which are all that are required. 

Since the surface S passes through P:(0, 0, o), equation (2) must 
be satisfied by its coérdinates. Accordingly, a=o. Again, the slope 
of the curve on S at P, which is the section of the surface made by 


R 
| 
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the (x, z) plane, is the value of p there, and consequently this is zero. 
Similarly, g=o. Thus, by virtue of our choice of axes, the series (2) 
reduces to 


(3) z= (rx? + 2sxy ty?) 4+ .-.-, 


CURVATURE OF A PLANE SECTION OF S AT P 
Let C be the curve of intersection of any plane zPR through the 
z axis with the surface S, and let the angle between this plane and the 
(x, z) plane be A. Obviously, then, the curve C is a plane curve and 
PR, the tangent line to it at P, lies in the (x, y) plane, with the angle 
xPR equal to X. The lines Pz and PR are at right angles to each other, 
and we choose these lines as céordinate axes in their plane zPR. Since 


(4) x = Rcosh, y = Rsink, 
the equation of C in terms of the codrdinates (R, z) is, by virtue of (3): 
(5) z = + 25 sind +ésin?A)+---. 


We recall at this stage the theorem of the elementary differential 
calculus that the curvature of the curve whose equation is y=/(x) 
at a point is the value of 


(6) k= 


at that point. In equation (5), z is expressed as a function of R, and, 
at P, z=R=o. The quantities 7, s, # are constants, inasmuch as they 
represent the values of the second derivatives of z with respect to x 
and y at P; and ) is a constant for the particular section C chosen. 
Consequently, the curvature of C at P is the value of 


d*z 


dR? 
(7) k = ’ 


dz 2 3/2 
when z= R=o. 


Differentiating (5), we obtain 


dz 
(8) Riv + as sind cosd + ésin’d) + 


d*y 
dx? 
dy 2 3/2 
|| 
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and 
d*z 


The terms of each of these series beyond those written involve R to 
at least the first power. At P, then, where z= R=o, 
dz d*z 


10.) 
ap 


= rcos?*’A + 2ssinX cosA + ésin?X. 


It is these values, as we have seen, that we must use in (7) to ob- 
tain the curvature of C at P, and substituting these, we have 


(11) k = rcos?\ + 2ssin\ + ésin? X. 


In ordinary work in the elementary differential calculus, the value 
of the curvature of a curve at a point is taken as the absolute value 
of (6) and no attention is paid to the fact that the numerator, 


7a? May be negative as well as positive. However, in surface theory, 


the algebraic value of the curvature is important, and our definition 

of the curvature of the plane section C of the surface S at the point 

P, as embodied in (7), places due regard on the algebraic value of 

TR We append a schematic series of figures here to show, somewhat 

more clearly, what is meant. 


z z z z 
c 
P 
c 
2a 2b 2c 


Fic. 2 
Figure 2a indicates the form of the curve C in the neighborhood 


of P, if, for the value of \ defining C, - >o. Figure 2b indicates the 
— 
same thing if —><o. Finally, if for that value of \ defining C, 


dR? 
d’z dz 
aR? =oand aR” ° at P, the form of C is as shown in one or the other 


of the figures 2c. 


=. 
~ 
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VARIATION OF CURVATURE WITH X 


Equation (11) expresses the algebraic value of the curvature at 
the point P of that plane normal section C of the surface S, made by 
the plane through the z axis whose angle with the (x, z) plane is X. 
Varying the value of \ from o° to 180° continuously corresponds to 
twisting the normal plane continuously around the z axis from the 
(x, z) plane, (A=o°), through its position of coincidence with the 
(y, 2) plane, (A=90°), and finally back to the (x, z) plane for which 
A=0° or 180° (from 180° to 360°, the sections obtained as \ varies 
from o° to 180° are repeated). 

The first fact to notice is that if for all values of \ from 0° to 180°, 
the curvature (11) is always positive, or always negative, the surface 
in the neighborhood of P lies on one side of its tangent plane. If, 
however, as \ varies in this range, the curvature changes sign, part of 
the surface lies on one side of the tangent plane to S at P and part on 
the other. The latter situation is the one shown schematically in 
Figure 1. 

Suppose that we fix our attention to the curvature &, at P of the 
section for which A\=\,; that is, 


(12) ky = rcos? + 2s sin cos Ay + sin? Aj. 


Consider the curvature ke at the same point, of the section whose 
plane is at right angles to that of the first section; that is, that for 
which A=\2=, + 90°. Substituting this value for \ in (11) yields: 


(73) ke = rsin?\ — 2s sind cos\ + # cos? 


The sum of (12) and (13) is: 
(14) 


a constant. To put this interesting result in words: The algebraic sum 
of the curvatures at a non-singular point of an analytic surface of any 
two normal plane sections at right angles to each other is constant. One- 
half of this constant is called the mean curvature of the surface at the 
point. 

The next question to be raised in regard to (11) is whether, as \ 
varies from o° to 180°, the curvature & reaches extreme (maximum 
and minimum) values; and, if so, for what values of \? A maximum 

ak 


dk 
or minimum value of & is obtained whenever — =o and —=o. 
dx 


Differentiating (11), we obtain: 
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dk 


— = (t — r) sin 2\ + 25 cos 2h. 
dx 


(15) 


If this derivative is set equal to zero, and solved for A, we find: 


2s 
(16) tan 2\ = ” 


There are two values of \ lying between o° and 180° which satisfy 
equation (16), and their difference is g0°.! That is to say, if A: is one 
solution of (16) lying in value between o° and 180°, then Ai1+ 90° is 


« 


a 


FIG. 3 


the other, that one being chosen which lies in value in the same range. 
This result may be stated in the form of a theorem: At every non- 
singular point (not umbilical) of an analytic surface there are two normal 
sections of the surface, at right angles to each other, for which the curva- 
tures attain extreme values. 

Combining this theorem with the preceding one, in which we have 
seen that the sum of the curvatures for two normal sections at right 
angles to each other is a constant (cf. (14)), we must conclude that 
at one of these sections the curvature is a maximum and at the other 
it is a minimum. 

* To show that d*k/d\*o0 when dk/d\X=0, we note that, since 
2 cos — r) — 2s tan 2d], 


its value when dk/d\=0 is obtained by substituting the value of \ of equation (16) in 
this last equation. After simplification, this reduces to 


+ — #)? + 45%, 
which can vanish when and only when r= and s=o. A point at which this occurs is 
an umbilical point; and, referring back to equation (11), we see that the curvatures of 
the normal sections to a surface at a point of this type are all equal. Such points will be 
excluded from this discussion. 


1 If r=t and so, the values of \ are 45° and 135°. 


| 
A % 
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If the curvatures & of the plane normal sections be plotted against 
\ for o° SA< 180°, the curve would have the characteristic form shown 
in Figure 3. If the curvature changes sign, as \ varies, the correspond- 
ing curve crosses the \-axis twice. 

It has been shown that the difference in the values of \ for which 
k is a maximum and that for which it is a minimum is 90°. Moreover, 
the algebraic sum of the ordinates (&) for any two values of \ whose 
difference is go° is a constant. 

The directions on the surface at P in which these extreme values 
of the normal curvature occur are called the principal directions. The 
reciprocals, p; and pe, of the extreme values of the curvatures, k, and 
ko, are called the radii of principal curvatures: 

I I 
Pi P2 
The total or Gaussian curvature of the surface at the point is defined 
as the product hike; and one-half the sum: 4 (ki +2) is the mean 
curvature of the surface at the point. 

Many further interesting properties of the curvature of a surface 
can readily be deduced from the analytical method here developed. 
However, it is not necessary for our purpose to delve into these. It is 
well to point out before leaving the matter, that, though, initially, 
we chose the x and y axes as any two mutually perpendicular lines 
in the tangent plane to S at P, we may now, in the light of our re- 
sults, choose them in the two principal directions. Analytically, this 
means that the solutions of (16) are \=o° and A=9g0°; that is, s=o 
in the series expansion (3). 

The equation (16) leads to the results: 


. + 2s + (r — 12) 
(17) sin 2\ = —, cos2\ = ———»> 
w w 
where 
(18) w= V(r — + 


The positive signs in (17) yield the value of \ for which the curvature 
attains one extreme value, and the negative signs yield the other 
value of \, differing from the first by 90°, for which the curvature 
attains the other extreme value. If we set either one of these values 
of \ for \; in (12) and (13), the &; and ke so obtained will be the values 
of these extreme curvatures. 
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Let us note that the difference of equations (12) and (13) may be 
written 


(19) ky — ke = (r — #) cos 24 + 2s sin 2k, 


and, using the results of (17), we find that the extreme curvatures 


satisfy the relation 
r — t)? + 4s? 
E +4 |- 


(20) ki =H 


APPLICATION TO EQUIPOTENTIAL SURFACES 


The torsion balance deals with the equipotential surfaces, which 
are the surfaces perpendicular to the lines of force of the earth’s 
gravitational field. Let the equation of the equipotential surface 
through the center of gravity of the suspended system of a torsion 


balance be 
(21) U(x, y, 2) = constant, 


with respect to the usual coérdinate system used; namely, the z axis 
is the direction of the force of gravity at the center of gravity of the 
suspended system (taken as the origin) directed positively downward, 
the x and y axes are the lines through the origin, tangent to the equi- 
potential surface, bearing north and east respectively. The value of 
gravity at the point is, then, 
(22) aU 
22 g 
02 
The torsion balance observes the values of 


dg og 


2 
(23) Ox dydz 
and 
(24) y 8U 
24 = zy > ’ 


* It will be noticed that the mathematical analysis does not differentiate between 
k, and ke, as to which is the maximum and which the minimum, and cannot de so. This 
becomes obvious when we consider that if, initially, the z axis had been oriented in the 
opposite direction, as it well might have, the maximum curvature would become the 
minimum and vice versa. This corresponds with a reflection of the curve of Figure 3 
in the \-axis; for what was originally z is now —z and the equation (10) changes sign. 
In short, equation (20) states that the numerical value of the difference between the 
principal curvature is w. In the numerical work used in torsion balance calculations, 
however, the curvatures assume, tacitly, numerical algebraic values and the cor- 
responding \ to each of these can be assigned. This will be indicated briefly in the last 
paragraph of the paper. 
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at the point in question. The first pair of equations deals with the 
gradient of gravity and those of the latter pair are referred to as the 
curvature quantities. 

Our problem is to find the values of p, q, 7, s, and ¢ for our surface 
defined in the form (21). This equation defines z implicitly as a func- 
tion of the two independent variables x and y. Recalling that the 
center of gravity of the suspended system bears the same relation 
to the equipotential system through it as does the point P to the 
surface S previously described, we conclude that 


4 Oz 
— = p=o, and — =q=o0 
Ox oy 
Since, however, 
aU oz dU aU 
(25) — —=0, —+——= 
Ox Oz Ox oy oy 
dU aU 


it follows that at the origin, ">. 
x 


Similarly, since at the origin, »=q=o, we find that 
aU 

a. 

Ox? dz Ox? 

aU aU a% 

dy? ds ay? 


— =o 
Oxdy AxdY 


(26) b. 


Hence 


(27) 


and 


* The method of obtaining these derivatives of s with respect to x and y when z 
is defined implicitly, as in (21), is fully discussed in Goursat-Hedrick, Mathematical 
Analysis, Vol. 1 (Ginn & Co., 1904), p. 42. 


By definition 
aU 
s=—;, and —=g. 
Ox? Oxdy dy? Oz 
Ua 
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(28) 


where we have used (22), (23), and (24). The first of these is the differ- 
ence of (26a) and (26b) and the second is equivalent to (26c). 

We must express all our results concerning curvature in terms of 
Us, and U.,, which are the quantities observed by the torsion balance. 
Equation (18) becomes 


(29) w=—V Us? + 4U2,? 
§ 


and equations (17) and (20): 


> + 2U zy + Us, 
(30) sin 2A = =» —— 
V Us? + U,? + 4U.,? 
and 
(31) (ki — ke) = + —V + 4U;z,’, 
g 
respectively. 


If, in (31), we use the positive sign, as is customary, then kz be- 
comes the curvature which is algebraically less in value than &,, and 
the value of \ corresponding to this least curvature is defined without 
ambiguity in sign: 


2U, 
(32) sin 2A = a 
V Us? + 
and 
(33) cos 2A = : 


JV U;? + 4 


This is the value of \ usually used in torsion balance work. Also the 
value of 


(34) R =| g(ki — ke) | = + 
leads to 

Bn 
(35) R= — 


cos 2X 


g 
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ABSTRACT 


The successes which have been obtained with electrical methods of prospecting in 
detecting metallic ore bodies have led to repeated attempts to use them also in oil 
prospecting, both in structural work and in locating the oil itself. These attempts have 
not met with as much success as the application of electrical methods to the location 
of ore bodies; however, the interest in the possibilities of electrical prospecting as ap- 
plied in oil work has been aroused again of late, due to the perfection of the resistivity 
and potential-drop-ratio method in regard to the determination of the depth to geologic 
bodies. 

The object of the writer is to give a summary of the whole field of the resistivity 
methods, with particular reference to the recent developments. The factors affecting 
the resistivity of rocks, and methods for the determination of resistivities of rocks and 
formations are described. Then follows a description of the various surface-potential 
methods. First, the resistivity methods proper are discussed, with reference to electrode 
arrangement, apparatus, and methods of interpretation used. Second, the potential- 
drop-ratio methods are treated, also with regard to the electrode arrangements, ap- 
paratus, and methods of interpretation. The opinions of various authors are discussed 
who have expressed their views about the possibilities of locating oil directly by elec- 
trical methods, and examples are presented of results which have been obtained thus 
far with resistivity methods, both in structural prospecting and in attempting to locate 
the oil directly. 


A. PRINCIPLES AND History OF RESISTIVITY PROSPECTING 


In prospecting for oil structure, the four major geophysical 
methods—gravitational, magnetic, seismic, and electrical—have been 
widely used. Both the seismic and electrical methods have a distinct 
advantage over the gravitational and magnetic methods: the possi- 
bility of controlling the depth of penetration. This is an important 
factor in the interpretation of the results, as not only the physical 
characteristics but also the depths of geologic formations are ob- 
tained. 

Thus, the seismic and electrical methods as applied in oil prospect- 
ing have a number of features in common; on the other hand, they 


1 Read before the Association of the Oklahoma City meeting, March 25, 1932. 
(Series of publications No. 43, Dept. of Geophysics, Colorado School of Mines.) 


2 Professor of geophysics, Colorado School of Mines. 
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also differ in a number of respects. The seismic method is essentially 
a dynamic method (velocity measurements) while the electrical 
method is essentially a static method (determination of potentials 
and electromagnetic fields), disregarding for the moment that the 
element of time enters into the latter through frequency and phase. 
The seismic and electrical methods each have a dual advantage; 
namely, that two types of the energy may be observed; in seismic 
methods, the refracted and reflected energy; in electrical methods, 
the potential distribution and the electromagnetic field of the ground 
currents. As far as seismic methods are concerned, the variety of 
application is very nearly exhausted with the foregoing. In electrical 
methods, added advantage is obtained by varying the type of appli- 
cation of primary energy: first, in regard to frequency; second, by 
employing either galvanic or inductive coupling. 

Notwithstanding such advantages, electrical methods, as a whole, 
have not found such wide application in oil work as seismic prospect- 
ing. One reason for this is their more complicated field technique when 
the horizontal radius of operation is increased. Another reason is the 
limit in depth penetration when the frequency is increased. A third 
reason is the following: in most geophysical methods, the sought for- 
mations must have a thickness comparable with their depth. This 
holds for the seismic refraction method and the electrical-potential- 
methods. Noted exceptions are the seismic reflection-method, in 
which even a comparatively thin formation may still be noticed, 
provided it has good reflecting characteristics, and the inductive 
electrical method, in which a thin layer may also be determined, pro- 
vided it excels in conductivity. However, there again, the electrical 
methods are at a disadvantage, as the depth penetration, according 
to practical experience, of the inductive method is not nearly as great 
as that of the seismic reflection method. A fourth reason, finally, for 
the more limited application of electrical methods is purely geological ; 
namely, that the factors which can modify the electrical characteris- 
tics of formations in both horizontal and vertical direction (such as 
pore volume and electrolytic conductivity of moisture) are more 
numerous than similar factors affecting their elastic characteristics. 
Theoretically, all these obstacles could be amply compensated if it 
should be possible to perfect electrical methods to such a degree as to 
enable the direct location of petroleum. However, the reasons that are 
a great handicap to such effort are again purely geological in nature; 
namely, the variable factors influencing the electrical characteristics 
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of formations as already stated are those which in most cases mask 
the effect of the oil. 

On the other hand, it must be admitted that notable advances 
have been made in approaching the solution of the problem, at least 
for favorable conditions. Furthermore, remarkable progress has been 
made in reducing the difficulties which have been outlined in the pre- 
ceding paragraph; that is, attempts have been made to increase the 
depth penetration, and to perfect the technique in such a way as to 
obtain a more interpretable type of indication. 

It is quite logical that the resistivity and surface-potential methods 
offered the greatest possibilities in this respect; for the first necessity 
in attempting an increase in depth penetration is a decrease in work- 
ing frequency. For inductive methods, this would materially reduce 
the current strength induced in subsurface conductors; for electro- 
magnetic methods, it would embody a complication in methods of 
observation (replacement of the telephone by the galvanometric- 
rectifier method, or, if the telephone is to be retained, the use of a 
frequency converter). 

From what has gone before it is obvious that an increase in depth 
penetration should be concomitant, if possible, with a greater pre- 
cision in determining effects from such depths. Therefore, efforts have 
been made to find such electrode arrangements which would give the 
most distinct indication of formation boundaries, as the original 
Wenner 4-terminal, and similar electrode assemblies do not always 
give readily interpretable curves. Probably the most successful ac- 
complishment along this line has been the perfection of the potential- 
- drop-ratio method which required an appreciable modification in the 
technique of the resistivity measurement which had been employed 
up to that time. The potential-drop-ratio method permits the de- 
termination of small differences in apparent resistivity with a great 
deal of precision, requires no connection of the ratiometer to the 
power electrodes, and gives well interpretable curves. 

There is a great deal of literature on the resistivity method (see 
reference list at the end of this article), chiefly referring to the appli- 
cation of the method in mining and civil engineering. In the geophysi- 
cal laboratory of the Colorado School of Mines, extensive studies have 
been made on the application of resistivity methods, not only to the 
problems already mentioned, but also to structural oil prospecting. 
The Gish-Rooney equipment was perfected, model experiments were 
made, a Racom outfit was acquired through the courtesy of the 
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Swedish American Prospecting Corporation, and a number of en- 
couraging results were obtained on geologic structure. 

The object of the writer is to give a preliminary account of some 
of the results obtained, and to present a review of the entire field, 
covering both technique and results obtained by other authors. 

As our activities are not limited to research in electrical prospect- 
ing alone, the work has been carried through a comparatively great 
length of time, and various associates of the writer have been en- 
gaged in the work. Grateful acknowledgment is made for the con- 
tributions of data to T. A. Manhart, C. D. Keen, J. A. Malkovsky, 
D. H. Griswold, to the Swedish American Prospecting Corporation 
for part of the equipment, and to Harry Aurand and R. Clare Coffin, 
of the Midwest Refining Company, for permission to use some of the 
data and experience which this company has obtained in resistivity 
work. 

As the name indicates, resistivity methods are electrical methods 
in which the variation of resistivity is measured on the surface. Nu- 
merous methods are available for measuring resistivities, either in the 
well known resistivity bridges for D. C. or A. C., or by the use of 
ohmmeters, or finally, by separate measurements of voltage drop and 
current. As for known or constant current, the resistivity follows 
immediately from the potential distribution, all methods for the 
determination of potential differences or potential ratios on the sur- 
face, the so-called surface-potential methods, may also be included 
under the heading of resistivity methods. 

It is not surprising that attempts were made rather early to use 
the simple idea of resistivity measurements in prospecting. As early 
as 1900, Brown and McClatchey applied for a patent on resistivity 
measurement in this country, and at the same time Daft and Wil- 
liams, in their English and American patents, suggested the use of 
potential-difference observations for resistivity work. W. Petersson 
(ref. list No. III,) describes a qualitative resistivity method for the 
location of ore, consisting of a buzzer and telephone, the intensity of 
the sound received being an indication of ground conductivity. The 
method was used successfully in Sweden in 1906. 

The early attempts at the resistivity method were, however, only of 
a qualitative nature, as neither potentia! differences for fixed electrode 
spacings, nor even equipotential lines were observed. 

The credit for making the first systematic studies of ground re- 
sistivities goes to C. Schlumberger, who began work on the method in 
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about 1913. The first year in which he applied the method of resis- 
tivity mapping was 1920, and numerous structural studies for com- 
mercial purposes were carried out by this company in subsequent 
years. 

Schlumberger obtained a French patent on his resistivity method 
on September 15, 1925. 

Considerable interest was attracted at the same time to an identi- 
cal method, described by Gish and Rooney, December, 1925 (ref. 
list No. IIIz), which was used in determining earth resistivities for 
the study of earth-current and related magnetic phenomena. Their 
method was based on Wenner’s method of measuring earth resistivity, 
published in rors (ref. list No. II;). 

The Gish-Rooney method was taken up with considerable interest 
by a number of mining companies; primarily, the Michigan School 
of Mines has been largely responsible for the advances of this method 
in mining. The papers written by W. O. Hotchkiss, J. Fisher, and 
W. J. Rooney (ref. list No. IIIs) in 1929, report on the results ob- 
tained in the Michigan copper and iron country. 

At about the same time, the U. S. Bureau of Mines became inter- 
ested in resistivity work, and codperated with A. S. Eve and D. A. 
Keys in the experimentation with electrical methods. Their studies 
have continued up to the present time, and the results are recorded 
in numerous papers, written by F. W. Lee, A. S. Eve, D. A. Keys, and 
J. H. Swartz (ref. list Nos. II, Ie, ITT3,4,9,19,29,30,33,37)« 

A number of fundamental papers also appeared on the theoretical 
foundation of the method, written by J. N. Hummel, W. Weaver, 
D. O. Ehrenburg, Lancaster-Jones, Tagg, R. J. Watson, T. Roman, 
L. J. Peters, J. Bardeen, S. Stefanesco, and C. and M. Schlumberger 
(ref. list, section IT). 

Meanwhile, in 1928, an important application of the resistivity 
method had been found, the determination of depth to bedrock in 
dam sites. Papers by I. B. Crosby and E. G. Leonardon (ref. list 

12,18,20,25, 32,35) deal with this application. 

As far as the use of resistivity methods in oil work is concerned, the 
Schlumberger Company had proved the practicability of the method 
in structural studies since 1921. Work was then begun in the Pechel- 
bronn oil region and continued until 1926. In Roumania, commercial 
work was carried out, among others, for the Steaua Romana in 
1923-1926. Salt domes were located in the Alsace region in 1926-1927. 
Resistivity prospecting for oil structures was begun in the United 
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States by 1925, when the Schlumberger Company was engaged to 
work for the Roxana Petroleum Corporation and the Shell Company 
of California. The work for the latter continued to about 1929. 

A number of oil companies have carried on extensive studies with 
their own electrical resistivity and potential equipment, such as the 
Sun Oil Company, the Pure Oil Company, and the Midwest Refining 
Company. 

In addition to the Schlumberger Company, several other geo- 
physical prospecting companies have since taken up studies, such as 
the McCollum Exploration Company, the Swedish-American Pros- 
pecting Corporation, the Radiore Company, the Elbof Company, the 
Geophysical Service Inc., and the International Geophysics, Inc. 
Very little has been published on the successes of structural resistivity 
prospecting; practically everything is by members of the Schlumberger 
Company (ref. list No. ITTz,10,12,16,18,20). 

One of the outstanding developments in resistivity-prospecting 
methods in late years was the perfection of the potential-drop-ratio 
method. This was suggested by three authors at almost the same 
time. Although a 3-contact ratic arm bridge was devised by A. B. 
Edge as early as 1925 and was used in northern Rhodesia during 
that year, and although an A. C. potential ratiometer was used, up 
to the end of 1929, by the I.G.E.S. in Australia, it was not until 
January, 1931, that the method was first described by Edge (ref. list 
No. IV,). Already before this publication appeared, Koenigsberger 
(ref. list No. IV:) had suggested a potential-ratio method in 1930. 
Two months after Edge’s article had been published, H. Lundberg 
and Th. Zuschlag published their description of the Racom (ref. list 
No. IV;). 

The potential-drop-ratio method, on account of its sensitivity 
and accuracy, gives very good results when suitable instruments are 
used, and it is believed that it is only in its initial stage of develop- 
ment. 

There follows now a discussion of (1) the resistivity of formations 
and methods for its determination; (2) a description of the technique 
of resistivity and potential-drop-ratio methods; and (3) a discussion 
of the results obtained in structural work. 


B. RESISTIVITY OF FORMATIONS AND ROCKS 


With the exception of the few rocks which contain some sort ofa 
metallic mineralization (mostly sulphidic), and the conductivity of 
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which is, therefore, of a metallic nature, the conductivity of all the 
other rocks and formations is of a purely electrolytic character. There- 
fore, all these rocks are highly resistant when dry; or, in other words, 
the conductivity of a rock is proportional to the amount of water 
which it contains, and to the quantities of ionized salts dissolved in 
the water. 


I. FACTORS AFFECTING ROCK RESISTIVITY 


From what has been previously stated, it is obvious that the re- 
sistivity of the formations which are of importance in electrical pros- 
pecting is affected primarily by the following factors: (1) the resis- 
tivity of the mineral constituents; (2) the resistivity of the water 
filling the pores; (3) the pore volume; (4) the shape of the pores; (5) 
the temperature; and (6) the pressure at the particular depth of the 
formation under investigation. 

As far as the resistivity of the mineral constituents is concerned, 
which make up most of the rocks encountered in oil resistivity pros- 
pecting, this factor is of comparatively small influence. As has been 
stated before, only rocks with a metallic mineralization have a com- 
paratively low resistivity, while the resistivity of the mineral matter 
in most barren sedimentary and igneous rocks is so high (of the order 
of 10°-10 ohms cm.~*), that they may be considered as insulators 
for the purposes of geophysical prospecting; that is to say, in most 
sedimentary rocks, the influence of the pore volume and of the con- 
ductivity of the waters filling the pores is so much more predominant 
a factor in determining the conductivity of the rock, that the influence 
of the resistivity of the mineral matter may be disregarded in com- 
parison. 

The resistivity of the water filling the pores in such rocks as are 
of importance in resistivity oil prospecting is of paramount importance. 

Following Sundberg (ref. list No. Ip,10), the impregnating waters 
may be divided into the following groups. 


A. Surface waters (waters from the surface of the ground down to and in- 
cluding the ground water). 


(a) Fresh water, such as rain, snow, lake, river, and ground water, and 
the normal moisture in the upper soil. The percentage of dissolved 
substances ranges from o.o1 to 0.1 at an average. These dissolved sub- 
stances are ordinarily carbonates and silicates and small amounts of 
chlorides and nitrates of calcium, magnesium, sodium, and potas- 
sium. The resistivity of such surface waters varies at an average be- 
tween 30,000 and several hundred thousand ohms cm~*. 
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(b) Mineral waters with varying concentrations. Such waters are ordi- 
narily locally limited, and are not commonly encountered on the sur- 
face of oil-bearing structures; they are not of as great importance 
for our problem as those previously mentioned, and the deep wat- 
ers to be discussed now. 


B. Deep waters. 
(a) Connate waters. 

These waters are ordinarily encountered in oil fields. As the name indi- 
cates, their origin dates back to the time when the formations were 
laid down in sedimentary basins. They contain chlorides of sodium, 
potassium, magnesium, and calcium. They differ from sea water, 
however, by the abundance of calcium chloride and by the absence 
of sulphates. They are distinguished by a very wide range of concen- 
trations, from dilute to concentrated solutions. Consequently, their 
resistance also varies within wide limits, to as low as a few ohms cm-*, 
Mine waters. 

These waters are ordinarily solutions of metal sulphates, but they 
also contain carbonates of sodium, calcium, and magnesium. They 
vary widely in resistivity which may be as low as 30 ohms cm~. 


The most important chemical constituent of waters of all kinds 
which determines their conductivity is the chlorine content. Figure 
1 b illustrates the relation between resistivity and chlorine content for 
three ranges of concentration. 


b 


a3 
Fic. 1.—Factors controlling resistivity of rock waters (after Sundberg). 


a: Relation between resistance factor P=p,/p and percentual volume » of water in 
rocks. (p,=resistivity considering grain arrangement; p=resistivity not considering 
grain arrangement.) 


b: Relation between chlorine content of waters and their specific electric resistance. 
c: Relation between temperature and specific electric resistance for electrolytes. 


The next outstanding factor determining the resistivity of rocks 
and formations is the pore volume. The pore volume depends on the 
size and shape of the grains, and their mutual arrangement. The 


| 
| c 
| 
| | 
j 
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porosities of such formations as are ordinarily encountered in oil 
prospecting cover a wide range, approximately between 10 and 45 
per cent. 

The pore volume, of course, is equal to the maximum amount of 
water which is possible in a rock. This condition is generally assumed 
to exist below the ground-water level, except where the water has 
been replaced by gas or oil. Above the ground-water level, there are 
considerable variations, which account for the varying values of re- 
sistivities that are often encountered near the surface in resistivity 
work. The moisture content of the soil is generally very low near the 
surface on account of evaporation, then increases rapidly with depth, 
then decreases and reaches a maximum in the ground-water level.' At 
the immediate surface, the moisture, and therefore the resistivity, 
varies considerably with the meterological factors. During rainfalls, 
the resistivity will usually be very great on account of the purity of 
waters; during droughts, highly conductive waters migrate upward, 
resulting in a decrease of the resistivity of the formations concerned. 

Due to the fact that the conductive medium, the water, in a rock 
is not everywhere of the same section, but varies from a very thin 
section where the grains are touching each other to a greater section 
where grains are not in contact, the relation between pore volume and 
rock resistivity is not altogether simple; it may, however, be worked 
out for certain simplified assumptions; for instance, for the case that 
the grains are spheres of equal size or for the case that the pores have 
the shape of three mutually perpendicular and parallel systems of 
tubes. Figure 1 a shows the relation between a resistance factor P, 
which is the ratio of the actual water resistivity depending on porosity 
to the water resistivity, without regard to porosity (Chapter B, II, d), 
and the water volume of a rock in percentage. 

Finally, the temperature and the pressure influence the electro- 
lytic resistivity of rocks. Figure 1 c shows the influence of tempera- 
ture; an increase in temperature decreases the resistivity consider- 
ably; the relation shown in the figure holds for any kind of an electro- 
lyte, and indicates that a rise in temperature by only 25° C. from zero 
is sufficient to increase the conductivity of a rock approximately 100 
per cent. 

The effect of pressure is much smaller and is illustrated in Figure 
2. The effect is different on concentrated solutions from what it is on 
solutions of smaller concentrations. At any rate, with the maximum 


1 For details, see ref. list III, 38. 
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effect on a concentrated solution, a pressure corresponding with a 
depth of burial of nearly ro kilometers would bring about a change 
in resistivity of only ro per cent. 

More details on the effect of water on rock resistivities may be 
obtained from K. Sundberg’s papers (ref. list I _10). 


II. DETERMINATION OF RESISTIVITIES 


The experimental determination of the resistivity of rocks or for- 
mations is by no means a very easy matter. This is chiefly due to the 
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Fic. 2.—Relation between pressure and resistivity of electrolytes (after L. H. 
Adams and R. E. Hall). 


fact that a sample can not well be kept in the same condition after 
being removed from its location as it has been in before it was dis- 
turbed, chiefly in reference to its mechanical constitution and its 
moisture content. Therefore, experiments to test the resistivity of 
rocks are usually made, whenever possible, on outcrops, or by means 
of special equipment, in wells (Schlumberger method). If there is no 
possibility of avoiding the testing of samples, the tests should be 
made in a portable apparatus immediately after the sample has been 
taken, or else, the resistivity of a formation may finally be computed 
from the pore volume of the constituent rock and the conductivity 
of the impregnating water (Sundberg method). 


C. A. HEILAND 


a. DETERMINATION OF RESISTIVITIES ON SAMPLES 


A simple arrangement for the determination of resistivities of 
solid rock specimens is shown on the right side of the panel illustrated 
in Figure 3. 

The arrangement consists of a frame supported on springs; one 
round electrode is connected to the bottom of this frame, the other 
is attached to the bottom of a spindle. The sample to be tested is 
placed between the electrode on the spindle and that on the frame, 


Fic. 3.—A.C. and D.C. resistivity bridge (C S M). 


with some tin foil between sample and electrodes. Mercury electrodes 
may be used on drill core samples instead of the solid metal plates. 
This arrangement may be placed in a Wheatstone bridge. 

The samples should, preferably, have a simple geometric shape. 
If / is the length, s the section of the sample, R the observed resistance 
and p the resistivity to be determined, p= S 

The resistivity is usually expressed in ohm-cm (see ref. list I, 12). 
In some publications ohm- meters are used; . :sistivities expressed in 
the latter units are one hundredth as large as resistivities expressed 
in the former units. 
The bridge shown in Figure 3 has on the left side the controls for 
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high frequency tests, and on the right are the meters and switches for 
D.C. testing. The bridge is A.C. operated, and consists of an oscil- 
lator, the frequency of which may be varied within wide limits in the 
high frequency range. A pickup coil with variable coupling transfers 
the energy of the oscillator to the testing circuit, which is essentially 
an A.C. Wheatstone bridge. Non-inductive resistors are used; the 
switches for the 10, 100, 1,000, 10,000 and 100,000 ohm ranges are 
seen on the panel. The same resistors are used when the bridge is 
operated on D.C. The binding posts on the extreme right edge of the 
panel are for making connections to the testing battery or a low fre- 
quency oscillator or commutator.! 

The bridge is only for small samples of regular shape, such as 
diamond drill cores and the like. For larger specimens of consolidated 
rocks, we use a different bridge which consists of a frame with two 
spindles to hold the specimen, and four arms spaced at equal intervals 
which carry contact points at their end. The Gish-Rooney (or Wen- 
ner) method is applied in testing the specimen, and the Gish-Rooney 
equipment as used in the field can be connected directly to this bridge. 

M. W. Puilen (ref. list No. Ig) uses also direct and alternating cur- 
rent in his resistivity tester. The resistivities are measured by com- 
paring the resistance of the specimen with a known resistance, that 
is, by switching the current from one circuit containing this speci- 
men over to the other circuit containing the known resistance. The 
electrodes used were chiefly mercury pools, held in place upon the 
specimen by wax dams if the specimen had an irregular shape, or by 
a bakelite cup and a paper collar for drill cores. To eliminate surface 
leakage of current across the specimen, guard rings were applied. In 
using continuous current, the observed resistivity values were found 
to change with time on account of the variable polarization. There- 
fore, discontinuous current was found to be more satisfactory, which 
was obtained by shortening the regular 60-cycle A.C. from a light 
circuit through a potentiometer (for obtaining different test voltages) 
and using a copper oxide rectifier to make the use of a direct current 
galvanometer possible (Fig. 4 a). 

For large specimens of irregular shape, not only the Gish-Rooney 
four-electrode method may be used, but also two like electrodes of 
regular geometric shape may be employed, provided their distance 
is small as compared with the dimensions of the sample. The theory 


1 This resistivity bridge was developed largely after the bridge used by the Radiore 
Company. 
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of such an arrangement (Sundberg, ref. list No. Is) shows that for 
square electrodes with sides of 0.37 cm. the measured resistance is 
equal to the specific resistance of the specimen. 

The bridges described before may only be used for solid speci- 
mens. They are not applicable for soils. For testing of soil samples in 


Fic. 4a.—Circuit of A.C. rock-resistivity tester, operated from light socket (after 
M. W. Pullen). 


Telephone 


Battery 


Batlery Switch nterrupter 
Fic. 4b.—Leeds and Northrup soil-resistivity tester (after R. O. E. Davis). 


the field, the Leeds and Northrup soil bridge may be applied. The in- 
strument (known as the Wheatstone-Kirchhoff bridge) as shown in 
Figure 4 5, is a modified form of a Wheatstone bridge. As a source of 
energy, a buzzer operated by a dry cell is used; consequently, a tele- 
phone may be employed to indicate the balance. The soil to be tested 
is placed in a cup; this represents the variable resistance in one arm 
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of the bridge; the fixed comparison coils are in the other arm of the 
bridge; both connect to the variable slide wire and the secondary of 
the buzzer coil. The comparison coils are made up of three fixed re- 
sistances of 10, 100, and 1,000 ohms. An additional resistance of 100 
ohms may be used in series with the cup when its resistance is low. 
The balance is adjusted on the slide wire, with proper setting of the 
fixed resistances, and the reading on the slide wire is multiplied by 
the resistance of the comparison coil used. 

The use of the cup is very convenient as not only soils, but also 
electrolytic solutions may be tested in it. It may be standardized by 
using solutions of known concentrations. The results must be re- 


Fic. 5.—Shepard earth-resistivity meter. 


ferred to a normal temperature. The publication by R. O. E. Davis 
(ref. list No. I,) gives a number of tables on resistivities of soils and 
solutions which will be found valuable when using this bridge. 


b. DETERMINATION OF RESISTIVITIES ON OUTCROPS 


On account of the difficulties which are inherent in the accurate 
determination of the resistivity of samples and which have been dis- 
cussed before, determinations of the resistivities on outcrops, or, 
generally speaking, on the ground surface are ordinarily much more 
satisfactory. 

Quite a few methods may be used for this purpose. For approxi- 
mate tests, the Shepard earth-resistivity meter may be employed 
(Fig. 5). It is an instrument working on D.C.; the ohmmeter and the 
battery are mounted on one probe; the probes are made of iron, and 
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the cathode is larger than the anode, reducing polarization to a great 
extent. For this purpose, the electrodes have a large ratio of areas, as 
seen in the figure.' 

For more accurate work on the surface, the Megger and McCol- 
lum’s earth-resistivity testers, or the 4-terminal Gish-Rooney equip- 
ment is used, which will be described later. 

J. Koenigsberger (ref. list No. I;) uses two circular iron electrodes 
about 25 cm. in diameter, and about 5 mm. in thickness, provided 
with a contact substance which is either hematite powder or a clayey 
paste with an aqueous solution of FeCl;, or FeSO;, or of 10 per cent 
NaCl. Where the necessity arises of using larger electrodes, iron 
screens may be employed. As a source of energy, a buzzer is used at a 
frequency from 100-400 cycles; the electrodes are placed in a Wheat- 
stone bridge with a telephone as indicator. Koenigsberger gives the 
formulae to be used with this arrangement of the electrodes. 


¢. DETERMINATION OF RESISTIVITY IN WELLS 


This is the only method for obtaining accurate resistivity values 
for any formation below the surface. The equipment and technique 
has been perfected by the Schlumberger Company (ref. list No. 
Is,1:), and excellent results have been obtained with it. The arrange- 
ment is shown in Figure 6. It consists of four electrodes, two of which 
are supplied with current, while the potential difference between the 
two intermediate electrodes is observed; one current electrode is 
grounded at the surface. The arrangement of the three other elec- 
trodes always remains fixed at any depth. With the symbols used in 
the figure, the resistivity at any depth is given by the relation 


if i is the current read on a milliammeter in the current circuit and 
AV is the potential difference read on a potentiometer between the 
electrodes 2 and 3. Excellent results have been obtained by this 
method, and have occasionally resulted in the detection of oil and 
coal horizons which had been overlooked in drilling. The results, 
together with some outstanding examples, will be discussed later. 


1 See: S. Ewing: “Soil Survey Methods,” Oil and Gas Journal, April 21, 1932, p. 42. 
C. R. Weidner and L. E. Davis: “Relation of Pipe Line Currents and Soil Resis- 
tivity to Corrosion,” The Oil Weekly, December 4, 1931, p. 26. 
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d. DETERMINATION OF RESISTIVITIES FROM WATER ANALYSES ‘ 
AND POROSITY DETERMINATIONS 

If it is not possible to determine underground resistivities by the 
Schlumberger method, that is, if there is casing in the well at the 
depth to be investigated, or if the wells are dry, sufficiently accurate 
estimates of resistivities of formations may often be obtained by a 
method suggested by Lundberg (reference list No. I9,10), provided 
water analyses are available for the formations under investigation, 
which is often the case. If the water analyses can, in addition, be 
supplemented by a determination of the pore volume of the forma- 


Fic. 6.—Schlumberger’s method of measuring resistivities in wells (after C. and 
M. Schlumberger). 


tion from an investigation of drill cores, for instance, a fairly accurate 
determination of the resistivity of the formation can be made. 

If the porosity of a sample is known, the resistance factor P may 
be found from the diagram in Figure 1 4, if it is assumed that all the 
pores are filled with water; this is ordinarily true, as already stated, 
below the ground-water level. If the chemical composition of the 
waters at a certain depth is known, and if it is assumed, with suf- 
ficient accuracy, that only sodium chloride is in solution, the diagrams 
in Figure 1 6 may be used to determine the resistivity of the water; 
to convert the results into actual resistivities at a certain depth, the 
temperature correction illustrated in Figure 1 ¢ may be applied. 
The resistivity factor P is of the order of 50 to 100 for dense lime- 
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stones and sandstone, 20 to 40 for clays and sands, and 2 to 20 for 
porous sands, clays, and soils. 


C. Resistiviry METHODS 


In the two sections that follow, the principles, apparatus, and 
electrode arrangements used in the resistivity and potential-drop- 
ratio measurements are discussed first. After that, some outstanding 
results are described which have been obtained with both methods 
on geologic structure of various types, as well as in direct prospecting 
for oil. 

Though the purposes for which the resistivity and potential-drop- 
ratio methods are used are usually identical, these two methods differ 
essentially in electrode arrangement and technique of measure- 
ment. 

In the resistivity methods proper, current is usually supplied to 
two points, and the potential distribution between these two points 
is investigated by either one or two electrodes. The current is deter- 
mined which flows between the two outside current terminals, and 
the potential difference is observed which exists either between the 
two potential electrodes, or between the one potential electrode and 
one neighboring current electrode. The ratio of voltage and current, 
multiplied by a constant depending on the mutual arrangement of 
the electrodes, gives the resistivity of the surface formation, if only 
one formation exists within the reach of the whole contacting arrange- 
ment. If more formations are present the “apparent” resistivity is 
measured as a function of the electrode spacing, and thus, as a func- 
tion of the depth to the interfaces of the formations involved. 

In the potential-drop-ratio methods, current is again supplied to 
two points in the area under investigation; but this time, the in- 
vestigation of the potential distribution takes place outside of the 
pole doublet and is usually made in the vicinity of one electrode. The 
drop of the potential in the vicinity of this one electrode is intimately 
related to the existence of resistivity changes with depth. Experience 
has shown that the most interpretable results are obtained by not 
investigating the drop of the potential itself, but by determining the 
ratio of potential drops in two adjacent intervals. That is to say, the 
testing device in this method requires only 3 electrodes, and is in- 
dependent of any connection to the current circuit or to the source of 
energy, which gives these methods a distinct advantage over the 
resistivity methods proper. 
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We shall first discuss briefly the electrode arrangements and ap- 
paratus used in the resistivity methods proper. 


vy 6 ha 


ae 
Fic. 7.—Electrode arrangements in resistivity prospecting. 


I. ELECTRODE ARRANGEMENTS 

a. POSSIBLE ARRANGEMENTS 
Five possible electrode arrangements which can be used in re- 
sistivity work are described in detail, and formulae are given for each 
arrangement, in J. N. Hummel’s publication (reference list No. Ii) 
and Figure 7. They are summarized briefly as follows. 
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(rt) Wenner’s arrangement, also called simply the 4-terminal method. (This 
is the electrode arrangement most widely used.) The potential electrodes 
are arranged in line with the current electrodes, and the spacing of all 4 
electrodes is equal. That is, the distance between the outside terminals is 
3a if a is the electrode interval. 

Single-current-probe method: One current electrode is far apart from the 
other, and the two potential electrodes are used in the vicinity of one 
electrode. If their interval is a, they are so used in the vicinity of the cur- 
rent electrode that the distance of the current electrode from the potential 
electrode next to it is also always kept equal to a. The potential electrodes 
need not be in a straight line with the current electrodes. 

This third possible electrode arrangement is practically identical with 
that described in No. 2, except that the second current electrode is no 
longer in infinity, and that the potential electrodes are now in a straight 
line with the current electrodes. 

The fourth method may be called the single-potential-probe method. 
Again two current electrodes are used, but only one potential electrode, 
the position of which is varied between the two current electrodes. 

(5) The fifth method is identical with the fourth, but this time one current 
electrode is at infinity. This method may thus be called the single-current 
and the potential-electrode method. 


(2 


(3 


(4 


Two more methods are possible and have been used, both being 
modifications of the second and the first method already described. 

The first of these is the single-probe method by Eve and Keys 
and is similar to the second method (Fig. 7, ZZ a), with the exception 
that the interval between the two potential and one current electrode 
is no longer equal; only the spacing between the potential electrodes 
remains constant. 

The second of these additional methods is a modification of Wen- 
ner’s and Gish-Rooney’s 4-terminal method. An additional electrode 
is provided between the potential electrodes; measurements take 
place as before, that is on four electrodes. Two sets of results, how- 
ever, are now obtained: one (at the right), using the electrodes R and 
S; the second (at the left), using the electrodes P and S. This is called 
Lee’s method (see ref. list No. III, 37) of partitioning.' 


1 A third additional method may be mentioned here, which is the “Potential center 
displacement” method of S. H. Williston and C. R. Nichols, described by F. H. 
Lahee in the third (1931) edition of his Field Geology (p. 694). The electrode arrange- 
ment is that shown in IV in Figure 7. The object of the method is the measurement of 
the displacement of the potential center from its geographic center between the two 
primary electrodes which it would occupy in homogeneous ground. In surveying an 
area, the whole contacting arrangement is stepped forward as in resistivity mapping, 
the electrode base remaining constant. The depth of penetration is roughly equal to 
half the length of this base; an area may be covered with different exploration depths 
as in resistivity mapping. The potential center displacements are at a maximum over the 
lateral edges of resistant bodies; zero displacements between two opposite maxima 
would thus indicate the crest of an anticline, etc. 
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In Figure 7, these electrode arrangements are shown. Also the 
formulae are given for the computation of the resistivities of the 
ground if no unhomogeneity is within the reach of the measuring 
arrangement. Otherwise, the resistivity values obtained by the ap- 
plication of the formulae represent apparent resistivities only, from 
which the true resistivities, as well as the depth to the formation 
boundaries, may be obtained by applying methods to be described 
later. 

b. CUSTOMARY ELECTRODE ARRANGEMENTS 


Of the many electrode arrangements previously described, only 
a few are customarily applied and it all depends on the purpose for 
which a resistivity investigation is made, which method is used. 

There are two distinct applications of the resistivity method. 
First is the method of resistivity mapping. The object of this method 
is to obtain a contour map of the area, showing lines of equal re- 
sistivity. These lines represent the ground resistivity only to a definite 
depth, and, to facilitate the interpretation of the results, two or more 
contour maps may be made for the same area, representing the re- 
sistivity down to two depths of investigation. This can be accom- 
plished simply by working with two fixed electrode separations, as 
described hereafter. The second method is the so-called method of 
electrical drilling. The results are to give the vertical variation of the 
resistivity at one point only and the depth at which any changes 
occur. The only way to accomplish this is, of course, to survey along 
a line beginning at the point where the vertical resistivity variation 
is to be tested; and this line is the longer, the greater the desired 
depth penetration. It is thus seen that accurate results in this method 
depend altogether on the question whether or not the vertical differ- 
entiation of the resistivity carries unchanged in a horizontal direction 
below the last point of measurement. 


b 1, RESISTIVITY MAPPING 

In this method, usually the 4-terminal electrode arrangement is 
employed and the whole contacting outfit is carried over the area 
to be surveyed with a constant electrode separation.’ If it is desired 
to obtain two equi-resistivity maps covering two depths of penetra- 
tion, two electrode separations are used at each point of observation. 
The results of some equi-resistivity surveys are shown in Figures 26, 
28, 29, 30. 


1 As a matter of interest it may be mentioned that the equiresistivity method has 
recently been used for oil prospecting in shallow water (ref. list No. III, 41). 
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b 2. VERTICAL ELECTRICAL DRILLING 

For this purpose, several methods of those already enumerated 
may be used. 

(1) The method that has been most commonly used is probably 
the 4-terminal Gish-Rooney-Wenner method (Fig. 7 J). The elec- 
trodes are arranged always symmetrically around the point to be 
drilled electrically; that is to say, this point is located where the 
electrode S in the Lee partitioning method would be placed. By ap- 
plying this modification of the Gish-Rooney method, one readily 
obtains the distribution of resistivity at the north and at the south 
(or any other two opposite directions) from the midpoint; that is to 
say, it is possible to determine by this method if there is a dip of the 
formations, or if formations come in north of this point that may be 
absent in the south. 

(2) The method that has been used perhaps as frequently as the 
4-terminal method for electrical vertical drilling is the single-current- 
probe method, that is, the investigation of the potential distribution 
in the vicinity of one current electrode only. For this purpose, either 
the electrode arrangement shown in JJ in Figure 7 may be used (ER 
equal to 2a) or the modification shown in JJa. As indicated in Figure 
7 II, it is not necessary that the two potential electrodes be in line 
with the current electrodes; in fact, to eliminate completely the 
influence of the second current electrode, which is assumed to be at 
infinity, the line of measuring the potential should be located on a 
circle through with the distance EZ, EZ, as radius around £2; or, if 
that is not feasible, a close enough approximation may be obtained 
by running the potential measurements on a line at right angles to 
E,E:. This is the method employed by Ehrenburg and Watson (with 
EP =a, ER=2a) (reference list No. II). Gilchrist has applied (refer- 
ence list No. IIIs4) a modification of this last method by using two elec- 
trodes E, at either side of E;, and passing currents of equal strength 
through both circuits (Fig. 7 IZ 6). Instead of using only one current 
base E2F, at right angles to the potential line, any number of current 
bases, arranged in the form of a symmetrical star, may be employed. 


II. APPARATUS FOR RESISTIVITY WORK 


We now come to a description of the equipment used in resistivity 
work. In order not to make this paper unduly long, rather extensive 
use of figures will be made instead of giving lengthy descriptions of 
the apparatus. Further details may be obtained from the literature 
given in the reference list. 
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a. GENERAL: POWER SOURCES 


As a source of power, one Radio B battery, or two, may be suffi- 
cient, depending on the length of the current basis, of course. If very 
great depth penetration is required, a D.C. generator is used, driven 
by a gasoline engine. When using a Megger, a hand-driven D.C. 
generator is used, which is incorporated in the meter box. Due to the 
fact that polarization must be eliminated, if non-polarized electrodes 
are not used, a commutator is provided in the measuring circuit, 
which brings it about that in reality the current used is not strictly 


Potentiometer 
Miliammerer 


Electrodes making contact with Earth 
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EFarth 


Fic. 8.—Schematic diagram of Gish-Rooney circuit. 


D.C., but an alternating current of low frequency. Induction effects 
are, therefore, sometimes observed between the current and the 
potential leads; besides, not the ohmic resistance of the ground is 
obtained, but the impedance for the particular frequency employed, 
which for the commutators in Gish-Rooney outfits is approximately 
16, and nearly 50 in Meggers. 

The only way to avoid polarization and induction effects is by 
using non-polarized electrodes, which consist of cups filled with copper 
sulphate and provided with a permeable bottom, made either of 
porous clay or of wood. 

The induction between leads is minimized when using the po- 
tential measuring line at right angles to the current basis, as already 
described, 
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There is no reason why alternating current could not be used in 
resistivity work, the same as it is done in the potential-drop-ratio 
method. For instance, the current base could be supplied with A.C. 
of moderate frequency by a buzzer, or vacuum-tube oscillator, or by 


Fic. 9.—McCollum earth-current meter (after McCollum). 
EARTH-CURRENT METER 


5, milliammeter; 6, voltmeter; 9, commutator switch; 10, short-circuiting switch; 11, 
voltmeter switch; 12, commutator crank; 13, ammeter switch; 14, resistance switch. 


EARTH-CURRENT METER CONTACTORS 


1, trench contactor; 2, cantilever contactor; 3, single terminal electrode with extension 
rod; 4, single terminal electrode (disassembled). 


taking, if available, current from the light circuit, and feeding it to 
the ground through a transformer. The only difficulty would be to 
obtain an accurate reading of the current; however, the absolute value 
of the current is not of great importance when using the methods JJ, 
II a, or II b. To measure the potential difference between the potential 
electrodes, it would be feasible to have a compensator, and a lead to 
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a pickup coil coupled to the generator or the power leads; the object 
of this lead being to carry the reference e.m.f. and the generator phase 
to the point of observation. 


b. ORIGINAL GISH-ROONEY APPARATUS 


This apparatus has been described so often in the literature that 
it is hardly necessary to gointo more details. It contains a potentiom- 
eter for measuring the potential difference between the two potential 
electrodes, a milliammeter for measuring the current, and a double 


Fic. 10a.—First model of CSM Gish-Rocney outfit, built after MSM model. Hand- 
cranked commutator box which contains also batteries, shown in center; milliammeter 
at left and potentiometer at right. 


commutator, the purpose of which is to pass current to the ground 
in alternating direction, yet pass the current and the potential to be 
measured always in the same direction through the meters. Figure 8 
shows a sectional view of the apparatus. 


¢. MC COLLUM EARTH-CURRENT METER 


Although this instrument was developed altogether independently 
from the Gish-Rooney apparatus and for a different purpose, it is, 
nevertheless, of practically the same construction as the Gish-Rooney 
apparatus. It is shown in Figure 9. It has a voltmeter and milliam- 
meter like the Gish-Rooney apparatus, and also a commutator. It 
differs somewhat from the Gish-Rooney as a number of ranges are 
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provided for both the milliammeter and the voltmeter, and as non- 
polarized electrodes are provided (reference list No. Is). 


d. MODIFIED GISH-ROONEY APPARATUS 


A number of oil and mining companies have designed their own 
Gish-Rooney equipment, which differs slightly from the original in 
one way or another. The Michigan School of Mines has obtained 
much experience in the design of this equipment, and their apparatus 


Fic. 106.—Construction of commutator in instrument shown in Figure toa. 


was used as model for one of the instruments built in the instrument 
shop of the Colorado School of Mines (Fig. 10 a and Fig. 10 5). 

Figure 11 shows the latest form of a Gish-Rooney equipment 
built at the Colorado School of Mines. On the left side of the figure, a 
box is shown containing the potentiometer, milliammeter, switches 
and jacks, and below the panel, the commutator, which is driven by 
a small D.C. motor, operated from a storage battery. The battery 
box is seen between the reels and the meter box. 


é€. MEGOHMERS 


These instruments have been used for years for the testing of 
grounds and‘insulation in electrical power-plant-, railway-, and radio- 
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engineering. They differ from the Gish-Rooney equipment chiefly in 
that the potentiometer and the milliammeter are combined in one 


Fic. 11.—Latest CSM model of Gish-Rooney equipment with motor-driven com- 
mutator. 


: 


Fic. 12.—View and circuit diagram of Megger manufactured by Biddle. 


instrument, the ohmmeter, the moving coil of which has two wind- 
ings, the voltage and the current winding. They also differ from the 
Gish-Rooney equipment as the power is not taken from batteries, 
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but from a hand-driven D.C. generator, on the shaft of which a double 
commutator, similar to the one used in the Gish-Rooney apparatus, 
is mounted. 

The most important manufacturers of Megohmers or Meggers 
are: Evershed and Vignoles, Acton, England; James G. Biddle, Phila- 
delphia; and Herman H. Sticht and Company, New York. 

The accompanying figures show a number of Meggers, which are 
all very similar in principle. Details on the construction of these in- 
struments may be obtained by writing to the manufacturers and ask- 
ing for their cata)ogs. Figure 12 shows wiring diagram and sectional 
view of a Megger manufactured by Biddle, and Figures 13 a and 13 b 
show Meggers manufactured by the Sticht Company. 


III. INTERPRETATION OF RESULTS 


The fundamental principle underlying the interpretation of the 
resistivity results is comparatively simple. 

The depth reached depends on the spacing of the potential electrodes 
(or their distance from a current electrode in methods where station- 
ary current electrodes are used). 

In other words, if the observed resistivity values are plotted as a 
function of the electrode separation (in the 4-terminal method), or 
as a function of the distance of the potential electrodes from the near 
current electrode in single-probe methods, then the electrode separation 
or electrode distance at which marked changes in apparent resistivities 
occur is approximately equal to the depth to the resistivity discon- 
tinuities underground. 

This statement is made merely to illustrate the nature of the 
fundamental principle underlying resistivity measurements. In prac- 
tice, however, there are more difficulties, as the presence of more 
than one discontinuity of resistivity is considerably complicating. 

Actual experience shows that this simple depth rule works sur- 
prisingly well in mining when insulated bodies of good conductivity 
are encountered. In application to stratified ground, severe complica- 
tions arise, chiefly with the Gish-Rooney method, and particularly, 
when there is more than one layer present. 

In the practice of resistivity work, therefore, we use several lines 
of attack in the interpretation of results. As in other types of geo- 
physical work, both qualitative and quantitative methods may be 
applied. The qualitative method usually precedes the quantitative 
method in application, and the quantitative method is applied if the 
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Fic. 13a.—Model D M Megohmer, manufactured by Sticht. 


Fic. 13b.—Paragon Megohmer, manufactured by Sticht. 
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nature of the indications, derived from their qualitative analysis 
warrants it. 

Due to the nature of resistivity work, qualitative methods are 
usually applied in resistivity mapping. The area is covered with a 
4-terminal contacting outfit, and resistivities are measured; the re- 
sults are plotted and points of equal resistivity are connected by lines 
which resemble contour lines, and which are interpreted accordingly. 
The interpretation assumes a somewhat more quantitative nature if, 
at every point, instead of one electrode spacing, two electrode spacings 
are used, so that two contour maps may be drawn for any one area, 
covering two depths of penetration. Finally, considerations of a 
quantitative nature enter also when the contact plane of two forma- 
tions of different resistivities is considered, and if one, two, or three 
of the contacts are in one, and the remainder in the second formation. 
The theoretical computations for this case have been carried out 
by Tagg (ref. list No. IIo) and Hedstrom (ref. list No. IVz) and the 
results are shown in Figure 14. The curves illustrate how the apparent 
resistivity is influenced, for various values of resistivity ratios of the 
formations involved, by the distance of the center of the contacting 
arrangement from the boundary. The curves of Figure 14 are com- 
puted for four electrodes at right angles to the formation boundary, 
and show, therefore, four discontinuities. If the contacting arrange- 
ment is parallel with the fault plane, the apparent resistivity changes 
gradually as the boundary is approached. Tagg has also computed 
the apparent resistivity curves for this case (ref. list No. Iie). 

Outside of the conditions occurring in the immediate vicinity of 
a formation boundary, the method of resistivity mapping involves 
essentially only qualitative methods of interpretation, which consist, 
as stated before, of an interpretation of the equi-resistivity lines 
(Figs. 28 and 29), or of a resistivity profile (Figs. 26 and 30). In inter- 
preting a resistivity profile, it is important to take into consideration 
the influence of the electrode separation, and thus of the depth pene- 
tration, for which Figure 26 is an example. Finally, in the interpreta- 
tion of results obtained by resistivity mapping, tank experiments may 
prove to be of great value, inasmuch as they are very easily performed. 
The geologic bodies are usually represented by metallic conductors, 
and these are arranged at various depths or with varying angles of 
dip and strike as the case may require. Then the 4-terminal contact- 
ing arrangement is carried across the body at the surface of the water 
or solution in the tank, and the apparent resistivities are measured 
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Fic. 14.—A: Arrangement of 4 electrodes crossing vertical boundary (after Tagg). 
B: Apparent resistivities for various electrode positions (after Tagg). C: Tanks used 
for model resistivity experiments. 
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and plotted. Figure 15 shows, as an example, the results of tank ex- 
periments made at the Colorado School of Mines with the objective 
to study the effect of dip on the equi-resistivity profiles. 

The quantitative methods of interpretation are chiefly used in 
interpreting the results of electrical vertical drilling. The quantitative 
methods may be of a direct and of an indirect nature. In the direct 
methods, we determine the depth to formation boundaries directly 


| | | 


Fic. 15.—Studies of effect of dip of formation models in water tank, using 4- 
terminal contacting arrangement with constant electrode separation of 8 inches. Solid 
lines: Apparent resistivity when line of electrodes is parallel with strike of formation. 
Dotted lines: Apparent resistivity when electrode line is at right angles to strike. 


from the curves. These direct methods are applicable in simple cases, 
that is, when only one formation boundary is present or possibly two. 
Otherwise, indirect methods are much safer; their principle is the as- 
sumption of certain formations in definite depths and a theoretical 
computation of the results to be expected. These results are then 
compared with the data obtained in the field, and the assumptions 
are modified until a satisfactory agreement between theoretical and 
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field data is obtained. In doing so, all available geological data must, 
of course, be considered. It may then happen that the correct inter- 
pretation is the one which does not show as good agreement between 
the theoretical and the field data as another assumption which is less 
probable for geologic reasons. It should again be recalled, in this 
connection, that the possibility of quantitative interpretation in 
electrical vertical drilling depends altogether on whether the forma- 
tions involved retain their resistivities in a horizontal direction. The 
determinations of theoretical effects, such as involved in the indirect 
methods of interpretation, may not only be made by means of com- 
putation, but also by means of model experiments. For this purpose 
tanks may be used, filled with water or weakly electrolytic solutions, 
and the geologic bodies may be represented by metal plates, et cetera 
(for ore bodies, water tables, and faults), or by sand or clay layers. 
When working with metal plates it is convenient to use a horizontal 
tank instead of a vertical tank, that is, a tank in which the model 
ore bodies are moved horizontally away from the electrodes to simu- 
late changes in depth (Fig. 14 C). Instead of using tanks for the model 
experiments, pits may be dug into the ground and filled with such 
alternating layers of clays, sands, et cetera, as the case may require 
(reference list No. III;;).! 

Summing up the methods of interpretation applied in resistivity 
work, we have: 

I. Qualitative interpretation: resistivity mapping 
II. Quantitative interpretation: electrical vertical drilling 
a. Indirect methods 


1. Computations 
2. Model experiments 


b. Direct methods 


There exists a very extensive literature on the theory of quantita- 
tive interpretation in resistivity work. The references to this literature 
are given in section II of the bibliography attached to the end of this 
paper. It is believed to be fairly complete. A perusal of this literature 
will show that practically the only electrode arrangement for which 
a complete theoretical treatment is given is the 4-terminal Gish- 
Rooney method, and the papers which give this theory in probably 
the most complete manner are those by Hummel (ref. list No. Ilh2), 
Tagg (ref. list No. IIjo,14,15), Roman (ref. list No. Il,s), Peters and 

1 If sand and clay layers, et cetera, are used in a small laboratory tank, the latter has 


to be compensated to eliminate the effects of the walls (see T. A. Manhart, ref. list No. 
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Bardeen (ref. list No. IIs), and Ehrenburg and Watson (ref. list No. 

Very little is available on the theory of interpretation of results 
obtained by the single-probe method. It seems, however, that the re- 
sults obtained by this method are generally not nearly as difficult to 
interpret as those obtained by the 4-terminal Gish-Rooney method. 
For the single-probe method in its application, illustrated in Figure 
7, IJa, the thumb rule that the depth of discontinuities of resistivity 
is approximately equal to the arithmetic means of the two potential 
probe distances at which a marked change in resistivity is observed 
generally works well. This holds in particular if the distance PR is 
much smaller than E,P; otherwise, the rule that k=+/2bc should be 
applied. Thus, if Z,P is used equal to PR (method II), the depth of 
the geologic body becomes equal to the electrode separation at which 
its influence is noticed, or h=<a. 

It is beyond the scope of this paper to go deeply into the details 
of the theory of interpretation. However, several diagrams are given 
which represent the results of theoretical computations of the effect 
of vertical changes in resistivity on the 4-terminal method. 

Figure 16 illustrates curves of apparent resistivity for the 2-layer 
case, and for various ratios of the resistivity of the lower layer di- 
vided by the resistivity of the upper layer. The corresponding values 

of & mean the ratio: f 

PrP 
as ordinates, the abscissas being electrode separations; for the pur- 
poses of interpretation, however, it is more instructive to make the 
depth axis the ordinate and plot the observed resistivity values as 
abscissas, in accordance with the statement previously made that, in 
a general way, marked changes in resistivities occur at electrode sepa- 
rations a which are equal to the depth 4 of the resistivity discontinuity 
underneath. (Or else, this depth /# is by a constant ratio greater or 
smaller than a, which is taken into consideration in the diagrams, as 


. The observed resistivities may be plotted 


not h, but = is plotted.) 


Thus, we see from Figure 16 that for all resistivity ratios the 
interface between surface layer and bottom layer coincides approxi- 
mately with the inflection point in the curve. The abscissas, by the 
way, do not represent the observed apparent resistivities directly, 
but the ratio of the observed resistivities to the resistivity of the 
upper layer. For the reason just stated, the ordinates likewise do not 
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give depth directly, but the electrode separation or depth of penetra- 
tion divided by the thickness of the upper layer. As both the re- 
sistivity and the thickness of the upper layer are constants, this 
method of graphical representation does not change the shape of the 
curves. 

It is seen from the figures that in this case the boundaries of the 
formations underneath are reflected as points of maximum curvature, 
or as third derivatives of the apparent resistivity with respect to the 
electrode separation. 

The curve denoted by ‘‘k equal to one”’ is of very great importance 
in practical work; it is frequently encountered in applying resistivity 
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Fic. 16.—Curves of apparent resistivity over surface layer of thickness h for differ- 
ent ratios of resistivities of surface layer and infinite lower layer (after Hummel). 


methods in civil engineering work and related geological problems 
where the thickness of the overburden above bedrock is to be de- 
termined. In such cases the resistivity of the bedrock is ordinarily so 
much greater than the resistivity of the overburden that the curves 
approach the theoretical case where the lower resistivity may be 
considered as infinite. For a comparison with the results obtained in 
practice, see the curves shown in Figure 27. 

The diagrams shown in Figure 17 a and b represent the results of 
the theoretical analysis for the 3-layer case. In Figure 17 a, the second 
medium has the lowest resistivity; in Figure 17 6, the second medium 
has the highest resistivity. It is seen that in this case (3-layer prob- 
lem) the peaks in the curves represent approximately the depth to 
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Fic. 17a.—Curve of apparent resistivity for 
(after Hummel). 


case h=h, and ratio p: pi: p2=1:0.5:2 


Fic. 176.—Curve of apparent resistivity 
=1:3:0 (after Hummel). 


for case of k:=3h and ratio p:pi:p2 
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the center of the second medium. It is of interest to compare the 
theoretical curve of Figure 17 6 with the curve obtained in actual 
practice shown in Figure 34. 

One can readily see that already in the 3-layer case the analysis 
of the results obtained by the 4-terminal method presents consider- 
able difficulties, and it is not surprising that on that account the use 
of this method for electrical vertical drilling has been largely aban- 
doned and has been replaced by the various types of the single-probe 
method and by the potential-drop-ratio method. 

Only in the 2-layer case, a rigorous direct interpretation, and a 
determination of the depth to the interface are possible. How this can 
be done in favorable geological cases, has been shown by Tagg (ref. 
list No. TI0,14,16)- 

Tagg’s method is based on the following consideration. If we 
measure the apparent resistivity as a function of horizontal distance, 
the following quantities are given: the apparent resistivities as a 
function of electrode separation, and the resistivity of the upper 
layer. The ratio of the upper to lower resistivity, and the depth of the 
boundary are unknown. Theoretically, therefore, the equations can 
be solved by determining the apparent resistivities at two electrode 
p’—p 
p’ +p 
by h the depth to the interface, it follows that from the surface re- 
sistivity, and the apparent resistivity for only one electrode separa- 
tion, k as a function of k can be determined and graphically repre- 
sented. By taking then the apparent resistivity for a second electrode 
separation, a second function of # of k becomes available and can 
be represented by a curve, and the intersection of these two curves 
gives the depth. In practice, the more electrode separations used, 
the greater are the number of curves of h=/f(k), and consequently 
the greater the accuracy in the determination of the point of inter- 
section of the corresponding curves. 

In order to show how the interpretation is made, a practical ex- 
ample is demonstrated in Figure 18. In the field, the curve shown in 
the upper part of Figure 18D was obtained, indicating immediately 
(by comparison with Figure 16) that a poor conductor is underlain by 
a good conductor and that & is negative. The resistivity of the top 
stratum is taken from the curve. If the data are not as good as they 
are in the illustrated case, several separate determinations may be 
made with small electrode separation. Then for definite electrode 


separations. Designating by & the resistivity ratio 
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Fic. 18.—A to C: Apparent resistivities as function of depth to boundary and as 
function of k, for separations of 50,100, and 150 feet (after Tagg). D: Application of 
above diagrams in depth determinations. 
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separations, say for 50, 100, and 150 feet (at the points 2, 3, and 4) the 
ratio of the apparent resistivity to the resistivity of the upper medium 
is computed (Fig. 18D). 


Diagrams have been prepared previously for this ratio Pe as a func- 
p 


tion of 4, a being constant for each diagram (Fig. 18 A—-C). (The ap- 
pearance of the curves in these diagrams is much the same as those 
shown in the diagram of Figure 16, except that not a/h, but its 
reciprocal, 4/a is used, and that a is a constant for each diagram.) 
In the example shown in Figure 18D, the ratio of apparent to surface 
resistivity is 0.9 for the separation 50 feet; hence, the lower diagram 
of Figure 184A is used, and the corresponding values of # and k which 
are intersected by the line o.9 are noted and plotted (curve 2 in 
Figure 18D). Then the next point, a= 100, is taken, at which the ratio 
is 0.7; one proceeds as before, this time using the diagram in Figure 
18B, and plotting & as a function of # for the ratio 0.7. This gives 
curve 2 in the interpretation diagram. Curve 3 is obtained in the same 
manner. All three curves intersect at a depth of nearly 58 feet, which 
is the desired depth of the formation boundary. 

Under favorable conditions, this method works very well for the 
2-layer case, and is often applicable in practice when the depth of 
overburden is to be determined. Its application, however, depends 
altogether on an accurate determination of the surface resistivity, 
which is not always easy to accomplish—or else, the resistivities de- 
termined at the surface do not represent at all the resistivity of the 
upper layer, the thickness of which is to be determined. In other 
words, the problem then is no longer a 2-layer problem.':? 

Schlumberger uses the following method of interpretation which 
deviates in a number of respects from that of Tagg’s: 

1) As in Tagg’s method, the resistivity of the upper stratum is 
first determined from measurements with small electrode separations. 

2) The resistivity of the lower layers is then obtained from ex- 
posures or measurements in wells (Schlumberger’s “electrical cor- 
ing’), or from assumptions based on past experience. 

3) With these resistivities, a diagram is constructed showing the 
apparent resistivity as function of electrode separation in several 
curves for a number of values of thickness of the upper stratum. 


1 Compare also the discussion of Tagg’s paper at the meeting of the American 
Institute of Mining and Metallurgical Engineers, February, 1932. 


2 The extension of Tagg’s method to the three-layer case has recently been worked 
out by Manhart (ref. list No. Ij7). 
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4) To find the thickness of the upper stratum, all that is necessary 
is to place the observed resistivity curve upon the one which has been 
computed theoretically in the manner described, and to select the 
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Fic. 19.—Electrode arrangements in potential-drop mapping and potential-drop- 
ratio work. 


depth for which the observed and computed curves show the best 
agreement, 

The same method is applicable to the three-layer case. Poldini 
(ref. list No. IIIs) has recently demonstrated the application of 
this interpretation method in a number of cases. 
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D. PoTENTIAL-DRop-RATIO METHOD 


Although this method has been chiefly developed for the purpose 
of obtaining better interpretable curves in electrical vertical drilling, 
it may, like the resistivity method proper, be applied also in what 
may be termed potential-drop mapping, and phase-anomaly mapping. 

As the name implies, the principle of measurement is the deter- 
mination of the ratios of the potential difference between three points 
(Fig. 19, a). The instruments give the voltage difference A—B di- 
vided by the difference B—C. It is therefore possible, by successive 
occupation of the points A, B, and C, and then of B, C, D et cetera, 
to determine the potential drop along a traverse line. This leads to 
the plotting of potential-drop curves and phase-anomaly curves 
(Fig. 19 6). It is seen that the potential decreases and the phase angle 
increases, above conductive bodies and vice versa.’ However, this 
method of plotting potentials and phases does not furnish nearly as 
much quantitative data in respect to depth of formation boundaries 
as the equi-resistivity map method (for example, Fig. 26). The infor- 
mation obtained by a potential and phase profile is at the present 
stage of interpretative technique, comparable with the information 
conveyed by an equipotential line map. Electrical vertical drilling 
by means of the potential-drop-ratio method, however, furnishes a 
type of curve which is most readily interpreted with reference to the 
depth of formation boundaries (for example, Fig. 35). While the 
mapping of potential drops is chiefly applicable to the detection of 
abrupt changes in conductivity in a horizontal direction’ (steeply 
dipping ore bodies, faults, vertical formation boundaries, and the 
like), the potential-drop-ratio method is to be given preference when 
the depth to horizontal formation boundaries is determined. Figure 
20 shows this advantage of the potential-drop-ratio method by a 
comparison of the indications obtained by it with the type of curves 
furnished by: (a) equipotential line, (0) potential drop, and (d) ap- 
parent resistivity methods. It is seen that the interpretation of the 
potential-drop ratios becomes particularly easy if the curves shown 
under (c) are turned through go°, that is, if the geologic section is 
plotted in the usual manner, with the exception that the vertical scale 
is stretched to 3/2 of the original scale; and the electrical curve is 
superimposed upon it. We shall return to this method again when 
dealing with the interpretation of the potential-drop-ratio indications. 


1 See also Hedstrom (ref. list No. IV7). 
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Fic. 20.—Comparison of potential indications: 
a: Equipotential lines 
b: Potential drop 
c: Potential-drop ratio 
d: Apparent resistivity 
(after Lundberg & Zuschlag) 
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I. ELECTRODE ARRANGEMENTS 
ad. IN POTENTIAL-DROP AND PHASE-ANOMALY MAPPING 


The arrangement customary in this work is shown in Figure 19a. 
Current from an A.C. source, that is, either from a generator or froma 
buzzer, is supplied to the ground at the points E,and £;. The potential 
drops are measured between the power electrodes, at the points A, B, 
C, D, E, et cetera, and the bridge arrangement occupies first the points 
A, B, and C; then the points B, C, and D; then C, D, and E, et cetera. 
The procedure is accurate but laborious. For the reasons just stated, 
that is, that the potential-drop indications are more pronounced for 
vertical formation boundaries than for horizontal ones, this method 
of potential mapping has not found very much application in oil 
work. Only vertical electrical drilling by means of the potential-drop- 
ratio method has gained importance. 


b. ELECTRODE ARRANGEMENTS FOR ELECTRICAL VERTICAL DRILLING 


The two types which have been used to date are illustrated in 
Figure 19 c and d. The electrode set-up used by Koenigsberger differs 
from the arrangement employed by the Swedish American Prospecting 
Company in that the spacing of the current electrodes is changed, 
while in the Swedish American arrangement the current basis re- 
mains constant. Furthermore, in Koenigsberger’s arrangement the 
current basis is short and only twice as long as the interval of the 
search electrodes A, B, and C while in the Swedish method the ex- 
ploring electrodes are used in the vicinity of one current electrode and 
the other current electrode is practically at infinity. Finally, in the 
Swedish American arrangement the interval between the exploring 
electrodes remains fixed, while in Koenigsberger’s method this interval 
increases with an increase in the length of the current basis. Hence, 
for purposes of interpretation, the depth reached in Koenigsberger’s 
method is a function of the interval a, while in the Swedish method 
the depth is proportional to the distance of the center probe, B, 
from the power electrode E. In Koenigsberger’s method the measure- 
ments are made in the extension of the current basis; in the Swedish 
American method, the traverses may be laid so as to radiate in all 
directions from the near power electrode; toeliminate the influence of 
the far current electrode completely, the most advantageous arrange- 
ment is at right angles to the current basis (strictly on the circum- 
ference of a circle with the radius £,E_ about the far current elec- 
trode). 
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II. APPARATUS 
a. GENERAL: POWER SOURCES 


Potential-drop-ratio measurements may be made with both 
D.C. and A.C. Koenigsberger has used both types of current, while 
the Swedish American employs only A.C. If D.C. is used, a number 
of B-batteries are sufficient, if the required depth penetration is not 
great. Otherwise, a gasoline engine D.C. generator is more advanta- 
geous for greater depth of exploration. The current should be inter- 
rupted by a commutator similar to the one employed in the Gish- 
Rooney apparatus, so that effects of polarization and stray currents 
are eliminated. As will be stated later, the Gish-Rooney equipment 
itself may be employed for potential-drop-ratio determinations (with 
slight modification) for simple problems, although the accuracy 
furnished by the potentiometer customarily employed in this ap- 
paratus is not sufficient for many purposes. 

For the potential ratiometers as used by the Imperial Geophysical 
Survey and the Swedish American Prospecting Company, alternat- 
ing current of about 500 cycles is required, which for small depths of 
penetration is furnished by a buzzer operated from a storage battery 
(Fig. 23) while for greater distances and unfavorable contact con- 
ditions, a 500 A.C. generator is preferable. It must be borne in mind, 
however, that the use of a stronger power source alone with average 
frequencies does not bring about a greater depth penetration, as the 
latter decreases with the frequency employed; hence, for great depth 
penetrations, the frequency has to be considerably lowered, which 
then also would require a change in the dimensions of resistances and 
inductances in the ratiometers, or else commutated D.C. should be 
used. 


b. KOENIGSBERGER’S AND OTHER D.C. METHODS 


Koenigsberger indicates in his publication (ref. list No. IV) in a 
general way a number of arrangements which he has used for the 
measurement of potential-drop ratios, but does not go into details 
with reference to apparatus employed. When using D.C., the potential 
differences A—B and B—C may, of course, be measured independ- 
ently by a potentiometer and their ratio may be formed by computa- 
tion. For this purpose, a regular Gish-Rooney equipment may be 
employed, with an electrode arrangement identical with that used 
in the single-probe method (Fig. 7, type J7a). The potential differ- 
ence between the points B and C is measured in the usual way; then 
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the electrode in C is left there and the electrode which was at B is 
moved to a next point, say D, which has the same distance from C, as 
B has from C. The curves shown in Figure 33 have been obtained in 
this manner. The method is not very accurate, as the surface and con- 
tact resistances enter considerably. 
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Fic. 21 
a: Detailed circuit diagram of Racom compensator 
b: Schematic circuit diagram of same 


Koenigsberger uses also a modification of this method: The central 
electrode is left permanently connected to the potentiometer and the 
other electrodes, B and C, are connected alternately with it by means Ly 
of a commutator. The ratio of the deflections gives the potential 
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drop ratio; or else, it is possible to leave one electrode, say C, station- 
ary, and change the position of the other until a balance is obtained. 

A still better arrangement would be the bridge shown schemati- 
cally in Figure 216, when modified for direct current: the variometer 
V would be left out, and the telephone T would be replaced by a sen- 
sitive galvanometer.! To overcome the contact resistance at the 
probes, two settings of the resistances R, and Rz would be required. 
If commutated D.C. is used, the commutator could be incorporated 
in the bridge. This would necessitate a lead from the bridge to the 
power source. In order to eliminate this lead, very low frequency 
A.C. could be supplied to the power stakes, and a rectifier could be 
used with the galvanometer. 

In addition to D.C., Koenigsberger has also used A.C. with 
bridges similar to the arrangements now to be described. 


b. SWEDISH AMERICAN RACOM 


The name Racom is an abbreviation for ratio compensator. The 
principle of the instrument is a compensation of the difference of the 
potential P4_z and Ps_c between the three points A, B, and C, by 
adjusting resistances in the two circuits AB and BC to balance 
on a telephone. Alternating current of about 500 cycles is used, and 
for small depths is furnished by a buzzer, operated from a storage 
battery, shown in Figure 23. A schematic circuit diagram is given 
in Figure 21, 6; as seen from the figure, the current in the circuit AB 


is equal to Ae and the current in the circuit BC is equal to 
RatR, 


-ee where Ry, and R¢ are the contact resistances of the stakes 
A and C. As the resistances R; and R2 are so adjusted that the tele- 
phone is silent, we obtain the equation for the potential-drop ratio 
Pa-p_RatR, 
Pec 
second setting of R; and Ry. The phase differences are compensated 
by a suitable setting of the variometer, V. If the observed currents are 
too feeble to be heard on the telephone an amplifier is employed, as 
shown in Figure 22. A schematic diagram and a detailed circuit dia- 
gram are given in this figure. The amplifier is a 2-stage, resistance- 


. The contact resistances may be eliminated by a 


1 A similar arrangement is used for pipe-line corrosion studies by Schlumberger 
(ref. list No. III). 
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Fic. 22.—a: Wiring diagram, Racom amplifier. Courtesy Swedish American Pros- 
pecting Corporation. b: Schematic diagram. 


Fic. 23.—View of Racom with amplifier and buzzer. 
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impedance-coupled amplifier, incorporating a tuning of the choke for 
the elimination of harmonics so as to enable a better setting of the 
bridge. 

Figure 23 shows a view of the Racom, with amplifier underneath 
the instrument box, and a buzzer on the side together with the power 
and the exploring electrodes. 

For more details on the operation of the Racom and examples of 
a number of surveys, see Lundberg and Zuschlag’s publication on 
the subject (ref. list No. IVs). 
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Fic. 24.—I,E.G,S. ratio arm bridges. 
Fic. 24a.—Series ratio arm bridge. Fic. 246.—Parallel ratio arm bridge. 


d. 1.G.E.S. RATIO COMPENSATORS 


The potential ratiometers used by the Imperial Geophysical Ex- 
perimental Survey in Australia are very similar in design to the 
Swedish American bridge just described. A schematized view of their 
bridge is given in Figure 19 a. It is seen that this bridge differs from 
the Racom type only by the fact that condensers are used instead of 
inductances for the determination of phase shifts. The wiring dia- 
grams for a series and a parallel ratio arm bridge are shown in Figure 
24. Figure 25 gives a view of the panel of one of these bridges. An 
amplifier is also used in conjunction with these bridges, and in the 
field is carried by means of straps on the operator’s back. For details 
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on the operation of these ratiometers and their theory, the publica- 
tion by Edge and Laby (pp. 50-54, and 268-274) (ref. list No. IVs) 
should be consulted. 


POTENTIAL ®RATIOMETER% 


reo 


Fic. 25.—View of I.E.G.S. ratio compensator. 


III. INTERPRETATION 


Methods of interpretation of potential-drop-ratio results have 
been worked out in a general way for Koenigsberger’s and for the 
Swedish American methods. 

The results of the computations show that Koenigsberger’s 
method is, in a number of ways, at a disadvantage as compared with 
the Swedish method, because his ratios are only opposite but not 
equal for good and poor conductors which have the same conductivity 
ratio. What is meant by this is the following: If we form, for a 2-layer 
problem, the ratio of the conductivity of the surface layer to the con- 
ductivity of the lower layer, then the potential-drop ratio in Koenigs- 
berger’s method is 3.57 if the conductivity ratio is one (homogeneous 
ground). However, if the conductivity ratio is 10:1, the potential 
ratio decreases only to 2.48, whereas, in the converse case (if the con- 
ductivity ratio is 1:10) the potential ratiois 8.70 (Fig. 40 c). Thus, it 
is seen that in his method the potential-drop-ratio responses are not 
symmetrical around the 3.57 position, and that the influence of good 
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conductors below is much greater than the influence of poor con- 
ductors below. This would make this method a rather valuable ac- 
cessory to the Racom, if it is desirable to emphasize in the results 
the effects of good conductors, such as ore bodies or formations 
which are filled with salt waters. 

Figures which give the general values for the potential ratios in 
the Racom method have not been published, but it appears from the 
diagrams given by Lundberg (Fig. 20) that for the average con- 
ductivity ratios the potential-drop-ratio responses are approximately 
the same for a good conductor below as they are for a poor conductor 
below; for, as seen from Figure 20, for the assumed conductivity ratio, 
the potential ratio is 1.5 for a good conductor above and o.5 for a 
good conductor below; or, the deflections in the Racom method are 
equal and opposite for equal and opposite conductivity ratios. Hence, 
when plotting in the Racom method the potential ratios on’a vertical 
axis and the formation boundaries in a horizontal position, a poor 
conductor manifests itself by a deflection of the ratio curve toward 
the right, and a good conductor appears as a deflection of the ratio 
curve toward the left. For a poor conductor below, the ratio is greater 
than one, and for a good conductor below, the ratio is smaller than 
one. 

The chief advantage of the potential-drop-ratio method lies in the 
ease with which the potential ratio indications, if plotted in a suitable 
manner, may be interpreted with reference to the depth of formation 
boundaries. Although the potential-drop ratios, and hence also the 
relations between depth of boundary and maximum of indication, 
depend on the conductivity ratios of the formations involved, yet 
several empirical rules have been found to hold very well for the 
average conductivity ratios encountered in the field. 

In Koenigsberger’s method, the depth corresponding with a ratio 
maximum is approximately 14 a for conductivity ratios up to 1:10 
or 10:1. In the Racom method, the depth to the formation boundary 
is, under similar conditions, equal to 24 r. (The significance of the 
symbols is explained in Figure 19.) Hence, when we turn the ratio 
traverses through go° and plot the curves vertically on a geologic 
section, we have to stretch the vertical scale of the section to 3/2 in 
the Racom, and double it in Koenigsberger’s method. Or, if we want to 
leave the scale of the geologic section as it is, we have to reduce the 
distances in Koenigsberger’s method to 4 and the distances in the 
Racom method to 3 of their original values. These relations are well 
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explained in the practical examples illustrated in Figures 33, 35 and 
40. 
E. OPINIONS REGARDING DrIrREcT ELECTRICAL 
PROSPECTING FOR OIL 


In view of the recent developments of the technique of resistivity 
and potential-drop-ratio-measurements which permit (as has been 
shown and can be verified by an inspection, for example, of Figures 
34 and 35) the accurate identification of certain geologic horizons 
under favorable circumstances, the question becomes, of course, acute 
again as to the extent to which oil deposits may be located directly by 
electrical measurement from the surface. 

Very much has been written about this subject, probably altogeth- 
er too much at a time when field technique and interpretation of elec- 
trical prospecting were not nearly far enough developed to enable 
anyone involved in the discussion to state definitely that certain 
things were altogether impossible. The writer believes that several 
authors who were involved in the animated discussion at that time 
wish now, in view of the recent developments, that they had not 
voiced their opinion in the matter. 

Another factor which influenced the exchange of opinions in the 
problem of the possibility of locating oil directly was the fact that 
most of the involved parties were interested in competitive geo- 
physical consulting companies. Hence, the writer doubts if very 
much could be gained by repeating in detail at this time all the argu- 
ments which have been advanced for and against the direct electrical 
location of oil. Only the main points will be discussed. At any rate, 
it is emphasized that it is not the object of this paper to take sides in 
an argument which is not as yet ripe for an altogether unbiased dis- 
cussion. We shall not try to decide at this time if a direct location of 
oil by electrical prospecting is possible. We shall merely state, by 
discussing actual results obtained, what has been done and can prob- 
ably be done with the present state of development of instrumental 
technique and interpretation. The whole discussion in this matter 
seems to have been started when the representatives of the Elbof 
electrical prospecting company claimed that they were in a position, 
by means of their method, to locate the oil directly. The chief factor 
which probably encouraged doubts and discussions of their statements 
was the fact that the type of indications which they obtained by 
means of their method, namely the deflection of current lines, did not 
lend itself very readily to such far-reaching interpretations. In fact, 
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even under favorable geologic conditions, the deflections of the cur- 
rent lines were usually small and therefore the claimed effects could 
be readily disputed. In addition, this company made their claims 
on the basis of some work which was carried out at comparatively 
great depth and under complex geologic conditions, so that it was 
easy for their opponents to go so far as to state that there was no 
telling whether the observed deflections came from salt water or from 
oil. 

The opponents of the Elbof theory, notably Ambronn, then went 
to the other extreme in trying to prove that it was theoretically im- 
possible, under any circumstances, to detect oil directly with alternat- 
ing current. Ambronn claimed that the conductivity ratio required 
for electrical prospecting by means of alternating current should be 
at least 1,000 for such a change to be detectable, a figure which has 
been proved to be much too high by actual experience obtained since 
with the potential and resistivity methods. In addition, Heine and 
Hunkel proved theoretically that the formulas used by Ambronn for 
his proofs were in error. 

Lately, the Elbof company has also taken up, in addition to their 
electromagnetic method, the study of the potential-drop-ratio and 
resistivity methods; the results obtained in several areas have been 
published by Gella and Koenigsberger (ref. list III:3; [V2). Unfor- 
tunately, in both publications, no curves are given and this is prob- 
ably the reason why these more recent claims of the Elbof company 
and Koenigsberger, which are based on reliable and favorable tech- 
nique, have passed almost unnoticed. The figures given for the re- 
sults obtained with the potential drop ratio at the oil field of Oberg 
in Germany are given in such manner that the reader who is not 
familiar with the locality can not form an opinion as to the extent to 
which the claims advanced are justified. The writer has made an at- 
tempt to use as many of Koenigsberger’s data as possible and to 
superimpose them upon the geologic section; the results are repre- 
sented in Figure 40 and are discussed at the end of this paper. 

While working with the potential-difference method, the Elbof 
company found that in certain areas where the oil-bearing strata are 
overlain by porous strata such as sands and sandstones, the gas had 
in many places migrated into them from the strata below, had dis- 
placed the salt water and had thus increased considerably the thick- 
ness of the strata acting as poor conductors. This may make it possible 
in some places to apply the electrical method to advantage, but it is, 
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on the other hand, a condition which is by no means representative of 
all oil reservoirs. 

A rather serious argument has been advanced lately by Hedstrom 
(ref. list No. IIIg) against the direct electrical prospecting for oil. 
This author claims that from theory and experiment it follows that a 
formation which is practically an insulator does not show in the 
electrical results any differently from what it would if its conductivity 
were only 1o times smaller than that of the surrounding formation; 
hence, a poorly conductive sand or limestone with a conductivity 
only about ro times less as compared with that of the surrounding 
strata, would not be distinguishable from the effect of an oil forma- 
tion. 

This argument, however, does not preclude the possibility of 
direct location of oil in an area where, by correlation with geologic 
data, a poor conductor is known to be oil bearing, and where it would 
thus be possible to trace this oil-bearing horizon in various depths 
of structure, and to determine where the filling of the pores with oil 
ceases and where the oil is replaced by salt water. If the argument 
just advanced is correct, it would seldom be possible to distinguish 
between a sand filled with gas or a dry sand, and a sand filled with 
oil. 

One might finally say in commenting on Hedstrom’s argument, 
that with a greater perfection of technique of electrical prospecting 
we may some day be able to distinguish the surface effects of poor 
conductors with a greater precision and obtain an electrical log in 
which the formation resistivities could be represented by their true 
and not their apparent resistivity, so that an oil sand would show 
somewhat like the sands illustrated in the electrical logs demon- 
strated in Figure 37 and Figure 38. 

However, as already stated, it is beyond the scope of this paper 
to take sides in this argument, and we will merely confine ourselves 
to a discussion of the results thus far obtained in attempting to locate 
the oil directly at the end of the following section (F II). 


F. RESULTS OBTAINED WITH RESISTIVITY AND 
POTENTIAL-DRop-RATIO METHODS 
I. IN STRUCTURAL WORK 


Several outstanding examples are here described which illustrate 
both the possibilities and the limitations of the resistivity and po- 
tential-drop-ratio methods in determining geologic structure. 


C. A. HEILAND 


@, DETERMINATION OF THICKNESS AND CHARACTER 
OF SURFACE FORMATIONS 

The chief potentialities of both the resistivity and the potential- 
drop-ratio methods lie in the determination of character and depth to 
formation boundaries of surface formations. In general, as in any other 
geophysical method, the thickness of geologic formations has to 
compare favorably with their depth, otherwise they may pass un- 
noticed in resistivity indications; on the other hand, conditions for 
recognizing formation boundaries are much more favorable in the 
potential-drop-ratio method than they are in the resistivity methods. 
In addition, a reason why these methods lend themselves so readily 
to the determinations of depths and thicknesses of surface formations 
is the fact that they often differ much more in their resistivity from 
the underlying series than these lower formations will differ between 
themselves. 

A problem which often arises in practical geology is the determina- 
tion of the thickness of glacial overburden. Both the method of re- 
sistivity mapping and the method of vertical electrical drilling can 
be applied to advantage, depending on the problem. 

Thus far, the majority of examples for the determination of the 
thickness of surface formations have become known from applications 
in mining and civil engineering, but they would also be of importance 
in certain types of oil-structural work. 


I. BY EQUI-RESISTIVITY MAPPING 

A rather convenient way to map bedrock contours in a fairly 
large area, or to determine the depth to bedrock along a traverse 
line, is to carry a 4-terminal contacting arrangement with fixed elec- 
trode separation over the area and to observe and plot apparent re- 
sistivities in the manner shown in Figure 26. The results obtained with 
two electrode separations are illustrated. The solid contour represents 
the depth to bedrock, the dashed line is the indication obtained with 
the correct electrode separation, and the dotted lines show the ap- 
parent resistivities with too small electrode spacing. The geologic 
section was glacial drift above limestones and conglomerates. This 
figure shows that with an electrode separation larger than the greatest 
depth to bedrock, the bedrock contour may be mapped with a good 


deal of accuracy. 
2. BY ELECTRICAL VERTICAL DRILLING 


This method has been applied very frequently in late years in 
determining the depth to bedrock in tunnel and dam sites. Figure 27 
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Fic. 26.—Resistivity mapping with two electrode separations to determine depth 
of cover (after Lancaster-Jones). 
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Fic. 27.—Results of electrical vertical drilling to determine thickness of cover 
(after Leonardon). 


A: Diagram in case of fairly homogeneous overburden. Accuracy fairly satisfactory. 
B: Form of results when overburden is composed of different layers. 
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shows two resistivity profiles with the corresponding drilling results. 
In interpreting the resistivity curves, reference should be made to 
the curves shown in Figure 16 for the case in which the resistivity of 
the lower layer is practically infinite. The two curves shown in Figure 
27 also illustrate the effect of a homogeneous and an inhomogeneous 


overburden. 
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Fic. 28.—Resistivity map of the Aricesti anticline. Work carried out in 1923 
(after Schlumberger). 


3. BY POTENTIAL AND PHASE MAPPING 


This method will probably find little application in oil structural 
work, although it can give valuable indications in determining the 
character of surface formations in suitable areas under a shallow 
cover. Figure 19 shows, as an example, the effect of a gabbro intrusion 


with two ore bodies on the contact zone upon the potential-drop and 
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Fic. 29.—Resistivity map of fault in Triassic in Lorraine. Study made in 1928 
(after Schlumberger). 
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phase curves. It will be noticed that the potential decreases above a 
good conductor, while the phase angle has a maximum, and vice versa. 


b. LOCATION OF STRUCTURE 


Both resistivity and potential-drop-ratio methods have found 
already several remarkable applications in locating geologic structure 
favorable for the accumulation of oil, such as anticlines, salt domes, and 
the like. They have also been applied very successfully in extend- 
ing known geologic information into unknown territory, which is an 
application always recommended for any kind of geophysical method. 
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Fic. 30.—Resistivity profile across buried anticline, taken with electrode spacing a 
equal to 200 feet (after Hubbert). 


I. BY RESISTIVITY MAPPING 


The Schlumberger Company has located and mapped geologic 
structure of several types by the resistivity method. As examples, the 
anticline illustrated in Figure 28 and the fault illustrated in Figure 29 
have been selected. The first figure shows an equi-resistivity map of the 
Aricesti anticline in Roumania which appears in the resistivity map 
as a body of lower resistivity. Very instructive is the location of a 
fault in Alsace-Lorraine by the method of resistivity mapping, as it 
illustrates how definite resistivity values may be identified with 
definite geologic horizons. The marls have the lowest resistivity, then 
follow the limestones, and then the sandstones in this area. The re- 
sistivity contours follow truly the trace of the fault which has dis- 
placed the formations under discussion in a horizontal direction. 
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Figure 30 shows another good example of structure located by 
means of a resistivity profile. The traverse was made with a 4-terminal 
contacting arrangement using an electrode interval of 200 feet. As this 
interval is greater than the greatest thickness of the glacial till in the 
area, the observed apparent resistivities represent chiefly the varia- 
tions in resistivities below the glacial drift. Thus it was possible to 
get the effect of the anticline composed of Devonian limestone of high 
resistivity, flanked by Mississippian shales and sandstones of lesser 
resistivity. Hubbert states that the causes of the smaller variations 
in the resistivity peak are unknown; it is probable, howe' ever, that they 
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Fic. 31.—Results of electrical vertical drilling with single-probe method (II) and 


corresponding geologic section in an area in New Mexico. Courtesy of Harry Aurand, 
Midwest Refining Company. 


are either due to variations in the composition of the glacial till, or 
due to the fact that this till occurs as filling in erosion channels on 
top of the structure. 

2. BY ELECTRICAL VERTICAL DRILLING 

There follows now a discussion of some examples of structural 
resistivity work, most of which were obtained recently and may not 
be found elsewhere in the geophysical literature. They illustrate well 
the great possibilities of the method of electrical drilling with the 
resistivity and potential-drop-ratio method under favorable con- 
ditions. 

Figure 31 shows the results of several resistivity traverses which 
were very carefully measured every 1o feet by means of the single- 
probe method (No. II of Hummel in Figure 7 or the Ehrenburg- 
Watson method). The geologic section consists at the surface of al- 
ternating shale, gypsum, and limestone beds with an average thick- 
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ness of 400 feet underlain by salt. The resistance of the gypsum is 
high, and that of the shales is low. Where the effect of the salt has 
been picked up, its resistance seems to be fairly low, such as shown 
in well No. 48, which is somewhat contrary to the experience obtained 
on salt domes. All resistivity traverses shown in Figure 31 have a simi- 
lar appearance as far as the very high resistance surface layers are con- 


haddaty Grp 42 
3 
apparen t 


Resistivity 


Single Probe 
Method 


= 
Water Sand 


Fic. 32.—Results of electrical vertical drilling using single-probe method (IIa) and 
corresponding well log at Berthoud dome, Colorado. (Measurements made by Keen, 
Bishop, and Johnson.) 


cerned. Then they all show a pronounced “low,” corresponding with 
the influence of the red beds and shales close to the surface. After that, 
all curves go up again, depending on the amount of gypsum present 
in the section. Thus, while the curves of Nos. 45 and 46 are practically 
flat, they are steeper in 49, 44, 43, and steepest in 48. Just what takes 
place in No. 48 is not very well known on account of lack of knowledge 
of the well log; but it seems that the formations present at that local- 
ity, probably gypsum, are very poor conductors. A striking similarity 
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prevails between the curves of 43, 44, 45, and 46. Particularly, 43 is 
very similar to 44, and 45 very similar to 46. Corresponding with the 
bottom of the first red-bed layer, all the curves mentioned have the 
same appearance, which becomes less pronounced toward the north- 
east, but may yet be recognized very faintly in curve 49. Therefore, 
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Fic. 33.—Results of potential-drop-ratio observations and corresponding structure 
at Berthoud dome, Colorado. 


the writer has used that particular part of the curve as horizon 
marker for the resistivity correlation. It is seen that the resistivity 
contour thus obtained is not everywhere in accord with the logs, but 
that it checks very well with the structure correlation made on top 
of the salt beds. The next two figures, 32 and 33, show some results 
which were obtained with the single-probe and the potential-drop-ratio 
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methods on the Berthoud dome, not far from Denver. Some oil and 
gas have been produced from this dome. The geologic section consists 
of Cretaceous shales and sands. The curve shown in Figure 32 was ob- 
tained by using the method No. JJa (Fig. 7), and is plotted together 
with the log of a well next to which the measurements were made. 
Pending further work in this area, the present interpyetation is that 
the peaks in the curve are produced by the water sands, and that these 
sands contain water of artesian origin, which is encountered in 
geologic structure of many types in eastern Colorado near the Front 
Range of the Rocky Mountains. If nothing else, this curve shows that 
certain resistivity peaks appear at certain depths which can be cor- 
related with definite geologic horizons, and which can be traced 
throughout the structure. 

This makes it possible to determine structural contours in this 
area as shown in Figure 33. It should be understood that these results 
were obtained with an equipment which was not at all adequate for 
the purpose. Before the regular Racom equipment (already de- 
scribed) became available for field work, we made some preliminary 
field tests with the Gish-Rooney apparatus modified for the deter- 
mination of potential-drop ratios. Hence, the results are masked to 
some extent by the variations in surface resistivities which can be 
eliminated to a much greater degree in the Racom. Nevertheless, the 
results are reliable enough to justify at least a preliminary inter- 
pretation. The potential-drop-ratio curves were turned go°® and 
plotted vertically. It is seen that certain peaks in the curves repeat 
themselves in every traverse, and these are indicated in solid black. 
Then the structural contour as known from the well logs was super- 
imposed upon the potential-drop-ratio profiles, and was plotted on a 
scale 3/2 that of the potential profiles, in accordance with what has 
been previously stated about the interpretation of potential-drop- 
ratio measurements. It is seen that the known structural informa- 
tion and the peaks in the potential-drop-ratio profiles correspond 
fairly well. In this connection, it should be remembered that the in- 
dications obtained from a potential-drop-ratio log in electrical verti- 
cal drilling are altogether different from the indications obtained in 
a resistivity log. In the former, formation boundaries are marked by 
a peak in the drop-ratio curve (Fig. 20). 

Although the results just described are not as good as they might 
have been if more adequate equipment had been used they proved 
definitely that the potential-drop-ratio method could be used to ad- 
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vantage in eastern Colorado, in close conjunction with geologic data, 
for the purpose of obtaining data about geologic structure. 

Figure 34 shows again some resistivity curves, obtained on four 
traverses in the San Fernando Valley in California. The results of 
this work were published by J. J. Jakosky (ref. list IIIs:) in different 
form; it appears from the description of the work that some sort of a 
single-power-probe method was used. It may be possible that the 5- 
electrode system described by Jakosky' in a later paper was employed, 
on which no data are at present available, and which may be some sort 
of potential-drop-ratio method with an electrode arrangement similar 
to that shown in Figure 19 or may be the partition method (Fig. 7, Ia). 
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Fic. 34.—Results of four resistivity traverses (single-probe method) and corre- 
sponding geologic section in San Fernando Valley, California. (Plotted from data pub- 
lished by J. J. Jakosky.) 


As stated by Lundberg and Zuschlag (ref. list No. [Vs) apparent re- 
sistivities may also be readily computed from the potential-drop ratios. 

In Figure 34, the results obtained by Jakosky are shown, plotted 
in millivolt units which are directly proportional to the apparent re- 
sistivity. The geologic conditions are as follows. The larger part of 
the area is covered by a blanket of alluvial strata with a maximum 
thickness of nearly 80 feet. Below, in the eastern part of the area, are 
inclined strata of the Fernando group of the Pliocene. These forma- 
tions are cut off by a fault in the western part of the area which has 
brought strata of the Cretaceous Chico formation in contact with the 
Pliocene in the east. These Cretaceous strata do not contain water, 
but the upper parts of the Pliocene are filled with highly conductive 
waters down to a depth where the impervious strata indicated as 


1 J. J. Jakosky, “Geophysical Examination of Meteor Crater, Arizona,” Geophysi- 
cal Prospecting (Amer. Inst. Min. Met. Eng., 1932), pp. 63-98. 
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such in the figure occur. As the Pliocene strata are cut off in the west 
by the fault, a conductive triangle, so to speak, occurs between the 
points A, B, and C. Below the line BC, the strata are less conductive, 
and the same is true of the strata above the water table. Due to the 
absence of water west of the fault in the Cretaceous section, these 
strata also are poor conductors. 

The resistivity profiles /7, and are similar to each other, but 
differ greatly from the profile I. They indicate a high resistance peak, 
corresponding with the Quaternary top layer and the Fernando strata 
above the water table. Then the resistivities drop to very low value, 
corresponding with the depth of the conductive triangle which is in- 
volved in every case. Thus, the horizontal part of the resistivity curve 
indicating the part of this triangle traversed by each traverse, is very 
short in traverse JV, longer in J//, and still longer in traverse JJ. As 
soon as the impervious strata are reached, the resistivity curves go up 
again, this rise occurring at the shallowest depth in traverse JV, being 
deeper at JJ], and the deepest in traverse J]. In other words, the point 
where, on each curve, this rise begins, gives approximately the depth, 
and thus the dip of the impervious formations. The agreement with 
the geologic section goes still further. After the impervious layer has 
been traversed in each section, the recorded resistivities are in pro- 
portion with the resistivities of the section below this impervious 
series. Thus, if we take a certain point on each curve, for example, 
those corresponding with a deflection of 20 millivolts, and connect 
these points, the resulting line gives directly the dip of the strata be- 
low the impervious layers. 

The resistivity curve obtained at the traverse J indicates forma- 
tions of much higher resistivity than the curves obtained on the other 
traverses. The peak corresponding with the dry layer above the 
water table which was obtained on the other traverses is also missing. 
A distinct break occurs in this curve which is explained by the fact 
that the apparent resistivity drops when a depth is reached which is 
equal to the horizontal distance (at that depth) of a better conducting 
medium from the traverse line. We, therefore, interpret this break 
in the traverse J as the effect of the strata east of the fault plane. 

In Figure 35, the results of two potential-drop-ratio traverses 
across an area in the Spring Coulee district in Alberta are shown. 
The geologic section is Cretaceous below glacial drift, the Cretaceous 
being alternating shales and sandstones of the St. Mary and the Fox 
Hills formations, The measurements are made with a Racom equip- 
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ment by the Swedish American Prospecting Corporation. As seen from 
Figure 20, in electrical vertical drilling by the Racom method, forma- 
tion boundaries are represented by peaks in the ratio curves. Plotting 
the data as done in Figure 35, that is, ratios higher than one, on the 
right, and ratios lower than one, on the left, a comparison with Figure 
20 shows that a peak on the right means a poor conductor below that 
peak, and a peak on the left means a good conductor below the peak. 
Thus it is readily seen that the sandstones are poorer conductors than 
the shales, which is in accordance with an observation which was 


A.P Con. Well 


Wig St. Marys 


Fic. 35.—Results of electrical vertical drilling and corresponding geologic section 
at Spring Coulee, Alberta, Canada. Potential-drop-ratio method was used. (From 
pamphlet of Swedish American Prospecting Corporation.) 


previously made when discussing the results represented in Figure 31. 
It is also seen from Figure 35 how well the electrical results check the 
well data and how the electrical vertical drilling by means of the 
potential-drop-ratio method may be used to advantage in carrying 
the structural data obtained from wells into unknown territory. 


3. RESISTIVITY CORRELATION IN WELLS 


In discussing the results which have been obtained on structure 
with the resistivity and potential-drop-ratio methods the experience 
obtained by measuring resistivities in wells should by no means be 
overlooked. In the first place, the results obtained by determining 
the resistivities im situ at the depths reached by the drill, give very 
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important information for the interpretation of resistivity measure- 
ments made at the surface. Secondly, the systematic “electrical cor- 
ing” makes it possible to correlate wells by their “resistivity log” in 
much the same manner as it is possible by the sample log, and gives 
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Fic. 36.—Electrical correlations between two wells, Oklahoma City, 1929. (Waters 
No. 3 and Mackey No. 1.) 


not only a possibility for structural correlations, but at times enables 
one to recognize beds which may otherwise be overlooked. 

Figure 36 represents two resistivity logs obtained with the Schlum- 
berger method of electrical coring in two wells 34 mile apart near 
Oklahoma City. In this figure, the resistivities are plotted in log- 
arithmic scale. It is seen that the two curves are very similar in almost 
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every minute detail. All other formations except the shales exhibit 
rather irregular characteristics which makes it possible to correlate 
definitely the shale horizons in the two wells. Water sands are char- 
acterized by very low resistivity; the resistivity of the water sand at 
the bottom of the holes is almost nil. The Tonkawa limestone shows 
by a pronounced and characteristic resistivity peak. 


II. DIRECT LOCATION OF OIL 


There are two possibilities of locating oil directly or of distinguish- 
ing oil-bearing formations from barren strata: first, by resistivity 
measurements in wells; and second, by resistivity measurement at the 
surface. The former method has proved its merit almost everywhere 
it has been applied, but the second possibility has not yet been estab- 
lished as a commercial method, and only with exceptionally favorable 
conditions have direct indications been obtained from oil-bearing 
formations where they were known to exist at shallow depths. 

The possibility of obtaining a direct indication from an oil-bearing 
horizon in a well does not only mean the possibility of tracing this 
same oil formation in another well; but, as experience has shown, it is 
possible to determine where this oil formation becomes barren or 
changes into salt water in other wells (which is of practical importance 
in deciding the depth of water shut-off), as well as how thick and 
how prolific an oil sand is, the information on the productivity being 
a valuable guide which can not be obtained from the drill record. 


I, RESISTIVITY MEASUREMENTS IN WELLS 


Figure 37 illustrates without further comment how much more 
resistant oil-producing formations are than non-productive forma- 
tions. This figure is an example of a resistivity log taken in the Semi- 
nole area, the two sands being saturated with heavy oil. It is seen that 
the resistivities measured are much higher for oil-bearing beds than 
for highly resistant, compact sedimentary formations. In case of 
doubt the ambiguity can be eliminated by simultaneous records of 
porosities, the technique of which is described in Schlumberger’s 
paper on electrical coring (ref. list I,,). 

Figure 38, representing several electrical logs taken in the Mara- 
caibo field, not only illustrates very well the possibility of resistivity 
correlations of the oil beds directly in adjacent wells, but also shows 
how the resistivity of an oil bed changes if it passes from the pro- 
ductive to the non-productive stage. The upper tar sand and the first 
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productive horizon both are seen to decline rapidly in resistivity 
toward the edge of the productive zone. 
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Fic. 37.—High resistivity of two producing sands, Wewoka field, Seminole area, 
Oklahoma, 1929. (Peney Reed No. 1.) (After Schlumberger.) 
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Fic. 38.—Electrical correlations between wells in the Maracaibo field, Venezuela. 
(After Schlumberger.) 
2. ELECTRICAL VERTICAL DRILLING IN DIRECT PROSPECTING FOR OIL 
If it is possible to recognize oil-bearing beds in wells in such indis- 
putable manner, is it not possible to get the same or similar results by 
measurements from the surface? The results obtained thus far in ap- 
proaching a solution of this problem will now be discussed. 
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Fic. 39 a.—Geological section in Allen County, Kentucky. 
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Fic. 39 6 and c.—Resistivity traverses (partition method) at same locality (after 
Lee and Swartz), showing differences of resistivity between east and west at station 2 
and station 3. 
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The writer knows of only two! instances where it is claimed that 
it was possible to recognize a direct effect of the oil in the resistivity 
data (in this paper we shall not deal with the claims which have been 
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Fic. 40 a and b.—Results of some potential-drop-ratio measurements, and corre- 
sponding geologic section, in Oberg oil field, Germany. (Plotted after data published by 
J. Koenigsberger.) c: Potential-drop ratios as function of conductivity ratio. 


made in that respect for the electromagnetic Elbof method). The first 
of the two instances in question has been published by Lee and Swartz, 
the second by Koenigsberger (ref. list Nos. I1Is7, and IV2). 

1 As a possible third instance may be mentioned the measurements carried out by 
Swartz and others to determine the depth of rock asphalt in very shallow depth (ref. 


list No. IIIz9, p. 114). The observations were made in the same area in which indica- 
tions from oil-bearing strata in 5-800 feet depth, to be described later, were obtained. 
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Koenigsberger has published results obtained along three potential- 
drop-ratio traverses in the Oberg field in Germany, of which only the 
ones measured on the extreme west and the ones measured on the 
extreme east side of the field lend themselves to graphical representa- 
tion. The location of these traverses is shown in Figure 40 a. By 
designating with Koenigsberger the interval between the power 
electrodes by 2a, results are available in this author’s publications 
for intervals a equal to 100, 200, and 300 meters for the southern 
part, and for the intervals of 50, 100, and 200 meters for the northern 
part of the traverse on the west side of the field, while the distances 
used were 50, 100, and 200 meters in both traverses on the east side 
of the oil field. The productive zone is located on the west side of the 
area, and the oil is taken from a depth of approximately 300 meters 
from strata of the lower part of the middle Jurassic, while in the 
central part of this area the productive zone is nearly 500 meters in 
depth. At the locality at which the eastern traverses were made, the 
lower Dogger strata presumably contain salt water instead of oil. 

From the foregoing figures for the distances used, it is seen that 
potential-drop ratios are given only for four (or three, respectively) 
depths in the western and for three depths in either traverse on the 
east. Koenigsberger states in his publication that the depth corre- 
sponding with the potential-drop-ratio indications is h= 14a, and the 
observed potential-drop-ratio values have been plotted in the middle 
part of Figure 40 in accordance with this depth rule. As for normal 
ground the ratio is 3.57, values greater than this ratio have been 
plotted at the right, values smaller than this value were plotted at the 
left. As seen from Figure 40 c, the indications obtained from poor 
conductors are not nearly as pronounced as the indications from good 
conductors below; in order to equalize the indications somewhat, the 
logarithmic scale was employed to draw the curves shown in Figure 
40 b. 

Koenigsberger, in his article, now makes the statement that “‘the 
influence of insulating formations is quite evident from the fact that 
the ratio drops below 3.6 in the traverse on the west side.” 

The writer does not see how the observed values can possibly be 
correlated with the productive area, as the electrode spacing used 
was in all cases too small to penetrate down to the oil-bearing forma- 
tions. On the other hand, it is true that some oil occurs in this field 
above the productive zone in the strata of the lower Dogger (that is, 
in depths less than 300 meters) and it is also possible that the diffusion 
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of gas from the productive layers into the upper strata, which has 
been previously mentioned, has replaced all salt water and has made 
the strata poor conductors. At what depth this zone begins would be 
difficult to state as the curves are not complete enough to draw such 
conclusions. If this assumption is correct, it follows that the electrical 
characteristics of the Wealden can not be much different from those 
of the Dogger, as the curve observed in the northwest profile in the 
east is not much different from the curve observed in the west trav- 
erse. It actually appears from geologic reports that the Wealden is 
also oil bearing in this area and contains a heavy oil which is not 
exploited, and it is further seen from the results that, if such an oil 
horizon exists below the locality of the northwest part of the east 
profile, this horizon is at shallower depth at this point than it is in the 
west. The southwest part of the east traverse does not show an in- 
fluence of poor conductors at the depths reached. 

As a whole, the data published by Koenigsberger appear to be far 
too meager to permit very far-reaching conclusions. They do not seem 
to indicate that it would be possible with his method to differentiate 
a poorly conductive sedimentary stratum devoid of oil, from a 
stratum which is saturated with oil. It is possible that Koenigsberger 
has used greater electrode separations and thus greater depths of 
penetration in addition to those for which results have been published, 
and it is possible that the results obtained at greater depth sub- 
stantiate his statements. However, the published data do not cor- 
roborate them, and the writer doubts very much that, if the illustrated 
results had been obtained in areas in which the presence of oil was not 
known, anybody would have seen anything unusual about them. 

The second instance where a direct location of oil by a resistivity 
method has been described offers much more concrete information. 
The results in question were obtained on a site in Kentucky, the 
geology of which is shown in Figure 39 a. Lee’s partition method was 
used, and the curves obtained when going to either west or east of 
both stations 2 and 3 are illustrated in Figure 39 6. The depth of the 
oil zone from the surface was in this place, however, much less than 
the depth of the oil zone on which Koenigsberger did his work; it 
was only 250-300 feet, while the depth of the productive oil forma- 
tions in Koenigsberger’s experiments was about three times as great. 
The oil zone occurs in the form of an oil sand between wells 2 and 3. 
Actually, east from station 2, the observed resistivity values were 
greater than those on the west, and the converse holds for station 3; 
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toward the east, the resistivities become smaller than the values 
obtained toward the west. 

This is the only place at which, to the writer’s knowledge, it has 
been definitely established that the oil can be located directly by 
surface resistivity measurements. 

After this article was written, Swartz (ref. list No. IIIs7) pub- 
lished the results of further resistivity measurements in Kentucky. 
In the recent observations, the depth of the oil-bearing sands was 
much greater than in the determinations discussed above; it varied 
between 500 and 800 feet. It was found that generally the resistivity 
curves rose where oil and gas were present and fell where dry terri- 
tory was encountered. One casé is cited where predictions as to the oc- 
currence of hydrocarbons made on the basis of resistivity determina- 
tions were checked later by drilling with striking accuracy. Further- 
more, in the area of the Legrande pool, the boundary of the dry and 
producing area could be outlined from the resistivity data and checked 
the information obtained from wells. The striking fact about this 
boundary is that it is not marked by edge water, but by a decrease in 
porosity with corresponding decrease in the oil and gas content: 
another fact worth mentioning is that the best indications were ob- 
tained from a gas sand only 4 feet thick in a depth of about 640 feet; 
that is, the sand detected was only 0.6 per cent of the depth. 

The author concludes that “for the partitioning method the shield- 
ing effect of such low resistivity beds as are encountered in oil and gas 
fields are negligible as far as the detectability of underlying gas and 
oil horizons is concerned,”’ and continues as follows: “It is feasible to 
locate oil and gas zones, as well as salt water zones, directly by re- 
sistivity measurements on the surface of the ground. Whether such 
beds can be detected under all conditions and at all depths is, of 
course, not at present known. So far it is safe to say only that for 
depths up to 1,000 feet and under the physical and stratigraphic con- 
ditions here obtaining, such beds may be determined with excellent 
accuracy.” 


G. SUMMARY AND CONCLUSIONS 


Resistivity and seismic refraction methods have several physical 
characteristics in common. They both involve the advantage of con- 
trolling the depth of penetration and thus a greater definiteness of 
interpretation than obtained in other geophysical methods. On the 
other hand, they are largely subject to the same limitations as other 
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methods, namely, that the depth of geologic bodies must compare 
favorably with their size or thickness. The resistivity methods have 
been applied for some time in oil work, but interest in them has been 
aroused of late due to the perfection of technique which makes it 
possible to obtain a type of indication which may be interpreted 
more readily (the potential-drop-ratio method). 

Although the determination of actual resistivities of samples 
is not a factor as important in resistivity work as the determination 
of related physical properties in other geophysical methods, remark- 
able progress has been made in the design of apparatus which permits 
the determination of resistivities on samples in the field and in the 
laboratory, and on outcrops of geologic formations. 

The requirements on instruments and technique are not exces- 
sive in resistivity and potential-drop-ratio work; for moderate depth 
of penetration, moderate power sources are sufficient. The equip- 
ment for the determination of resistivities and potential-drop ratios 
is readily portable and easily operated. In both resistivity and po- 
tential-drop-ratio work, two distinct fields of application may be 
distinguished; first, the method of resistivity mapping, and second, 
the method of electrical vertical drilling. For the former, the 4- 
terminal Gish-Rooney equipment with one or two fixed electrode 
separations is usually employed, while for the electrical vertical 
drilling any one of the methods in resistivity work, and three poten- 
tial electrodes near one power electrode in potential-drop-ratio work, 
is employed. 

The results obtained thus far indicate that it is possible to locate 
readily geologic structure with the method of equi-resistivity mapping 
at not too great depth, and that detailed information may be obtained 
at any one point by the method of electrical vertical drilling with a 
fairly great degree of accuracy down to a depth of 1,500 feet with the 
present state of the technique, under favorable circumstances. 

While remarkable progress has been made in identifying and cor- 
relating indications due to oil beds in wells, the location of oil directly 
by measurements from the surface is still in its very early stages of 
development, and there is only one case where the direct location of 
oil in very shallow depth from the surface by electrical measurement 
has been definitely established. 

The resistivity methods depend probably more than any other 
geophysical method on a very good correlation of the geophysical 
data with known geological results, and their chief application is 
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given where the problem is to extend the structural information in 
a known area into adjacent territory. 

The chief obstacle which now faces the direct electrical location 
of oil is that, with our present technique, we can not detect a perfect 
insulator any better than a stratum, the conductivity of which is only 
approximately 1o times smaller than that of the adjacent beds. How- 
ever, there is no reason why in the near future the methods can not 
be perfected to such an extent that, when a resistivity indication has 
been definitely established to be due to an oil bed in a certain depth, 
this indication can be traced into adjacent territory with the object 
of determining whether or not the oil is replaced by salt water. It is 
doubtful, however, whether or not the reverse will be accomplished 
so soon, as it is difficult, at the surface, to distinguish a dry resistant 
formation from an oil-bearing bed. In problems of the type just dis- 
cussed, that is, in tracing a formation with the object of determining 
whether its contents change from oil to salt water, the potentialities 
of the resistivity and potential-drop-ratio methods will probably be 
utilized to the greatest advantage when the resistivity indications 
are combined with results obtained from other geophysical methods 
which would give the structural data alone, such as the seismic re- 
flection method. 
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ABSTRACT 

Factors which have adversely influenced opinions and general impressions regard- 
ing the value of geolectric methods as aids in the search for oil are thought by the au- 
thor to be in considerable part of subjective rather than objective origin. However, 
geoelectric methods entered the field relatively late and with inadequate theoretical 
equipment for the best interpretation of results. Furthermore, sources of error which 
doubtless vitiated results have been in some cases overlooked because these did not 
come into serious account in the previous experience of exploration for ores. The prin- 
cipal sources of error for the resistivity method are pointed out. Although electro- 
magnetic and resistivity methods have shown distinct promise in this field, it is, in the 
author’s opinion, not possible with the data now available to make a reliable compari- 
son, in terms of potential results per dollar, between these and gravimetric or seismic 
methods. 


Published pronouncements on the success of geophysical methods 
as aids to the geologist in the search for oil have been numerous and 


various. It is not intended here to add to the number and it would 
be quite impossible to add to the variety. A careful quantitative 
appraisal of any one method would require the consideration of many 
factors and would accordingly call for a statistical analysis of a great 
mass of data. If all the data now resting in the files of various organ- 
izations were pocled, they would probably be adequate for such a 
study, at least in the case of those methods which have been used 
most extensively. Unless and until such a study can be made, a con- 
siderable diversity of opinion, as to the degree of success of a given 
method under the various conditions thus far encountered, may be 
expected to continue. Such studies of individual methods would con- 
stitute an important preparation for the much more difficult task of 
evaluating the relative merits and determining the particular province 
of each of the different methods. And then there would remain for 
consideration the advantage which accrues from the joint use of two 
or more methods. 


1 Manuscript received, July 21, 1932. 


2 Department of Terrestrial Magnetism, Carnegie Institution of Washington. In- 
troduced by F. H. Lahee. 
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Statistical studies of these questions may seem so laborious as to 
be uneconomical. However, the variables which are not under man’s 
control are here so numerous that the statistical along with all other 
means of attack must be used. It is feared that the importance of 
considering geophysical problems and data from a statistical view- 
point is not fully appreciated. 

The physicist coming from his problems in the laboratory where 
accessory factors are largely under his control, on his early encounters 
with geophysical problems is likely to be either over-optimistic or 
pessimistic, according to temperament. Although his technique and 
thorough knowledge of the underlying principles of method and ap- 
paratus are needed, yet he is likely to be handicapped at the outset 
by a lack of that “‘statistical sense” which is so necessary in the suc- 
cessful geophysicist and which is generally acquired only with time 
and experience. Unfortunately, some persons, otherwise well qualified, 
seem unable to acquire this sense or judgment and as a result continue 
for years to harass editors and to clutter up the literature of geophysi- 
cal science with pet theories and general conclusions based on in- 
adequate data or even trivial manifestations. 

One in whom this statistical sense is lacking is likely either to 
extol or to condemn a method after too brief a trial. In either case he 
may be unaware not only of the many phenomena of the earth which 
act as disturbing influences, but also of the fact that many suitably 
distributed data, or samples, are required under such circumstances 
before a reliable interpretation can be made or a decision either for or 
against a method is justifiable. Both experience and temperament go 
to determine the extent to which a judicious balance may be main- 
tained between optimism and pessimism when evaluating methods or 
interpreting results. 

The element of uncertainty in the final interpretation of a geo- 
physical survey must be not only recognized but also admitted by the 
geologist and the geophysicist who share jointly the responsibility 
for the interpretation. The enthusiast can do real harm by leading 
the patron to expect too much of geophysical methods and thus add 
to the damage already wrought by self-styled “geophysicists” of 
either the unscrupulous or the incompetent variety. This probably 
applies with especial force in the case of geoelectric methods. The 
quack and the shyster seem to have a strong predilection for electrical 
vestments. 

Another unfortunate circumstance is that electrical trappings are, 
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in the minds of many laymen, endowed with mystical power. No 
doubt everyone who has worked with geoelectric methods has been 
amused by the curious spectator who loses all interest when he finds 
that the apparatus in use is not provided with a pointer and scale to 
indicate the number of barrels of oil which may be obtained from a 
well sunk at a particular spot. But it is not only the curious spectator 
who entertains this attitude. In recent years some hard-headed 
business men in an eastern oil field patronized a self-styled geophysi- 
cist who with a voltmeter and a few vacuum tubes connected in a 
meaningless circuit professed to perform this miraculous feat. The 
vacuum tube is the superb present-day garnishment for the doodle- 
bug and doubtless whets the appetite of those who have a taste for 
this morsel. 

The bonafide geophysicist will, of course, not make such pre- 
tentious claims. He will recognize in his task a parallel to that of the 
weather-forecaster and will not shrink from using that qualifying 
word, “probable.” He will speak only of indications, good, bad, or 
indifferent. These should be among the passwords, without which 
an entrance into the confidence of prospective patrons can not be 
won. 

Although some geoelectric methods may under favorable condi- 
tions detect oil directly because of its high insulating property, geo- 
electric as well as other methods are efficacious in the location of oil 
only to the extent that they disclose subsurface structural features 
which in turn may indicate oil-bearing structure. That geoelectric 
methods are capable of disclosing hidden structure had been defi- 
nitely shown before these methods were tried in the oil field, and 
since then a number of cases in which geoelectric surveys disclosed 
oil-bearing structure have been reported. Some of these are mere 
statements and have to be accepted on faith, but in several reports 
data are published which seem definitely to support the claims. 

It is, of course, obvious that structure will be revealed by any 
geophysical method only when that property of the structure upon 
which the method depends presents a certain degree of contrast to 
the surroundings. Thus two different methods, as, for example, the 
gravimetric and the electric, may reveal entirely distinct structures, 
or one method may give indications of structure while the other does 
not. It must not be assumed, however, that the method which gives 
no indication is not of positive value, for if the information which it 
supplies is considered along with that obtained by the other method, 
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then more definite conclusions as to the nature of the structure are 
likely to be justified. 

The electrical conductivity, or the resistivity, of earth materials 
is the electrical property upon which those geoelectric methods 
which have shown some promise in the oil field depend. Methods em- 
ploying radio waves would, of course, respond to differences in the 
dielectric constant and some entertain hopes for such methods. Per- 
haps, when the technique for sending directed beams of radio waves 
is sufficiently developed some success will be enjoyed by these 
methods. 

The resistivity of earth materials varies between great extremes. 
Thus of the materials which come into account in the oil field, salt 
water, at one extreme, may have a resistivity as low as 10 ohm- 
centimeters, whereas the value for the igneous bed rocks, the other 
extreme, may run high in the millions of ohm-centimeters. Oil itself 
has very high resistivity, but to utilize this property for direct de- 
tection requires a relation between the extent of the pool and the 
depth from surface which probably occurs only rarely. Although some 
claims of direct detection have been made, the evidence presented is 
not conclusive. 

The diversity in the resistivities of earth materials is far greater 
than that of any other property upon which geophysical methods de- 
pend. This is, however, not the advantage that it may at first thought 
seem. One reason for this is that the condition for detecting structure 
does not continue to improve indefinitely as the contrast in resistivity 
increases, but a limit is soon approached where “enough is as good as 
a feast.” Most of the benefit is obtained when the resistivities differ 
by a factor of roo. Another reason is that this great diversity increases 
the chance that inhomogeneities within a structure and even small 
local differences in constitution may act as disturbing factors. Of 
course, if the data are obtained in such a way as to constitute an 
adequate sampling, then by suitable treatment in which these dis- 
turbing effects are regarded as accidental errors of observation they 
can be largely smoothed out. Such disturbing effects are doubtless 
met with to some extent in all geophysical measurements. Their 
magnitude is a determining factor both as regards the number and 
the distribution of observations required in a survey in order that 
reliable indications of the more general features may be obtained. 
This is also a factor which plays a part in determining the optimum 
sensitivity of the measuring instruments. Thus, if the accidental 


| 
| 
i 
| 

ff 

| 

i 

i 

| 

I 

i 


USE OF GEOLECTRIC METHODS IN SEARCH FOR OIL 1341 


errors of observation are about one-tenth the magnitude of the dis- 
turbing effects, the measurements are doubtless of adequate precision. 
Greater precision may make the work unduly tedious and lead to 
confusion. This does not apply to systematic errors of measurement. 
Advantage may accrue if these are of an order of lower magnitude 
than that of the accidental errors. 

The various electric methods which have been in use for one or 
another phase of geophysical exploration may be classed as follows: 
radio methods, electromagnetic methods, equipotential methods, re- 
sistivity methods, and miscellaneous electric methods. It appears 
that of these the electromagnetic and the resistivity methods have 
enjoyed a measure of success in the exploration for oil. When this is 
taken along with the fact that, in comparison with gravimetric and 
seismic methods, these electric methods are late comers in this field, 
their present prospects seem good. 

The introduction of electrical methods into the more difficult 
phases of exploratory work has doubtless been hindered by the tardy 
development of a quantitative theoretical basis for the interpretation 
of the electrical survey data. In the case of seismic and gravimetric 
methods such a basis had been fairly well developed through their 
years of use in the basic studies of geophysics. The theoretical prob- 
lems presented by the electric methods are more complicated, but 
these have been attacked in recent years by several investigators and 
the results obtained are of considerable practical value. This is es- 
pecially true of those investigations which apply to resistivity methods. 
The detailed solutions here are restricted to the simpler ideal cases, 
but these provide norms with which survey results may be compared. 
It is to be expected that calculations will be extended to include other 
typical cases and that in the near future methods to facilitate such 
calculations will also be developed. 

From the theoretical considerations it appears that geophysical 
measurements alone will admit of a unique quantitative interpreta- 
tion only in special cases. Of the several interpretations which would 
be consistent with the measurements, some may be eliminated as 
being incompatible with facts known to the geologist. The number of 
possible interpretations is likely to be further reduced if surveys by 
more than one method are made. Thus, by the consideration of other 
information along with the survey data, it seems possible to decide 
upon one as the most probable interpretation. When there is added to 
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this the data from a single test hole, the structure of the region covered 
by the survey should be determinable in many cases. 

Interpretation of some aspects of geoelectric surveys has often 
been made with the aid of empirical rules, for which no satisfactory 
theoretical support has yet been found. A few remarks on one of 
these, for which the writer bears some responsibility, may be of inter- 
est. In that form of development of the resistivity method where four 
electrodes are set in line at equal intervals, it was at the outset tenta- 
tively considered that the distance between adjacent electrodes was 
a rough measure of the depth of earth involved in the measurement. 
This rule seemed to do better service than was expected. Such checks 
as it was possible to make strengthened confidence in it, and even to 
the present time some geophysicists with extensive experience in 
geoelectric work continue to use it even though they are aware that 
the theoretical considerations of recent years give no general support 
to such a rule. The observational data now available to most workers 
are probably inadequate to determine whether there is a conflict 
here between fact and theory. The theory seems sound, except for 
the assumption that earth materials are isotropic as regards the 
property of resistivity. There seems to be considerable evidence that 
this assumption is not justified, but that the resistance to electric 
flow in a horizontal direction is different from that in a vertical direc- 
tion. A more general theory, embracing eolotropic materials, is, be- 
cause of its difficulty, not likely to be developed unless overwhelming 
evidence is adduced to show that the simpler theory is inadequate. 
Hence at present it would seem best to use the theory now developed 
as a guide in the interpretation of resistivity surveys. 

One who is confronted with the task of selecting an electric method 
for use in the search for oil will find his chief difficulty in making a 
choice between electromagnetic and resistivity methods. If he has 
had experience with one of these, he will probably do well to select 
that method, for in “the present state of the art” success depends to 
a considerable extent upon familiarity with the method and ap- 
paratus, and this is not gained in a day. The remarks which follow 
may help to elucidate this and, at the same time, call attention to 
some systematic errors which have not always been avoided in re- 
sistivity measurements. 

In 1922 the writer gave consideration to various possibilities of 
obtaining a measure of the resistivity of large masses of undisturbed 
earth. Such data were desired in connection with investigations of the 
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natural electric currents in the earth which were then being initiated 
in the Department of Terrestrial Magnetism of the Carnegie Institu- 
tion of Washington. The method of measurement which was finally 
decided upon should be regarded as a direct-current method.' It was 
thought that with this method it should be easier to carry out sur- 
veys, especially when great depths of earth are to be included in the 
measurements, and that the interpretation should be simpler than 
for a number of other possibilities which were examined. This method 
is an improvement and extension of that developed by McCullom? 
and used by him in electrolysis surveys and which was in turn based 
on an alternating-current method suggested by Wenner.’ The first 
preliminary survey, made in the early autumn of 1924 with the ap- 
paratus designed for large-scale measurements, showed that the 
method could be used for determining and locating structure under 
certain circumstances. Following the publication of a report on this 
work, the method was adopted in some form by several organizations 
which had previously made extensive use of equipotential methods. 
In response to requests in recent years, working drawings and other 
information required for the duplication of the apparatus have been 
supplied by the Department of Terrestrial Magnetism to a number 
of responsible individuals and research organizations. It appears from 
reports, published by some of these, that a few of the more essential 
features of the apparatus have at times been overlooked or mis- 
understood. The fact that this occurred, even where the work was 
in charge of capable physicists, indicates the need of the additional 
emphasis which a restatement of these features may give. 

There are situations where the detection of the desired structure 
is so simple that relatively crude apparatus will suffice. However, 
when the exploration is to extend to greater depths, in the hundreds 
and into the thousands of feet, and under conditions which are made 
difficult by the relationship between the conducting and non-conduct- 
ing strata, success is dependent upon careful attention to the follow- 
ing matters: (1) the polarization, contact potentials, and earth po- 
tentials of various origin of the pick-up electrodes must be largely 
eliminated; (2) the effects of induction between the two circuits must 
be made negligibly small; and (3) the insulation between the rela- 
tively high-potential energizing circuit and the low-potential pick-up 


1 Terr. Mag., 30 (1925), pp. 161-88; Bull. Nat. Res. Council, 11, Pt. 2 (1926), 
pp. 86-91; U. S. Patent No. 1, 813, 845 (1931). 

2 “Measurement of Earth Currents,”’ Elec. Ry. Jour. (November 5, 1921). 

+ Bull. Bur. Standards, 12 (1916), pp. 469-78. 
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circuit must be maintained in a much better condition than is required 
in the more common electrical measurements or else some provision 
must be made in the design to obviate the effects of leakage currents. 

In the apparatus developed at the Department of Terrestrial 
Magnetism by the writer, the contact potentials, etc., are eliminated 
by employing a double commutator which reverses the energizing 
current about twenty times per second and at the same time rectifies 
the pick-up current so that a direct-current potentiometer may be 
used in this part of the measurement. If this is to be satisfactorily 
accomplished, the commutator must be carefully made and properly 
adjusted. When the contact potentials, for example, are very large, 
0.5 volt and upwards, the forced vibration of the pointer of the 
potentiometer-galvanometer becomes so great that settings are diffi- 
cult. This difficulty is eliminated if a condenser of about 20 micro- 
farads capacitance is thrown in series with one of the pick-up lines. 
The condenser should not be required when the contact potentials 
are much less than o.5 volt unless the galvanometer has characteris- 
tics which make it unsuitable for this work or the adjustment or con- 
struction of the commutator is defective. If the latter be the cause, 
serious systematic errors may enter the measurements. 

The effects of induction between the two circuits can be very seri- 
ous, especially when measurements are being made with longlines. 
Some workers have attempted to correct for the effects by obtaining 
measurements at two different and known rates of commutation, but 
in the opinion of the writer this is not satisfactory either in theory or 
in practice. These induction effects may, however, be avoided if the 
potential between the pick-up points is measured only during that 
part of the cycle when the energizing current is constant. That this 
can be accomplished automatically if the commutator is designed and 
adjusted with this end in view has been demonstrated by many com- 
parisons made by W. J. Rooney.' It is, of course, obvious that if the 
commutator in the present form should be used with lines many times 
the length of any employed heretofore in surveys, the time constant 
of the circuit may become so great that adjustment would no longer 
be possible. It should also be mentioned that an improperly adjusted 
commutator can falsify the potential measurements by effectively 
shunting the potentiometer during a part of the cycle. 

Errors from defective insulation may at times be very great. Dur- 
ing early experience in this work, errors of several hundred per cent 


' Terr. Mag., 35 (1930), pp. 61-72. 
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were found on a few occasions. It is when these give rise to negative 
values in the calculated resistance that they are especially conspicu- 
ous. With some types of apparatus which have been used in resistivity 
surveys, such absurd results are not forced upon the attention of the 
observer, but the error is not less serious for his being unaware of it. 

The errors which may arise from defective insulation on the field 
cables are nearly always additive. It is only for some unusual positions 
of the cables that they cause the measured values to be too small. 
However, in all our experience, errors from this source have been 
inappreciable except on a few occasions when old cables were used in 
rainy weather. 

It is in the instruments that troublesome insulation defects may 
develop. The errors which result from these are in part dependent 
upon the relationships between the circuits in the instruments and 
upon the contact resistance of the electrodes. Since both these factors 
may be varied to some extent, this provides a means of detecting such. 
errors. Four different circuit relationships are readily obtained by 
changing the connections of the field cables with the instruments. 
The values of resistivity obtained from measurements with these 
different combinations may all differ if leakage effects are present. 
Their mean is in general not free from leakage error as has been as- 
sumed by some. These errors are completely eliminated from this 
mean only when the contact resistances of the pick-up electrodes or 
those of the energizing electrodes are equal, and this holds only if the 
condition of the insulation remains constant during the series of 
measurements. 

The uncertainty and inconvenience which attend attempts to 
eliminate these errors by such means can, however, be avoided by 
making suitable provisions in the apparatus. The device which has 
been adopted for this purpose is similar in conception to the “guard 
ring” which is often used, especially in instruments employed for 
measuring very small quantities of electricity. This consists essentially 
of an independent system of conductors so arranged as to intercept 
all possible paths by which current may leak from the energizing 
circuit to the pick-up circuit. If this guard system is connected to 
earth at a relatively neutral point, its potential is maintained at a 
value which differs so little from that of the pick-up circuit that this 
serves to reduce such errors to several orders’ lower magnitude than 
those which would result if this device were not employed. In thou- 
sands of tests made by W. J. Rooney under a wide variety of con- 
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ditions, no evidence of errors from this source has ever been found 
when the “guard” was in use. 

The foregoing rather technical remarks apply when a commutator 
is used in the resistivity measurements. Resistivity surveys have, 
however, been carried out without employing a commutator but by 
using, instead, simple hand-operated switches. The errors which may 
enter into measurements made by such a method include those which 
have been discussed for the case where a commutator is used. Of 
these, that type which arises from the various extraneous potentials 
which act upon the pick-up electrodes is likely to be much more 
serious with this method than with the commutator method. This is 
partly because, even with the best technique, the time required for 
a single setting of the potentiometer and the interval between read- 
ings for reversed current are both so long in comparison with the cor- 
responding cycle in the commutator method that the error in the 
result is likely to be considerably greater on this account alone. On 
the other hand, a single reading obtained by the use of a commutator 
is effectively a mean taken over a considerable number of cycles and 
on this account is of enhanced accuracy in cases where errors of a 
random nature may be appreciable. In the simple switching method, 
errors due to defective insulation may be controlled by the design 
and arrangement of the switches, whereas induction effects will de- 
pend upon the technique of measurement. As far as the last two 
sources are concerned, the results obtained with the simple method 
should be better than those obtained with an unsuitable commutator. 
The writer is, however, of the opinion that if the possibilities of the 
direct-current resistivity method are to be fully realized, suitable 
mechanical means of rapid commutation are required. 

The search for oil often requires exploration to greater depths than 
was required in the earlier applications of geoelectric methods, 
namely, in the search for ores. Views as to the capability of geoelectric 
methods for such deep exploration vary, whether these be based upon 
theory or practice. Obviously the details of structure will become less 
distinct the greater the depth from the surface. Thus, a feature of the 
hidden structure can be detected at the surface only when it is of 
sufficient extent. The necessary relation between depth and extent is 
dependent upon the accuracy of the measurements and upon the mag- 
nitude of the uneliminated influences of the disturbing features of the 
earth. The accuracy of the measurements can be controlled by suit- 
able design and operation of the measuring apparatus. The errors 
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arising from disturbing features of the earth, such as topography and 
inhomogeneities within structure, are to a considerable extent de- 
pendent upon the number and the distribution of measurements and 
upon the method by which these data are treated. In order that cor- 
rections may be made for the more general features of topography, a 
further development of theory is required. 

Although there is some evidence that features of structure located 
between 2,000 and 4,000 feet below the surface have been disclosed 
by resistivity surveys, this is by no means definitely established. 
General theory, as the writer understands it, does not exclude such 
a possibility. In fact, the likelihood of exploring to far greater depths 
seems so good to some of us that plans are being developed for a large- 
scale codperative undertaking in which it is hoped to obtain by re- 
sistivity measurements information about the structure of the earth’s 
crust down to a depth of about 50 miles. The results to be obtained 
from such a venture are, of course, not expected to be of practical 
value, except perhaps indirectly, but they should have important 
bearing on problems in several branches of the basic science. 

The principal points which it has been attempted to bring out in 
this article may be briefly stated as follows. Although evidence is at 
hand showing that some geoelectric methods are aids in the search 
for oil, yet it can not be determined at this time whether the results 
per dollar from these methods compare favorably with those obtained 
from gravimetric or seismic methods. 

An attitude of disfavor toward geoelectric methods has perhaps 
in considerable part sprung from subjective rather than objective 
sources. The widespread belief in the near-omnipotence of electricity 
has brought forth many impotent geoelectrical schemes and has pro- 
vided easy prey for these as well as for the outright faker. Of course, 
an unfavorable reaction followed. Furthermore, a bonafide geo- 
electric method can not qualify as the idol of these “believers.” 

Geoelectric methods also entered this field with the handicap of 
inadequate theoretical equipment and with practical experience in 
a class of problems not adapted to bring out some details of apparatus 
and method important in the more difficult work to be encountered 
here. These conditions doubtless added considerably to the cost and 
lessened the value of the results obtained from some of the geoelectric 
surveys. And withal there is the seniority of gravimetric and seismic 
methods to be considered. 


In principle it seems entirely feasible in many cases to determine 
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from the data of resistivity surveys, taken along with other informa- 
tion which can be supplied by the geologist and the petroleum engi- 
neer, the approximate depth and the features of the more general 
structure associated with petroleum resources. 

Reports of actual tests seem to bear this out, perhaps rather too 
well in some cases. This applies especially to some recent tests which 
came to attention!’ as this conclusion was being written. These con- 
sisted of the correlation of hundreds of drill-hole records with results 
of resistivity surveys made by five different parties in as many differ- 
ent oil fields in Russia. According to the abstracts, one of the articles 
contains the conclusion that “‘more tests are required to draw defi- 
nite conclusions on the degree of usefulness” of this electric method. 
This leaves one uninformed as to whether or not any indication that 
the method may prove useful was found. However, there is no am- 
biguity in the four other conclusions, which run as follows: ‘‘This 
method is of great importance in geological interpretation of the sub- 
soil and should be used on a large scale’”’; ‘‘whether the mechanical 
coring can entirely be substituted by this method or not can not be 
definitely decided at present.’’ One in which 180 drill-holes were in- 
vestigated concludes that this ‘must be recognized as one of the most 
useful and necessary methods of determining the correlations in the 
stratigraphical conditions of oil deposits.” In what appears to be the 
most extensive of the five investigations 220 drill holes were covered. 
In the report of this it is concluded that “the great advantage of 
electrical coring is established.” The full reports may restrict some of 
these conclusions more than is indicated in the abstracts. 

It is perhaps utopian to wish that these investigations had also 
included surveys by the other outstanding geophysical methods. An 
investigation of such scope should have considerabiy lessened the 
element of guesswork which is involved in making a choice from the 
several available methods. The writer is not impelled to venture such 
a guess here; he is, however, convinced that some ge electric methods 
will in due time be accepted as useful aids in the search for oil. 


! Geophys. Absir., No. 37 (1932), pp. 450-53. 
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MINERAL CONTENT OF SOILS' 
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ABSTRACT 


The presence of radon in the soil cannot be used as a criterion for locating oil 
fields, since the source of the radon lies in the radioactive minerals of the soil and 
not in the petroleum at greater depths. Formation contacts could be located by a 
radon survey even though obscured by top soil and vegetation, because there is 
quite likely to be a difference in heavy mineral character and content of two different 
formations. Likewise faults may be indicated if they bring two quite different forma- 
tions into juxtaposition, but there is always danger of misinterpreting variations in 
the radon content of a soil which may be due simply to local variations in the amount 
of radioactive minerals in one and the same formation. 


INTRODUCTION 


In recent years radioactive methods of geophysical prospecting 
have attracted some attention. Most prominent among these methods 
has been that of measuring the radon content of soils‘and attempting 


to deduce from these results the nature of the geological substrata. 
Radon is one of the constituents of the soil gas and is the immediate 
decay product of the metallic element radium. The radon itself is 
radioactive and its radioactivity per unit of weight is several times 
that of radium, thus enabling extremely small quantities of the gas 
to be detected. It is claimed that a fault will be evidenced by a higher 
radon concentration in the overlying soil,5 and that the radon content 
of the soil above an oil bearing formation is higher than elsewhere.*® 
Since there were no data available as to the success or failure of these 
radioactive methods in this country, an extensive study was made of 
the radon content of soils in regions whose geology and lithology were 


1 Manuscript received, June 23, 1932. 
2 Geologist, Gulf Companies. 
* Physicist, Gulf Research Laboratory. 


4 “Radon Content of Soil Gas,” H. G. Botset and Paul Weaver, Physics, Vol. 2 
(1932), pp- 376-85. 


F, Miiller, “Radioaktivitalsmessungen als geophysikalische Aufschlussmethode,” 
Zeits. fiir Geophysik, Vol. 3, No. 7 (1927), pp. 330-36. 


®°L. N. Bogoyavlensky, “Radiometric Exploration of Oil Deposits,” Bull. Inst. 
Practical Geophysics (Leningrad), No. 3 (1927), pp. 113-24. 
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quite well known. The authors acknowledge their indebtedness to 
A. E. Ruark, of the University of Pittsburgh, and Paul Weaver, 
Houston, Texas, for valuable suggestions inthe conduct of this work. 


APPARATUS AND FIELD METHODS 


The measurements were made in soils which were developed over 
sedimentary formations consisting predominantly of sandy clays with 


SIO 


Fic. 1.—Diagrammatic section of radon measuring apparatus. 


some slightly consolidated sands, all of Tertiary age. Most of the 
measurements were made in wooded areas, but some were in culti- 
vated fields, and some in natural prairies. 

The technique of making the measurements is a modification of 
the method developed by Ambronn.' The equipment used is shown 
diagrammatically in Figure 1. The hole into which the sampling tube 
is placed is made by first driving a steel rod into the ground about a 
foot and then continuing the hole with an earth auger. This is to 


1 R, Ambronn, Elements of Geophysics, English translation by Cobb (1928), p. 119. 
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avoid compacting the soil around the hole and thus reducing its 
porosity, which would tend, in clayey soils, to make it more difficult 
to obtain a sample of soil gas. The sampling tube (A) consists of a 
14-inch iron pipe, broken and coupled near the top to permit the 
insertion of nipples so that the tube can be extended when samples of 
soil gas from greater depths are desired. A wall packer (B) is welded 
on the outside of this tube about a foot from the bottom. The packer 
consists of a piece of tubing of the same diameter as the hole, around 
which an iron rod (C) of 14-inch square cross section was wrapped 
spirally and welded in place. This rod cuts into the walls of the hole 
and forms a tight seal against leakage from the atmosphere. The gas 
sample is pumped through drying tubes (D) by the piston pump (£) 
into the ionization chamber (F) of the electroscope. The ioniza- 
tion chamber is a cylinder whose bottom is a movable piston (G). To 
empty the ionization chamber of gas the outlet cock on the top is 
opened and the piston pushed up. The outlet cock is then closed and 
the inlet opened. The gas sample being pumped in thus exerts an 
outward pressure on the system and any leakage occurring is out- 
ward. This prevents the dilution of the gas sample with atmospheric 
air. The minute traces of radon in the soil gas, by virtue of radio- 
active decay, produce an ionization of the gas in the ionization cham- 
ber and discharge the electrometer fiber which has been previously 
charged to a potential of about 200 volts. The rate of discharge is 
proportional to the amount of radon present in the gas. The sensi- 
tivity of the instrument is such that one part of radon can be de- 
tected in about 10” parts of soil gas. 

Several profiles were run across a known fault in the Balcones fault 
region. The measurements made in these profiles are shown in Figure 
2. It is difficult to express the results in terms of quantity of radon per 
cubic centimeter of soil or per gram of soil mineral, because of the 
great variation in the porosity of the soil and because of another factor 
which may be even more variable, the so-called emanating power of 
the soil grains. Since the radon gas is the product of the decay of 
radium atoms, the radon must be formed on the surface of the soil 
grains or within the grains themselves. As the radon is a gas formed 
from a solid, only that part formed at the surface of the grains is able 
to get into the soil gas. The emanating power is the ratio of the 
amount of radon which gets into the soil gas to the total formed in 
the soil. This obviously depends, then, not only upon the size of the 
soil grains, but upon their shape (ratio of volume to surface) and upon 
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whether the parent radium has been deposited near the surface of 
the grain (secondary deposition), or whether it is distributed more 
or less uniformly throughout the grain. Naturally radium formed from 
uranium salts deposited chemically on the surface of the grains will 
have a large emanating power, while that formed in minerals in which 
the uranium was originally incorporated will have a lower emanating 
power. According to Satterly' the emanating power of a soil may vary 
from one-sixth to one-twentieth. Consequently, the ordinates of the 
curves presented here are expressed simply as electroscope readings, 
no account being taken of the variations in emanating power of differ- 
ent soils. These ordinates are divisions per second of movement of 
the electrometer fiber multiplied by 1,000. A calibration measure- 
ment on the electroscope made by using a known quantity of radon 
showed that an electroscope reading of 100, that is an actual rate of 
movement of the electrometer fiber of 0.1 division per second, was 
produced by 1,539X10~” curies of radon in the ionization chamber. 
The curie is the amount of radon in equilibrium with one gram of 
radium. If we assume a soil porosity of 33 per cent, 1,539 X 10-” curies 
in the electroscope would correspond to about 0.57107” curies of 
radon per cubic centimeter of soil. 

As will be seen from Figure 2, an attempt to locate the fault by 
radon measurements would be very misleading, since the radon values 
obtained at the points marked B, where there is no fault, are just as 
high as, or higher than, those obtained at A, where the fault is actu- 
ally located. 

In order further to study the radon content of soils, surveys were 
made over the South Liberty salt dome. This dome has at present 
very little elevation above the surrounding terrain. The center of the 
dome is crossed by a small river. The eastern edge of the dome is a 
producing oil field. It is claimed by Bogoyavlensky’ that an oil field 
will be indicated by a change in the radioactivity of the soil. To test 
this statement a profile was run across the producing section of the 
dome, starting from the outside edge of production and extending 
toward the center of the dome some distance beyond the producing 
area. The result of this survey is shown in Figure 3. The holes from 
which samples were taken for these measurements were 5 feet deep. 
The higher radon concentration in the producing area is obvious. 
After these results were obtained the profile was extended eastward 

1 J. Satterly, Proc. Cambridge Society, Vol. 16 (1911-12), pp. 336, 356, 514. 

2 Op. cit. 
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off the dome for at least 4 mile, and the radon content was still 
practically as high as was found in the producing area. This seemed 
to indicate that the presence of oil had little to do with the radio- 
activity of the soil. 

Knowing the association of radioactivity with heavy minerals,’ it 
was decided to make another profile parallel to and near the first and 
to take samples of the soil from the bottom of each hole. These 
samples were taken to the laboratory for quantitative determination 
of their heavy mineral content. 
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Fic. 3.—Radon profile across South Liberty salt dome. 


HEAVY MINERAL SEPARATION 


About 30 grams of each of these samples were ground to pass a 
130-mesh sieve. The heavy minerals were first separated by using 
bromoform in a separatory funnel. This bromoform had a density, as 
determined by the pyknometer, of 2.5. It was found, however, that 
the heavy mineral fraction obtained by this separation contained 
quartz and calcite. Consequently, this heavy mineral fraction was 
again separated, Thoulet’s solution (density 3.15) being used. This 
gave a very complete separation and was adopted as standard pro- 
cedure. 


1 E. H. Biichner, Jahrb. f. Radioakt., Vol. 10 (1913), p. 516. 
W. Waters, Phil. Mag., Vol. 19 (1910), p. 903. 
R. H. Strutt, Proc. Roy. Soc., Vol. 84 (1910), p. 377- 
A. Gockel, Die Radioaktivitat von Boden und Quellen (Braunschweig, 1914). 
St{Meyer and Schweidler, Radioaktivitét (Leipsic, 1927). 
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The following table taken from Gockel will indicate some of the 
heavy minerals associated with radioactivity. 


TABLE I 

Radium per Gram of 

Mineral Location Mineral X 10 Grams 
Zircon Kimberley 19.1 
Eudalite Greenland 12.6 
Orthite Sweden 236.0 
Gadolinite Hitteré 156.0 
Keilhanite Alve, Norway 452.0 
Niobite Connecticut 97.0 
Apatite Canada 14.6 
Cerite Sweden 30.0 


The average radium content of sedimentary rocks is about 1.4X10~” 
grams! per gram of rock. Thus it is seen how considerably the radio- 
activity is concentrated in the heavy mineral constituents. 

Petrographic examination of the heavy mineral residues from the 
South Liberty dome revealed the presence of zircon, rutile, brookite, 
tourmaline, magnetite, ilmenite, apatite, biotite, augite, and cyanite. 
Zircon, tourmaline, and magnetite are the most abundant in all the 
samples examined. If the radioactivity of the soil could be ascribed 
to the mineral composition of the soil itself rather than to any deep- 
seated source, zircon is the one mineral of all those observed that is 
the most likely to be radioactive. F. L. Hess* describes a considerable 
number of unusual rare-earth minerals occurring in some of the 
igneous rocks of the central mineral region of Texas. Since this area 
was probably the source of much of the material of the late Gulf Coast 
Tertiary, these minerals were undoubtedly present to some extent 
in the soils but not recognized in the heavy mineral residues examined. 
In these minerals are such rare-earth elements as thorium and ura- 
nium. Thorium is the parent of one series of radioactive substances and 
uranium is the parent of the radium series and probably also of the 
actinium series. Since these elements (uranium and thorium) are 
isomorphous with zirconium, it is only logical to look for a possible 
relationship between the radon content of the soil and the amount 
of zircon and other rare-earth minerals in the soil. Experimentally, 
only the total heavy mineral content of the soil was determined, since 
zircon was always present in large percentages in the residues and not 
all the possible sources of radioactivity could be differentiated and 
segregated. 

1 G. Kirsch, Geologie und Radioaktivitét (Berlin, 1928), p. 46. 


2 “Minerals of the Rare-Earth Metals at Baringer Hill, Llano County, Texas.” 
U. S. Geol. Survey Bull. 340 (1907), pp. 286-94. 
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Figure 3 shows almost a parallelism between the radon content of 
the soil and its heavy mineral content. This is one of many profiles 
across the South Liberty dome, all of which are generally similar. 
Such an abundance of corroborative evidence leads one to conclude 
that the radon content of soils is directly related to and due to the 
presence of radioactive minerals in the soils. 


CONCLUSIONS 


The presence of radon in the soil cannot be used as a criterion for 
locating oil fields, since the source of the radon lies in the radioactive 
minerals of the soil and not in the petroleum at greater depths. For- 
mation contacts could be located by a radon survey even though ob- 
scured by top soil and vegetation, because there is very likely to be 
a difference in heavy mineral character and content of two different 
formations. Likewise faults may be indicated if they bring two dif- 
ferent formations into juxtaposition, but there is always danger of 
misinterpreting variations in the radon content of a soil which may 
be due simply to local variations in the amount of radio-active 
minerals in one and the same formation. 
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South Louisiana Oil Scouts Association and Texas Gulf Coast Oil Scouts Asso- 
ciation. Oil and Sulphur Development in the Texas and Louisiana Gulf 
Coast Salt Dome Region. Review by Donald C. Barton ‘ 

South Permian Basin, West Texas, Structural and Stratigraphic Dev elop- 
ment of. By R. L. Cannon and Joe Cannon . 

Southward Extension of Bonham Clay, Texas. Geological Note by C. I. Alex- 
ander and J. P. Smith 

Southwestern United States, Decline of Great Basin. By J.E ‘dmund E ‘aton 

Spofiord, H. N. Pecan Gap Chalk: New Localities in Red River and Bowie 
Counties, Texas. Geological Note. 

Spooner, W. C. Salt in Smackover Field, Union County, Texas. Geological 
Note. 

Stamey, R. A. , Judson, Sidney A. , and Murphy, P. C. Ov erhanging Cap Rock 
and Salt at Barbers Hill, Cc hambers County, Texas : 

Stein, Paul. Leitfaden der Tiefbohrtechnik (Guide to Deep- Drilling Technique). 
Review by R. D. Ree ‘ 

Stewart, Irvine E. Memorial of Howard Walde Kitson 

Stille, Hans. Asymmetric Folds with Reference to German Salt Bodies , 

Stratigraphic and Structural Development of South Permian Basin, West 
Texas. By R. I.. Cannon and Joe Cannon. 

Stream Sculpture on the Atlantic Slope, a Study in the Evolution of Appala- 
chian Rivers. By Douglas Johnson. Review by Sidney Powers. 

Structural and Commercial Oil and Gas Possibilities of Central Valley Region, 
California. By Walter Stalder 

Structural and Stratigraphic Development of South Permian Basin, West Tex- 
as. By R. L. Cannon and Joe Cannon 

Structural Maps, Contoured from Dip Observ ations, Simple Graphical Method 
for Determining True Dip from Two Components and for Constructing. 
By John L. Rich 

Structure, Geologic, Limitations of Ground Water as Aid in Determination of 

Hidden. By E.K.Soper_ . 

Structure Contour Map of the Montana Plains. By C. E. Dobbin and C. E. 
Erdman. Review by MaxW.Ball . 

Studies of Geophysical Methods, 1928 and 1929. By i Gilchrist, 5. B. Mawd- 
sley, and others. Review by Donald C. Barton ‘ 

Stutzer, O. Erdél. Allgemeine Erdélgeologie und Uberblick iiber die Geologie der 
Erdélfelder Europas (Petroleum. General Petroleum Geology and a Short 
Account of the Geology of nena Oil Fields). Review by Parker D. 
Trask . 

Subterranean W ater, Factors Involved in Segregation of Oiland Gas from. By 
JanVersluys. . 

Sulphur and Oil Development in the Texas and Louisiana Gulf Coast Salt 
Dome Region. By Texas Gulf Coast Oil Scouts Association and South 
Louisiana Oil Scouts Association. Review by Donald C. Barton . 

Supplementary Membership List. The Association Round Table. 

Surface Layer, Seismic Weathered or Aerated. By O.C.Lester, Jr... 

Symposium on Occurrence of Petroleum in Igneous and Metamorphic Rocks 
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———. Editor’s Foreword. By Sidney Powers 
on Reservoir Conditions in Oil and Gas Pools 


Taranaki Oil Field, New Zealand. By Frederick G. Clapp 
Teas, L. P. Cameron Meadows and Iowa, Two New Coastal Louisiana Fields. 
Geological Note . 
Technique and Application of Resistivity and Potential- Drop- Ratio Methods 
in Oil Prospecting, Advances in. By C. A. Heiland . 
Texas, Boggy Creek Salt Dome, Anderson and Cherokee Counties. By H. J. 
McLellan, E. A. W endiandt, and E. A. Murchison : : 
—- , Bryson Oil Field, Jack County. By J. P. Bowen and James F. ‘Gibbs 
, Disseminated Oil in Pleistocene Water Sands of Corpus Christi Area. 
By W.Armstrong Price. 
———., Earth Temperatures of North-Central. By Virgil E. “Barnes 
, East, Recovery Problem, Water Encroachment in Bartlesville Sand 
Pools of Northeastern Oklahoma, and Its Bearing on. By D. R. Snow. 
, Gideon Oil Well No. 3, Luling Field, Caldwell County. Geological 
Note by E. W. Brucks. i 
—, High Island Dome, Galveston ‘County. "Geological Note by M. T. 
Halbouty 
— , Interpretation of Grain of. Geological Note by R.E. Rettger 
, Occurrence of Strata of Bend Age in Sierra Diablo, Culberson County. 
Geological Note by M. B. Arick s 
, Oil Fields in Igneous Rocks in Coastal Plain of. By E. H. Sellards | 
, Oil-Field Waters of North-Central. By VirgilE.Barnes. 
’ Overhanging Cap Rock and Salt at Barbers Hill, Chambers County. 
By Sidney A. Judson, P. C. Murphy, and R. A. Stamey 
—, Pecan Gap Chalk: New Localities in Red River and Bowie Counties. 
Geological Note by H. N. Spofford . 
—, Possible Silurian and Devonian Strata in Van Horn Region. Geological 
Note by Philip B. King 
, Regional Structure of C retaceous on Edwards Plateau of Southwest. 
By Lon D. Cartwright, Jr. . ‘ 
—, Shallow Salt-Type Structure in Permian of North-Central. By Urban 
b, ’ Hughes 
——, South, Frio Clay. Discussion by Thomas L. Bailey . 
South, Frio Clay. Discussionby F.H.Lahee_ 
Southward Extension of Bonham Clay. Geological Note by 
_ Alexander and J. P.Smith 
, West, Structural and Stratigraphic Devahement of South Permian 
Basin. By R. L. Cannon and Joe Cannon. R 
Texas and Louisiana Gulf Coast Salt Dome Region, Oil and Sulphur Dev elop- 
ment in the. By Texas Gulf Coast Oil Scouts Association and South 
Louisiana Oil Scouts Association. Review by Donald C. Barton . 
Texas District, West, General Report and Review to December 31, 1931. Com- 
piled by oil company scouts. Review by Sidney Powers ; 
Texas Oil Field, Interpretation of Bottom-Hole Pressures in East. ‘By E. v. 
Foran 
Texas Gulf Coast Oil Scouts Association and South Louisiana Oil Scouts Asso- 
ciation. Oil and Sulphur Development in the Texas and Louisiana Gulf 
Coast Salt Dome Region. Review by Donald C. Barton 
Textbook of Geology, Part I, Physical Geology. By yond R. Longwell 
Adolph Knopf, and Richard F. Flint. Review by K.C.H 
Thomas, Harold S. Proration at Oklahoma City Field, Culshene 
Thomas, W. A., and Fitzgerald, P. E. Occurrence of Fluorite i in Monroe Forma- 
tion of Vernon Township Pool near Mount Pleasant, Michigan. conan 
Note. 
Tiefbohrtechnik, Leitfaden der (Guide to Deep. Drilling Technique). ad “Paul 
Stein. Review by R.D.Reed . 
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Time-Equivalent Versus Lithologic Extension of Formations. Discussion by F. 
A. Meltonand J.Edmund Eaton 

Tioga Region, Pennsylvania, Gas in the. By G.H. Ashley and S. H. Cathcart. 
Review by Sidney Powers . ; 

Torsion Baiance, Accuracy of Determination of Relative Gravity by. By 


Donald C. Barton 

Trask, Parker D. Review of Erdél. Allgemeine Erdélgeologie und Uberblick iiber 
die Geologie der Erdélfelder Europas (Petroleum. General Petroleum 
Geology and a Short Account of the shasta! of ee Oil —— by 
O. Stutzer 


Udden, Johan August. Memorial by E. H. Sellards 

United States of America and Argentina and Oil Prospects in Australia, Report 
on Tour of Inspection of the Oil Fields of the. By W. G. W oolnough. Re- 
view by Donald C. Barton 

United States, Permo-Carboniferous Orogeny in South-Central. Discussion by 
W.A. J. M. van Waterschoot van der Gracht 

———., Southwestern, Decline of Great Basin. By J. Edmund Eaton | 

United’ States and Possessions, Comment on Petroleum of the. Discussion by 
H.B.Goodrich . 

——, Petroleum in the. By Ralph ‘Arnold and W ‘liam J. Kemnitzer. Review 
by 'Sidney Powers 

Upp, J. E., and Condra, G. E. Correlation of the Big Blue Series in Nebraska. 
Review by Robert H. Dott. 

Use of Geoelectric Methods in Search for Oil. By O.H. Gish | 

of Record Character in Interpreting Results and Its E fect on Depth 

Calculation in Refraction Work. By T.L. Allen. 

Uwatoko, Kunio. Genesis of Oil by High Radial Axial Pressure 


van der Gracht, W. A. J. M. van Waterschoot. Occurrence and Production of 

Petroleum in Germany P 
-. Permo-Carboniferous Orogeny ii in South-Central United States. Dis- 
cussion . 

———. The Pre-Carboniferous Exotic Boulders in the So-Called “Caney 
Shale” in the Northwestern Front of the Ouachita Mountainsin Oklahoma. 
Review by Robert H. Dott. 

Van Horn Region, Texas, Possible Silurian and Devonian Strata in. Geological 
Note by Philip B. King ‘ 

Van Tuyl, F. M., and McLaren, R. L. Occurrence of Oil in Cry stalline Rocks i in 
Colorado 

Vers la recherche de nouveaux ix champs pétroliféres (In Search of New Oil Fields). 
By Stanislas Zuber. Review by Charles Ryniker 

Versluys, Jan. I’actors Involved in Segregation of Oil and Gas from Subter- 
ranean Water 

Vermilion Bay Salt Dome, Louisiana, Seismological Discovery and Partial 
Detail of. By E. E. Rosaire and O. C. Lester, Jr. : 

Vernon Township Pool near Mount Pleasant, Michigan, Occurrence of Fluorite 
in Monroe Formation of. Geological Note by P. E. Fitzgerald and W. 
Thomas. 

Vicksburg Formation in Deep Test, Acadia Parish, Louisiana. Geological Note 
by M. T. Halbouty 

Virginia, Possible Origin of Graphit te in Some Ancient Quartzites, Slates, and 
Schists in. By Justus H. Cline ; 

Volcanic Plugs in Mendoza Province, Argentina, Oil Seepages and Oil Produc- 
tion Associated with. By F. H. Lahee ; 

von Zwerger, Rudolf F. Review of Das Erdil in : Deutschland (Petroleum j in 
Germany), by Fred S. Baumann . 


Wahlstrom, Edwin A., and DeFord, Ronald K. Hobbs Field, Lea eaten New 
Mexico 
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Washburne, Chester W. Petroleum Geology of the State of Sao Paulo, Brazil. 

Review by Helen M. Martin 

Water Encroachment in Bartlesville Sand Pools of Northeastern Oklahoma and 
Its Bearing on East Texas Recovery Problem. By D.R.Snow . 

Water Sands, Pleistocene, of Corpus Christi Area, Texas, Disseminated Oil in. 
By W. Armstrong Price. oe 

Weatherby, B. B. Review of Physics, by American Ph ysical Society 

Wegemann, Carroll H. Memorial of Salvador Ortiz Davila . 

Wendlandt, E. A., McLellan, H. J., and Murchison, E. A. Boggy Creek Salt 
Dome, Anderson and Cherokee Counties, Texas . 

Whiteside, Robert M. Geologic Interpretations from Rotary Well Cuttings 

Wilson, J. H., Lovering, T. S., Aurand, H. A., and Lavington, C. S. Fox Hills 
Formation, Northeastern Colorado. Geological Note 

Woolnough, W. 6. Report on Tour of Inspection of the Oil Fields of the United 
States of America and Argentina and Oil Prospects in Australia. Review by 
Donald C. Barton 

Wright, Frank J. The Older Appalachians of the South. Review by Sidney 
Powers . ‘ 

Wyoming and Montana, Character of Producing Sandstones and Limestones of. 


Year Book, 1932, National Oil Scouts Association of America, Inc. Review by 


Zuber, Stanislas. Vers la recherche de nouveaux ures siesta Review 7 
Charles Ryniker . 
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